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ABSTRACT 


This report presents the results of several 
analytical studies related to Space Shuttle naviga- 
tion. It is divided into two Volumes, the first 
dealing with studies related to the addition of 
NAVSTAR Global Positioning System user equipment to 
the Shuttle avionics suite, and the second dealing 
with studies of the baseline avionics suite without 
GPS. The GPS studies center about navigation accu- 
racy covariance analyses for both developmental and 
operational phases of GPS, as well as for various 
Orbiter mission phases. The baseline navigation 
system studies include a covariance analysis of the 
Inertial Measurement Unit calibration and alignment 
procedures, postflight IMU error recovery for the 
Approach and Landing Phases, on-orbit calibration of 
IMU instrument biases, and a covariance analysis of 
entry and prelaunch navigation system performance. 
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1 . 


INTRODUCTION 


For the past six years TASC has participated in the 
development and evaluation of the navigation software for the 
Space Shuttle. The earliest studies, summarized in Refs. 1 
and 2, involved development of error models for a variety of 
navigation aids and generation of detailed performance projec- 
tions for candidate navigation systems for the entry and landing 
mission phases. Studies that followed, summarized in Refs. 3 
and 4, expanded on the earlier studies to encompass navigation 
performance analyses for hangar and preflight phases (IMU cali- 
bration and alignment), for the ascent phase, and for the land- 
ing phases (from the termination of entry through touchdown). 

In the current study, TASC has revised its assessment of pre- 
flight IMU calibration and alignment to keep pace with changes 
in the software procedures, and has expanded the scope of its 
analytic support to encompass baseline navigation performance 
analysis for 

• Post flight IMU error recovery for ALT 

• On orbit calibration of IMU accelerometer 
and gyroscope biases 

• Entry and preland baseline navigation. 

As a major addition to these baseline navigation system 
studies, much of the current effort has focused on consideration 
of the NAVSTAR Global Positioning System (GPS) as a Space Shuttl 
navigation aid. This work grew out of the ascent studies re- 
ported in Ref. 4, and now includes on-orbit and deorbit/entry 
mission phases as well. The principal objective of the GPS 
study effort was to support the Mission Planning and Analysis 
Department (MPAD) in its work with the Space Shuttle/GPS Panel. 
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This involved numerous Panel and working group meetings to 
prepare for the generation of procurement specifications for 
GPS user equipment hardware and software and to prepare in* 
terface specifications between the GPS user equipment and the 
General Purpose Computer (GPC) of the Space Shuttle. One out- 
put of this effort was a series of working papers and recom- 
mendations of an informal nature impacting the Panel studies. 
Another output was the formal, connected performance studies 
presented in Volume 1 of this report. 


1 . 1 OVERVIEW 

The various studies contained in this report all re- 
late to Space Shuttle navigation during the several Orbiter 
mission phases, either with the equipment of the baseline nav- 
igation system, or with the addition of GPS user equipment to 
the Orbiter. Figure 1.1-1 summarizes these studies. Most of 
them were previously reported in Quarterly Progress Reports 
(Refs. 5 to 16), but they have, in some cases, been updated, 
reworked, and combined for this final report. 

1.1.1 GPS as a Space Shuttle Navigation Aid 

While it is not part of the baseline navigation sys- 
tem, GPS user equipment (consisting of antennas, preamplifiers, 
cabling, receivers, and associated interface equipment) is 
currently being planned for addition to the Orbiter avionics 
suite. GPS will have a major impact on Shuttle navigation 
performance. Once GPS itself is fully operational (Phase 111) 
the passive system aboard the Orbiter will provide a semi- 
autonomous source of highly accurate navigation data at almost 
all times, depending on acceptable antenna coverage. The only 
outages in coverage will come during reentry blackout and during 
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OMIT Al COAST 


- navigation with ert 

• SNBU SATtUITC NAVIGATION 

- ON ONMT IMU CAUMATION 




Figure 1.1-1 Summary of TASC Studies Relative to 

Space Shuttle Navigation (1976-1979) 

ascent, where the external fuel tank may block antenna coverage. 
During the GPS Phases I and II, satellite coverage will be 
quite limited, but will still support accurate Orbiter naviga- 
tion on orbit and at times of good satellite visibility during 
deorbit . 


The GPS user equipment differs from the baseline nav- 
igation equipment in that it will incorporate an extensive 
amount of internal software in the receiver/processor assembly 
(R/PA) and will interface on a computer- to-computer basis with 
the Shuttle General Purpose Computers. A multiplicity of is- 
sues related to software and to navigation system performance 
need resolution as part of the planning process to incorporate 
GPS user equipment aboard the Orbiter. Many of these were 
addressed during the course of this study. 
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Since the initial utilization of GPS as a Shuttle 
navigation aid is expected to occur before GPS is fully opera* 
tional , an extensive study of navigation during GPS Phases 1 
and II was performed. That study features signal coverage 
analysis as well as covariance analysis of integrated naviga- 
tion system performance. To look at a limiting case, a situ- 
ation was examined with only one GPS satellite assumed to be 
available. The on-orbit navigation performance of the Shuttle 
was estimated for this condition. 

Since navigation performance with Phase 111 GPS is 
expected to be very good at almost all times, the study here 
focuses on the exceptional circumstances associated with re- 
entry blackout, when the GPS signals will not be available. 

The other situation with potentially limited GPS signal avail- 
ability during GPS Phase III is Shuttle ascent, when the ex- 
ternal fuel tank may block receiving antennas. A study was 
performed to see if the satellite constellation could be use- 
fully augmented with ground-based GPS transmitters for the 
ascent and post-blackout mission phases. 

There is a considerable connection between the R/PA 
software and the hardware signal tracking loops. To shed 
further light on this interaction and on the performance po- 
tential of the tracking loops, a study of one mechanization of 
the tracking loops was performed. That mechanization is asso- 
ciated with the GPSPAC receiver (being developed for space 
applications by Magnavox) which was an early candidate for a 
place aboard the Orbiter. 

Finally, two special topics related to Orbiter use of 
GPS were studied. One is the relativistic effects on the Or- 
biter clock used in the receiver to make GPS measurements. 

Due to the very large changes in velocity and gravitational 
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potential experienced by the Shuttle, the relativistic effects 
are much larger than they would be for earth-based users. 

Since an atomic time standard is being considered for Shuttle 
R/PA application, It may be desirable to explicitly account 
for these effects in the R/PA software. The other special 
topic addressed is geometric dilution of precision (GDOP). 

GDOP is not so meaningful for Shuttle operations as it is for 
other applications, since the anticipated receiver configura- 
tion for the Orbiter is a two-channel sequential set, since 
four GPS satellites may not always be available during the 
initial usage of GPS on the Orbiter, and since the potential 
exists for use of an atomic time standard in the receiver. 
Nonetheless, it is useful as a measure of relative navigation 
performance potential. 

1.1.2 Baseline System Navigation Studies 

In addition to the GPS studies, several of the studies 
summarized in this report address the baseline navigation sys- 
tem. Four specific areas are addressed; three of them involve 
calibration of the 1MU, while the fourth centers on navigation 
performance . 

Preflight calibration and alignment of the IMU's has 
been studied in prior contracts, but this report contains the 
first (final) report on covariance analysis of the current 
cal/align procedure. The previous efforts utilized monte carlo 
techniques to generate an overall performance projection for 
the Approach and Landing Tests (ALT) and to identify potential 
risk areas. The current effort uses covariance techniques to 
develop a detailed error budget for the Orbital Flight Tests 
(OFT's, now called Space Transportation System, or STS, mis- 
sions) and to identify the major error mechanisms that limit 
calibration and alignment performance. 
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The ALT’S provided NASA’s first opportunity to eval- 
uate both the navigation performance achievable with the IMU 
and the true calibration and alignment uncertainties in a 
flight environment. TASC’s study of postflight IMU error 
recovery provided a preliminary assessment of the potential 
accuracy of recovering IMU errors from ALT’S. 

After lengthy periods on-orbit, the IMU gyros and ac- 
celerometers will need recalibration of their only observable 
characteristics in the weightless environment: their biases. 

The study of just how to best perform such a calibration recom- 
mends procedures, timings, and analytical techniques based on 
laboratory models of the instruments. 

The final baseline navigation topic is the assessment 
of navigation performance during entry, prior to the landing 
phase. This work, too, was begun under the previous contract, 
but is substantially revised here, based on more recent assump- 
tions of initial errors and sensor performance. A candidate, 
baseline navigation filter is examined in detail. 


1.2 ORGANIZATION OF THE REPORT 

This report is divided into two volumes. Volume 1 
contains all the GPS-associated studies, while Volume II pre- 
sents those studies addressing the baseline system. 

In Volume I, Chapter 2 summarizes results pertaining 
to Shuttle navigation with Phase I and II GPS, the test con- 
stellations of navigation satellites. Chapter 3 is a study of 
Shuttle navigation on orbit with a single GPS satellite. 
Chapter k discusses post-blackout navigation with the fully 
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operational. Phase III GPS constellation. Chapter 5 considers 
possible augmentation of the GPS satellite constellation with 
GPS ground transmitters for the ascent and post-blackout mis- 
sion phases. Chapter 6 contains an analysis of a particular 
version of a GPS receiver, the GPSPAC receiver, discussing 
tracking loop performance in the Shuttle environment. Chapter 
7 summarizes the most significant conclusions of Volume I. 
Appendix A considers the relativistic effects on the Space 
Shuttle GPS clock that must be part of the GPS user equipment. 
Appendix B is a tutorial discussion of the subject of dilution 
of precision as it applies to GPS. 

In Volume II, Chapter 2 reports the final results of 
the revised IMU calibration and alignment covariance analysis. 
Chapter 3 considers postf light IMU error recovery for ALT. 
Chapter 4 summarizes the study of on orbit calibration of IMU 
accelerometer and gyroscope biases. Chapter 5 discusses the 
baseline entry and preland navigation performance. Finally, 
Chapter 6 summarizes the most significant conclusions of 
Volume II. 
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2. NAVIGATION WITH PHASE I AND II GPS 


This chapter considers the problem of navigating the 
Space Shuttle with the developmental satellites of the Phase I 
and II GPS constellation. During the time in which much of the 
analysis was being done (Refs. 5 through 8), the concept of using 
GPS aboard the Space Shuttle was evolving rapidly, and the ana- 
lytical baseline was changing in parallel. Those study results 
still relevant to the Space Shuttle program are presented in 
this chapter, accompanied by some partly outdated, isolated 
segments believed to still have current, if limited, value. 

Section 2.1 provides a brief introduction to the remain- 
der of the Chapter. Section 2.2 is a lengthy summary of GPS 
satellite visibility with the Phase I and II GPS constella- 
tions. Sections 2.3 and 2.4 summarize navigation performance 
results for the deorbit and on-orbit Shuttle mission phases. 


2 . 1 INTRODUCTION 

Both the plans for the GPS itself and for its use aboard 
the Space Shuttle have been evolving over the three year period 
encompassed by this report. It is still expected that there 
will be a period, possibly of several years duration, in which 
the GPS user equipment will be aboard the Orbiter, but the 
fully operational GPS satellite constellation of Phase III 
will not be in place. The ability of the Shuttle to navigate 
during this period with GPS is open to considerable question, 
as there will generally be fewer satellites available than are 
sufficient to support accurate navigation. 
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This chapter assesses the navigation performance 
achievable before the Phase III GPS constellation is avail- 
able. Since the navigation performance depends so much on the 
limited GPS signal availability, that issue -- satellite visi- 
bility --is treated at some length. Visibility results for 
an example of the deorbit and on-orbit mission phases are used 
as inputs to statistical simulations of navigation performance. 
In addition, the effects of antennas on the visibility situ- 
ation are examined (although the antennas considered here do 
not match those currently planned for the Orbiter), and a 
parametric study of orbital visibility is made covering the 
constellation options for GPS Phase II. 

The navigation performance studies posit navigation 
filter algorithms for the Shuttle GPS receiver and evaluate 
those algorithms via a filter-model-versus-truth-model co- 
variance analysis approach. Sensitivity studies are done on 
some key assumptions. Both the deorbit and on-Oibit flight 
phases are considered. 

A previous report (Ref. 4) considered GPS-aided navi- 
gation in the ascent phase, but without regard to GPS receiv- 
ing antenna blockage by the external fuel tank of the Orbiter. 
The presence of the tank makes the use of GPS during ascent 
more difficult. Chapter 5 of this report considers ..ugm :'ng 
the satellite constellation with ground GPS transr tters to 
support the Shuttle ascent phase. 


2.2 PHASE I AND II GPS SATELLITE VISIBILITY 

GPS transmitting antennas are pointed at the earth 
and have beamwidths adequate to cover the entire earth and 
near-earth region, including the outer limits of Space Shuttle 
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travel. The L-band GPS signals have wavelengths of about 20 cm, 
long enought to penetrate almost all weather conditions, but 
short enough so that signal transmission is largely a line-of- 
sight matter. To users with full, spherical, receiving antenna 
coverage, signal availability is a simple matter of satellite 
visibility. If a satellite is visible from the user's antenna 
position, then its signal is available to the user's receiver. 
The only reason the satellite might not be visible is that the 
earth or some other body is blocking the signal path. 

The ability of GPS to aid in Space Shuttle navigation 
depends strongly on the number of signals available to the re- 
ceiver. Having four signals available implies good navigation, 
except for situations of singular or nearly-singular geometry. 
Having more than four :;ignals available increases the likelihood 
of finding a suitable subset of four. When fewer than four sat- 
ellites are available, however, the situation is more complex. 

When fewer than four GPS signals are available to the 
Shuttle receiver, the quality of navigation will depend on 

• The past history of signal availability 

• The quality of the clock (oscillator) in 
the receiver 

• The magnitude of the unaccountable portion 
of disturbing accelerations on the Orbiter 
(due primarily to venting on orbit or to 
1MU errors during dynamic flight phases) 

• The degree to which the satellite-to-Orbiter 
geometry has varied in the recent past 

• The type of navigation filter in the GPS 
receiver. 


*See Appendix B for a discussion of such situations. 
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Since navigation performance depends so strongly on 
signal availability, visibility studies are an important prel> 
ude and accompaniment to navigation performance studies. Such 
signal a\ailability studies are all termed here "visibility" 
studies, although they sometimes go beyond the simple concept 
of line-of-sight visibility into signal availability, and even 
into signal usability, if necessary. The main factors in the 
visibility analyses are: 

For line-of-sight visibility 

The positions of the GPS satellites 

The position of the Orbiter 

The location of the earth and other sig- 
nal blocks such as the vehicle structure, 
the external fuel tanks, and the reentry 
plasma sheath. 

When antenna effects are considered 

The locations of the antennas on the 
Otbiter 

The attitude of the Orbiter 

The gain patterns of the antennas. 

When signal dynamics are considered 

The capability of the receiver tracking 
loops 

The ability to switch from one antenna 
to another. 

2.2.1 GPS Satellite Positions 


The GPS is undergoing a phased development. In the 

JL 

final, operational, Phase 111, there are to be 24 satellites" 


*A proposed program restructuring would reduce the number of 
satellites to 18. 
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"evenly" distributed in a constellation that will provide good 
signal coverage everywhere on the earth's surface at all times. 
Once this constellation is achieved, any limitations on satel- 
lite visibility from the Space Shuttle will come from less 
than full coverage by the receiving antennas, or from antenna 
blockage by the external fuel tank during ascent or by the 
plasma sheath during reentry blackout. 

Schedules for GPS development and for Shuttle GPS 
user equipment development are still in flux, but it is expec- 
ted that there may be a considerable period of operation of 
GPS user equipment aboard the Orbiter before the operational 
constellation is achieved. During GPS developmental Phases I 
and II (Phase II began in 1979) there will be from four to six 
satellites "bunched together" so that four or more are visible 
at certain times and places, while fewer or none may be visible 
at other times and places. Generally speaking, good visibil- 
ity will be obtained from four broad areas of the earth (one 
of which is the Yuma, Arizona GPS testing grounds) for periods 
of a few hours each day. 

GPS satellite orbits are all nominally circular with 
radius somewhat over four earth radii, selected so that the 
satellite period is almost exactly one half of a sidereal day. 
Each sidereal day, the satellites make two complete revolutions 
while the earth completes one rotation on its axis, so the 
ground tracks of the satellites repeat. In fact, the orbital 
periods are adjusted so that the satellite ground tracks do 
repeat, taking into consideration a small retrograde movement 
of the satellite orbital planes brought about by the oblateness 
of the earth. Each ordinary, solar, day each satellite will 
appear over any given point on its repeating ground track just 
4 min, 3.4 sec earlier than it did the day before. Viewed from 


2-5 



THE ANALYTIC SCIENCES CORPORATION 


any point on the earth, the time when the Phase II satellites 
bunch together for multiple signal availability steadily gets 
earlier each day, taking a little over a year to come full 
cycle through a complete solar day. 

For most of the visibility studies reported here, the 
trajectory of the Space Shuttle is related to features fixed 
to the earth, but the time and date of the mission are com- 
pletely unknown. It is impossible to predict, then, what the 
visibility situation will be, even approximately. A complete 
range of visibility situations must be examined, parameterized 
on the time connection between the Shuttle and GPS satellite 
orbits. This time connection is termed "offset time" and runs 
through one full day before the visibility situation approx- 
imately repeats. That is, if the "starting time" of a mission 
under consideration were to vary over a span of one day, then 
a range of visibility situations would be encountered, with 
offset time being related directly to the mission starting 
time. Another useful interpretation of offset time is this: 
if the mission starting time is a fixed time of day, but the 
date of the mission were to vary over the span of one year, 
then the same range of visibility situations would be encoun- 
tered, with offset time now corresponding to the mission date 
in an obvious, if indirect, manner. The initialization of the 
parameter "offset time" for the studies reported here it gen- 
erally based on a nominal mission starting time and a satel- 
lite constellation epoch time that are no longer meaningful. 

No importance is ever attached to results for such arbitrary 
times; rather, the range of results is considered for vari- 
ations . 


All of the satellite orbits have the same nominal 
inclination angle to the earth's equator. Initially, this 
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inclination angle was planned to be 63 deg, and the satellites 
of the Phase I & II constellation have (or will have) this 
nominal inclination. Many of the results reported in Chap* 
ter 4 for the Phase III constellation also assume a 63 deg 
inclination. Recent planning for Phase III, however, calls 
for a 55 deg inclination angle to improve visibility from mid- 
dle latitudes while sacrificing some visibility from polar 
regions. Only minor differences in visibility are observable 
between these two inclination angles. 

The Phase III constellation has its 24 satellites 
arrayed as eight satellites in each of three orbit planes. 

The eight satellites are equally spaced, 45 deg from one to 
the next, in an orbit plane. The three planes are equally 
spaced, 120 deg from one to the next, in right ascension from 
the vernal equinox. The relative positions of the eight satel- 
lites in one plane, compared to those in another, are offset 
by 15 deg, so that each satellite has a unique ground track on 
the earth, although there are only three orbital planes in 
inertial space. Table 2.2-1 shows the initial right ascen- 
sions and initial anomalies of the 24 Phase III satellites. 

Both the initial epoch (the time at which the values the table 
are valid) and the reference point for the right ascensions 
are regarded as arbitrary for purposes here. 

The plans during Phase I were to array six satellites 
as three in each of two orbit planes. The three satellites in 
each plane were to be separated by 40 deg, the right ascensions 
of the planes were to be separated by 120 deg, and the relative 
positions of satellite anomalies in one plane versus those in 
the next t,ere to be 40 deg different. Table 2.2-2 summarizes 
this situation, and is the constellation assumed in this report 
whenever Phase I is referenced without qualific tion. 
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TABLE 2.2-1 

RIGHT ASCENSIONS AND INITIAL ANOMALIES 
FOR THE PHASE III CONSTELLATION 


ORBIT PLANE 1 * 
RIGHT ASCENSION 
0 deg 

ORBIT PLANE 2 
RIGHT ASCENSION 
120 deg 

ORBIT PLANE 3 
RIGHT ASCENSION 
240 deg 

SATELLITE 

NUMBER 

INITIAL 

ANOMALY 

SATELLITE 

NUMBER 

INITIAL 

ANOMALY 

SATELLITE 

NUMBER 

INITIAL 

ANOMALY 

1 

0 

9 

15 

17 

30 

2 

45 

10 

60 

18 

75 

3 

90 

11 

105 

19 

120 

4 

135 

12 

150 

20 

165 

i 

5 

180 

13 

195 

21 

210 

6 

225 

14 

240 

22 

255 

7 

270 

15 

285 

23 

300 

8 

315 

16 

330 

24 

345 


♦Arbitrarily chosen reference point, not necessarily referenced 
to the vernal equinox. 


TABLE 2.2-2 

RIGHT ASCENSIONS AND INITIAL ANOMALIES 
OF THE PHASE I CONSTELLATION 
(BASED ON 1978 PROJECTIONS) 


ORBIT PLANE 1* 
RIGHT ASCENSION 
-130 deg 

ORBIT PLANE 2 
RIGHT ASCENSION 
+110 deg 

SATELLITE 

INITIAL 

SATELLITE 

INITIAL 

NUMBER 

ANOMALY 

NUMBER 

ANOMALY 

1 

0 

4 

40 

2 


5 

80 

3 

80 

6 

120 


^Referenced to the location of the Greenwich 
Meridian at 2 hours GMT on 26 August 1977. 
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Figure 2.2-1 illustrates the ground tracks of the 
Phase 1 satellites. The circles on the ground tracks corre- 
spond to the satellite positions at the beginning of the 
Orbiter deorbit trajectory segment considered in the next 
section. 


R -2 7404 



Figure 2.2-1 Ground Tracks of the Phase I 

GPS Satellites 

The tolerances on the orbital parameters of the GPS 
satellites at the Phase 1 epoch are shown in Table 2.2-3. The 
variation in the orbital period is especially important. During 
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TABLE 2.2-3 

TOLERANCES ON THE GPS PHASE I 
CONSTELLATION AT EPOCH (Ref. 6) 


Period 

±2 sec 

Inclination 

±2 deg 

Eccentricity 

< 0.02 

Longitude of the 

±2 deg 

Ascending Node 

Angle Between Any 

±5 deg 

Two Satellites 


the approximately 580 days from the GPS Phase I epoch to the 
nominal mission dates considered in the following sections, 
the orbital period variations could cause satellite mean anom- 
aly variations of over 19 deg. 

GPS program plans were revised before all of the satel 
lites were launched. Table 2.2-4 summarizes the revised Phase 
I constellation, having 42 and 43 deg separations between satel 
lites in one orbit plane and a 47 deg difference in relative 
anomalies from satellites in one plane to those in the other. 
The initial anomalies of the table apply at the epoch time of 
midnight GMT on 21 March 1977. At this time the right ascen- 
sion of the Greenwich meridian was 178.4 deg. 

The four satellites actually launched to date were to 
nominal positions corresponding to those characteristic of 
Table 2.2-4, and they are identified by the "NS” numbers in 
parentheses in the table. Where only four satellites are con- 
sidered as a constellation in this chapter, it is these four. 
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TABLE 2.2-4 

RIGHT ASCENSIONS AND INITIAL ANOMALIES 
FOR THE REVISED PHASE I CONSTELLATION 


ORBIT PLANE 1* 
RIGHT ASCENSION * 
219.49 deg 

OP3IT PLANE 2 
RIGHT ASCENSION 
99.49 deg 

SATELLITE 

INITIAL 

SATELLITE 

INITIAL 

NUMBER 

ANOMALY 

NUMBER 

ANOMALY_ 

l(NS-l) 

-29 

4 

18 

2 

14 

5(NS-2) 

61 

3(NS-4) 

56 

6(NS-3) 

103 


^Referenced to the mean vernal equinox of 1950.0. 

2.2.2 Phase I Deorbit Visibility 

Satellite visibility during the deorbit segment of a 
Space Shuttle trajectory is considered in this section. Only 
one trajectory is used for analysis, and that is to a temporary 
landing site, but the results should be indicative of perform- 
ance on most deorbit trajectories from similar orbital inclina- 
tions. Using only one trajectory has the advantage that direct 
comparisons of results obtained under a variety of assumptions 
is not further complicated by trajectory differences. Naviga- 
tion performance associated with these visibility results is 
reported in Section 2.3. 

Deorbit Trajectory - A single, simulated Orbiter tra- 
jectory generated at JSC is used throughout this report for 
all studies of deorbit navigation. This trajectory was gen- 
erated as a candidate for an OFT-1 (STS-1) mission segment, 
and covers the span from just prior to the deorbit burn at 655 
kft of altitude, through touchdown at Edwards AFB runway 17. 
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The trajectory ground track is plotted on a world map in Fig. 
2.2*2, beginning in Australia and continuing to California. 

Also plotted are the nominal ground tracks of the Phase I GPS 
satellites over the same period. Figure 2.2*3 is a plot of 
Orbiter altitude versus time. Both figures show a nominal 
span labeled "blackout" referring to an assumed loss of L*band 
(GPS and TACAN) signals during reentry. Additional information 
relevant to this trajectory is found in Chapters 4 and 5 of 
this volume, in Chapter 5 of Volume 2, and in Refs. 1 through 
4, and 17 through 19. 


R-2741? 



Figure 2.2-2 Ground Tracks of the Shuttle Deorbit 

Trajectory and Phase 1 GPS Satellites 
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Figure 2.2-3 Altitude versus Elapsed Time 

for the Deorbit Trajectory 

The specific epoch assigned to this trajectory is 16 
March 1979 at 19:16:40.4 GMT, and this epoch is used for the 
Phase I and II visibility studies of this chapter. While the 
trajectory continues down to a landing, it is sometimes trun- 
cated at 20 kft of altitude in this report when the landing 
phase is not of interest. Besides including Orbiter position, 
velocity, and sensed acceleration as functions of time, the 
trajectory tape contains a time history of Orbiter attitude 
(and numerous other quantities) useful in assessments of visi- 
bility with specific antenna locations. 

Deorbit Visibility Results - In the visibility re- 
sults presented in this section, a satellite is assumed visible 
from the Orbiter if it is more than 10 deg above the Orbiter* s 
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horizon. No antenna limitations are considered. The figure 
of 10 deg, called the visibility mask angle or the earth mask 
angle, is somewhat arbitrary; values from 5 deg to 15 deg are 
reasonable, depending on such factors as antenna gain and at- 
mospheric attenuation of the GPS signals. An approximation to 
the ellipsoidal earth was used which yields horizon angle com- 
putation accuracy of 0.1 deg or better. 

Figure 2.2-4 shows how many satellites were visible 
at each minute from the deorbit epoch down to touchdown on 
runway 17 at Edwards AFB, assuming the nominal epoch times for 
the trajectory segment and for the Phase I constellation. As 
illustrated later, the touchdown time selected for the trajec- 
tory tape corresponds to a period of extremely poor GPS satel- 
lite visibility along the approach path. 
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Figure 2.2-4 


Satellite Visibility During the 
Deorbit Phase at Nominal Epochs 


During a portion of the deorbit phase of the test, 
there may be a blackout of radio communications, meaning that 


2-14 


THE ANALYTIC SCIENCES CORPORATION 


measurements to the GPS satellites are not possible, whatever 
the geometry for visibility. The blackout portion of the test 
is indicated in Fig. 2.2-4, based on the assumption of blackout 
conditions between 300,000 ft and 130,000 ft of altitude during 
deorbit. 


If the mission timeline were advanced so that the de- 
orbit epoch and touchdown occurred at some later time than the 
one given, then the satellite positions in inertial space would 
be different, as the satellites advance in their orbits. Also, 
the Orbiter position in inertial space (given on the trajectory 
tape) would have to be rotated with the earth so that the Or- 
biter would still land at Edwards AFB. Precisely these calcu- 
lations were carried out for a sequence of possible offset times, 
a sequence spanning one day in steps of 30 min. The calculations 
leading to the visibility plot of Fig. 2.2-4 were repeated for 
each offset time to examine satellite visibility as a function 
of deorbit epoch. Figure 2.2-5 presents all these results. 

A vertical line in Figure 2.2-5 shows how many satel- 
lites would be visible to the Orbiter as a function of elapsed 
time since the deorbit epoch, for one particular deorbit epoch. 
The leftmost vertical line corresponds to the originally quoted 
epoch, and thus reproduces the data of Fig. 2.2-4. Each succeed- 
ing vertical line corresponds to a 30 min offset of the epoch. 

As the figure shows, if the deorbit epoch were chosen to maxi- 
mize GPS satellite visibility, it would be about 16 hours later 
than the time quoted. Figure 2.2-6 shows the visibility for 
the deorbit phase if the deorbit epoch were offset 16 hours, 
and the figure includes a bar plot showing which satellites 
are visible at which times. Figure 2.2-7 is a second such 
plot showing reduced satellite visibility for an offset time 
of 21 hr 43 min. Both the good visibility segment and the re- 
duced visibility segment are used in the deorbit navigation 
studies of Section 2.3. 
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Figure 2.2-5 Satellite Visibility During the Deorbit 

Phase as a Function of Offset Time 
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Figure 2.2-6 Good Satellite Visibility During Deorbit 

for a 16 hr Offset Time 
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Figure 2.2-7 Reduced Satellite Visibility During Deorbit 

for a 21 hr 43 m-n Offset Time 
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Instead of considering the offset time as an advance 
in the Orbiter epoch, it is equally valid to consider it as a 
regression of the GPS Phase 1 epoch, or as a combination of 
the two. Because the satellite constellation at epoch was 
described in earth relative terms, the visibility results are 
quite the same if the Orbiter is late as they are if the sat- 
ellites are early; it is only the elapsed time from the sat- 
ellite epoch to the Orbiter epoch that determines satellite 
visibility. 

2.2.3 Phase I Orbital Visibility 

Visibility results for the Shuttle on orbit are re- 
ported in this section in a manner similar to that used for 
the deorbit results of the previous section. One orbital tra- 
jectory is baselined for visibility results that are later 
input to navigation performance analyses. Visibility results 
for the Shuttle over a range of other orbits, considering an- 
tenna effects explicitly, are presented in Sections 2.2.4 and 
2.2.5. 


Orbital Trajectory - The orbital trajectory used in 
the study of orbital navigation is extrapolated from a set of 
initial conditions (position and velocity) by assuming a simple 
earth gravity model (the same one used for GPS satellite orbit 
propagation) with no sensed acceleration. Table 2.2-5 sum- 
marizes the orbital initial conditions, while Fig. 2.2-8 plots 
a 90 min segment of the orbit ground track on a world map. 

Orbital Visibility Results - Figure 2.2-9 shows how 
many satellites were visible at 5 min intervals over the 12 
hour (8 orbit) period commencing at orbital epoch. Since the 
GPS satellites have periods of almost exactly 12 hours, which 
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TABLE 2.2-5 

ORBITAL TRAJECTORY PARAMETERS 


him 


ORBIT DATA 

30 March 1979 

INITIAL TIME 

13:15:38.8 GMT 

INITIAL POSITION* 

(-2081432, 20009848, - 8493118) ft 

INITIAL .'ELOCITY* 

(-21116, -7323, -12066) fps 

ALTITUDE AT APOGEE 

155.3 n mi 

ALTITUDE AT PERIGEE 

149.9 n mi 

ORBITAL INCLINATION 

37.9 deg 


*With respect to Aries-Mean-of-1950 coordinates. 





Figure 2.2-8 Shuttle Orbital Trajectory 
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Figure 2.2-9 Satellite Visibility During 

the Orbital Phase 


also corresponds to an integral number of Orbiter revolutions, 
the visibility patterns should almost repeat every 12 hours. 
The figure indicates that, in most instances, four or more 
satellites are simultaneously visible at some time during each 
orbit. Three or more satellites are visible about 50% of the 
time . 


As was done for the deorbit phase, the orbital re- 
sults were recomputed for a 24 hour sequence of offset times 
at 30 min intervals to assess visibility if the orbital epoch 
or the GPS Phase I epoch should change. Figure 2.2-10 (2 pages) 
presents these orbital visibility results in the same form 
that Fig. 2.2-5 used for the deorbit phase. 

Visibility on orbit, while it depends in detail on 
the offset time, is usually good for a significant portion of 
each orbit whatever offset time is chosen. In order to make 
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Figure 2.2-10 Satellite Visibility During the Orbital 

Phase as a Function of Offset Time 
(See key on Fig. 2.2-5) 
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use of the results of a detailed statistical simulation of the 
GPS satellite ephemeris and clock errors done for the deorbit 
navigation study it was necessary to assume an offset time of 
16^58 m 58f4 (not plotted in Fig. 2.2-10) for use in the orbital 
navigation study. For the most part, this assumption has little 
effect on the significant features of orbital navigation. The 
visibility situation for the above offset time is illustrated 
in Fig. 2.2-11. The lower half of the figure shows which sat- 
ellites are available at which times and also indicates, by 
the tic marks, a measurement schedule used in the orbital nav- 
igation study of Section 2.4. 


A*UM 




6 LAPSED TIME DURING ORBIT (mm) 

Figure 2.2-11 Orbital Satellite Visibility and 

Measurement Schedule 


2.2.4 Antenna Effects on Phase I Orbital Visibility 

During the evolution of the concept of Space Shuttle 
navigation using GPS, there was an experimental stage consid- 
ered as a part of a phased development plan. That phased plan 
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was later dropped, as was the experiment. Some of the visi- 
bility results developed for the experiment are reported here 
for their general insight into orbital visibility considering 
antenna effects. These results are isolated to this section, 
and do not form the basis for further analysis in this volume. 

Experimental Antenna Types and Locations - Two types 
of GPS antennas were considered for experimental use aboard 
the Orbiter. One type is an antenna being constructed by APL 
for use with the GPSPAC receiver under development by NASA 
Goddard and the DMA for use in satellites. The gain pattern 
of this antenna is shown in Fig. 2.2-12 for right hand cir- 
cularly polarized signals at the GPS frequency. (Since 
signals were originally planned to be weaker than signals, 
the I <2 performance was deemed more critical than the per- 
formance.) The coverage of this antenna may be approximated 
by the specification, which calls for a conical zone of cover- 
age with half-cone angle of 100 deg. This antenna would have 
to be located in the payload bay of the Orbiter or be deployed 
from the bay on a boom. 

The second type of antenna considered is a stripline 
antenna that would be mounted on the inside of a window into 
the crew area. The coverage zone of such an antenna is uncer- 
tain without testing. It is modeled here as having a conical 
coverage zone with a 45 to 75 deg half-cone angle. 

Figure 2.2-13 shows five locations considered for 
placement of the experimental antennas. Locations I and II 
are on pallets in the payload bay, appropriate for the APL 
antenna. Locations III and IV are alternate positions of a 
boom extending from the payload bay, also for the APL antenna. 
Location V is a window into the crew area, appropriate for the 
stripline antenna. 
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Figure 2.2-12 APL Antenna Pattern (RHC Gain for l^) 


ft -29 999 



Figure 2.2-13 Alternate Experimental GPS Antenna Locations 
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The locations I and II in the payload bay suffer from 
masking caused by the Orbiter structure. Figures 2.2-14 and 
2.2-15 show the antenna masking by the Orbiter for the APL 
antenna in locations I and II. Any direction from the antenna 
may be described by two angles termed "look angles." Refer- 
ring to the standard Orbiter body coordinate system with its x 
axis out the nose, the y axis out the right wing, and the z % 
axis out the bottom, the look angles are defined as: pitch 
look is the elevation angle of the antenna look direction above 
the xy body plane; yaw look is the angle the projection of the 
look direction from the antenna onto the xy plane is rotated 
about the positive z axis from the x axis toward the y axis. 

The figures consider all possible directions from the antenna 
and show which are masked by the vehicle itself (the hatched 
region at the bottom of the figure) and which, although not 
masked, are outside the 100 deg half cone angle limit for the 
APL antenna (the shaded region below -10 deg of pitch look 
angle). The Orbiter tail is evident in the figures on the left 
hand side; it is not symmetrically located about the 180 deg 
yaw look angle because the antenna locations are slightly to 
the right of the Orbiter centerline (i.e., have +y components 
of location in body coordinates). The figures were prepared 
from equivalent data in a different form (Ref. 20). 

Experimental Orbiter Trajectories - The missions assumed 
for the Shuttle GPS experiment were denoted as OFT-5 and OFT-6. 
OFT-5 was to have been a Skylab revisit launched on 3 February 
1980. Even as this work was progressing, the mission plans were 
changing. Originally the OFT-5 mission was to contain lengthy 
segments with the Orbiter in a bottom-down, earth oriented at- 
titude, and some of the results reported here assume such an 
attitude. The plans changed so that during most of OFT-5 the 
Orbiter would be stabilized in inertial space with its wing 
toward the sun or with its tail toward the sun. At this point, 
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Figure 2.2-14 APL Antenna Masking Due to 

the Orbiter -- Location I 



*10 dtg 
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Figure 2.2-15 APL Antenna Masking Due to 

the Orbiter -- Location II 


interest was expressed in OFT-6, since it was to have lengthy 
segments of the desirable, bottom-down, attitude stabilization 
that improves satellite visibility from topside antennas. Table 
2.2-6 shows the basic orbital parameters for the two missions. 
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TABLE 2.2-6 

ORBITAL PARAMETERS FOR OFT- 5 AND OFT- 6 


PARAMETER 

OFT-5 

^OFT-^6 ^ 

Altitude 

175 nm 

160 nm 

Eccentricity 

0 

0 

Inclination 

50 deg 

57 deg 

Period 

90.9 deg 

90.45 min 


Experimental Orbital Visibility Results - Figure 2.2-16 
shows the nominal satellite visibility from the Orbiter for OFT-5 
without regard to antenna considerations. The only visibility 
criteria is a 10 deg earth mask angle. To see how such full 
coverage visibility changes with offset time. Fig. 2.2-17 is a 
visibility density plot showing the number of satellites vis- 
ible as functions of both time on orbit and offset time. 




10° EARTH MASK ANGLE 



Figure 2.2-16 GPS Phase 1 Satellite Visibility 

During OFT-5 Orbit 
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2.2-17 GPS Phase I Visibility from OFT-5 
on Orbit as a Function of Elapsed 
Time and Offset Time 
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To confirm that important orbital visibility con- 
clusions are largely independent of offset time, some average 
visibility results were prepared. Figure 2.2-18 shows the 
average number of satellites visible (averaged over time on 
orbit) as a function of offset time. Figure 2.2-19 shows what 
percentage of time on orbit three or four satellites are vis- 
ible as a function of offset time. The relatively small vari- 
ations in these plots confirm that while visibility is a 
detailed function of offset time, the averages are much less 
sensitive to offset time. 



Figure 2.2-18 Average Visibility for OFT-5 as Function 

of Offset Time 
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Figure 2.2-19 Percentage Visibility for OFT-5 as 

Function of Offset Time 

Figures 2.2-20 and 2.2-21 show example visibility 
penalities associated with the APL antenna placed in locations 
I and II. These figures were prepared when the baseline for 
OFT-5 was a bottom-down attitude. The "top side" of the plot 
is the same in both figures and is identical to Fig. 2.2-16, 
for it represents full coverage visibility considering only 
earth masking with the 10° mask angle above the Orbiter's 
horizon. The shaded portions of the figures represent sat- 
ellites "lost" to the two antenna placements due to Orbiter 
masking (the antenna beam pattern was never a factor). Both 
locations offer nearly full coverage, with placement I being 
slightly superior as it is farther from the tail and suffers 
less from tail masking. 
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APL ANTENNA PLACEMENT I 



012345 678 

ELAPSED TIME ON ORBI f (hrs) 

( SLACK REGION IS VISIBILITY LOST FROM FULL COVERAGE ) 


2.2-20 Example of Orbiter Masking Visibility 
Penalty for Bottom-down Attitude -- 
Location I 


APL ANTENNA PLACEMENT II 



ELAPSED TIME ON ORBIT (hrs) 


( BLACK REGION IS VISIBILITY LOST FROM FULL COVERAGE I 

Figure 2.2-21 Example of Orbiter Masking Visibility 

Penalty for Bottom-down Attitude -- 
Location II 
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Table 2.2-7 compares full coverage visibility with 
that from locations I and 11 and from a location V window 
antenna of three possible different half-cone angles. This 
table still assumes that the OFT-5 mission would be flown 
bottom-down. 


TABLE 2.2-7 

ANTENNA EFFECTIVENESS COMPARISON FOR PRIOR 
OFT-5 BASELINE (BOTTOM-DOWN) 


ANTENNA 

AVERAGE NUMBER 
SATS VISIBLE 

X TIME 
NONE VISIBLE 

7. TIME 
3 OR MORE 

X TIME 
4 OR MORE 

IDEAL 
(APL IV) 

2. 52 

12 

51 

26 

APL I 

2*39 

13 

49 

22 

AFL II 

2.27 

15 

45 

20 

75° HCA 

1.16 

38 

12 

2 

60° HCA 

0.69 

54 

3 

0 

45° HCA 

0.37 

70 

0 

0 


HCA • HALF* CONE* ANGLE WINDOW ANTENNA (V) 


Before proceeding with visibility results, the re- 
lated matter of kinematic ranges should be discussed. Any GPS 
receiver will have limitations on the signal dynamics that it 
can handle. These limits may usually be expressed in terms of 
the range from the receiver to the satellite and the first 
three derivatives of this range: range-rate, range-accelera- 

tion, and range- jerk. Table 2.2-8 gives the extreme values of 
range (and its first three derivatives) from the Orbiter to 
the visible satellites for a 12 hr segment of the OFT-5 orbit. 
The effect of vehicle maneuvers are not included in the values 
given in the table. Due to unbalanced attitude thrusting or 
to high attitude rates, the range jerk might, for very short 
periods of time, be significantly larger than the tabulated 
values, but the other variables should be largely unaffected. 
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TABLE 2.2-8 
KINEMATIC RANGES 


KINEMATIC VARIABLE 

OFT-5 ( 
EXT! 

DN ORBIT 
REMES 

MINIMUM 

MAXIMUM 

Range (10 6 ft) 

65.1 

87.6 

Range Rate (10 ft/sec) 

-25.9 

26.0 

2 

Range Acceleration (ft/sec ) 

-5.3 

41.2 

Range Jerk (10*^ ft/sec^) 

-49.9 

49.8 


Figures 2.2-22 through 2.2-24 show what nominal visi- 
bility may be expected on the OFT-5 orbit with the stripline 
window antenna in location V. Various values for the antenna 
half -cone angle are considered. These figures should be com- 
pared with Fig. 2.2-16 showing visibility for the same orbital 
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Figure 2.2-22 GPS Phase I Satellite Visibility from OFT-5 

Orbit with a 45° Half-Cone Window Antenna 
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Figure 2.2-23 GPS Phase 1 Satellite Visibility from OFT-5 

Orbit with a 60° Half-Cone Window Antenna 
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Figure 2.2-24 GPS Phase I Satellite Visibility from OFT-5 
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segment without regard to antenna limitations. Again, a bottom- 
down attitude of the Orbiter is assumed. Note that no more than 
two satellites can ever be seen at once with a 45 deg half-cone- 
angle antenna, and that 70% of the time none are visible. 

Figure 2.2-25 shows a six hour segment of visibility 
on orbit for OFT-6 as a function of both elapsed time and off- 
set time. Omitting some of the detailed results, Table 2.2-9 
summarizes visibility for OFT-5 and OFT-6 for antenna loca- 
tions I and II, now assuming that the Orbiter attitude on 
OFT-5 will be space stable while on OFT-6 it will be earth 
stabilized bottom-down. The conclusion to be drawn is that 
the space stable attitudes cut average visibility by 50% and 
reduce the amount of time 4 or more GPS satellites are visible 
by a factor of 3 to 4, so that bottom-down, earth stabilized 
attitude is desirable for a GPS demonstration. 


2.2.5 Phase II Orbital Visibility 

This section summarizes the results of a parametric 
look at GPS satellite visibility from the Space Shuttle on 
orbit. The satellite constellation used is an approximation 
to the GPS Phase II constellation of six satellites, with re- 
sults also for a four satellite subset to illustrate a minimum 
visibility situation. Full coverage antennas are assumed for 
the Orbiter, and a 5 deg earth mask angle is used. This study 
is separate from the other results of this chapter in that 
these particular visibility results were not used to generate 
navigation performance results. 

A conceptually straightforward way to organize a study 
of orbital visibility would be to fix a satellite constellation 
and then vary the orbit of the Space Shuttle. In general, 
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Figure 2.2-25 GPS Phase I Visibility from OFT-6 
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TABLE 2.2-9 

VISIBILITY COMPARISON FOR CURRENT 
OFT- 5 AND OFT- 6 MISSIONS 
(FOR NOMINAL ORBITS OF 12 HR DURATION) 


T-1478 


MISSION 

ANTENNA 

AVERAGE NO. 

X TIME 

X TIME 

X TIME 

(ATTITUDE) 

LOCATION 

SATS VISIBLE 

NOME VISIBLE 

3 OR MORE 

4 OR MORE 

OFT- 5 

Ideal 

2.5 

12 

51 

26 

Wing Co 

I 

1.2 

41 

21 

8 

Sun 

II 

1.2 

42 

20 

6 

Tail to 

I 

1.3 

39 

22 

8 

Sun 

II 

1.2 

40 

22 

7 

OFT- 6 

Ideal 

2.7 

11 

53 

30 

Boccom co 

I 

2.5 

13 

50 

28 

Earch 

II 

2.3 

16 

46 

23 


there are six parameters necessary to describe the Shuttle 
orbit, far too many to consider in this study. That number may 
be reduced by two by considering only circular Shuttle orbits, 
thus eliminating orbital eccentricity and argument of perigee 
from consideration. An additional parameter may be removed by 
fixing the orbit altitude at the typical value of 150 nm, as 
altitude does not strongly affect visibility. Orbital inclina- 
tion is a significant factor in the visibility patterns, so it 
must be retained as a variable parameter. Similarly, the right 
ascension of the Space Shuttle orbit with respect to the GPS 
satellite orbits is quite significant and must be retained as 
a variable. 

The initial anomaly of the Shuttle in its orbit is a 
significant factor in the visibility. When visibility is also 
parameterized as a function of time in orbit, however, and that 
time is carried over several complete revolutions of the Orbiter, 
then the resulting visibility patterns are nearly cyclical. 

The general pattern features under these conditions are not 
strongly dependent on the initial Shuttle anomaly, so it may 
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be removed as a free parameter. Thus, only two free parameters, 
inclination and right ascension, remain besides the independent 
variable of time in orbit. 

The Phase II visibility results were generated in a 
manner similar to that discussed above, but are organized dif- 
ferently. The concept of restriction to circular Shuttle 
orbits at 150 nm altitude was assumed. Shuttle orbital in- 
clination was varied over the discrete values of 0, 30, 45, 

60, and 90 deg. The initial anomaly of the Shuttle orbit was 
fixed at 0 deg. , so the Orbiter is always crossing the equator 
from south to north at time zero. 

The absolute right ascension of the Orbiter was not 
varied, always being assumed to be 0 deg with respect to the 
Greenwich meridian at time zero. In order to vary the Shuttle 
orbit right ascension with respect to the GPS satellite orbit 
right ascensions, the concept of "offset time," previously 
built into the software, was used. By specifying the initial 
longitudes of the ascending nodes of the GPS satellite orbits 
with respect to the Greenwich meridian, the right ascensions 
with respect to the vernal equinox vary with the satellite 
epoch time. Offset time is the span from the satellite epoch 
time to the time zero of the "time-on-orbit" parameter, so a 
positive offset time implies that the satellite epoch precedes 
the time zero. 

Table 2.2-10 gives the satellite constellation ini- 
tial longitudes of the ascending nodes and initial anomalies 
at the satellite epoch time. This constellation of six satel- 
lites is unique to this section; it approximates the Phase II 
constellation. The satellite inclinations are assumed to be 
63 deg. A reduced constellation of four satellites, omitting 
satellites 2 and 4, is used to simulate a possible reduction 
of the Phase II constellation. 
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TABLE 2.2-10 

SATELLITE CONSTELLATION FOR SAMPLE ORBITAL VISIBILITY RESULTS 


SATELLITE 

NUMBER 

INITIAL LONGITUDE 
OF THE ASCENDING NODE 
(wrt GREENWICH, deg) 

INITIAL 

ANOMALY 

(deg) 

1 

0 

0 

2* 

0 

45 

3 

0 

90 

JL 

4“ 

240 

30 

5 

240 

75 

6 

240 

120 


*These satellites are omitted for the four 
satellite version of this constellation. 

When the offset time parameter of the results is zero, 
the Shuttle begins its orbital travel directly beneath the 
Number 1 GPS satellite. As the offset time varies through 24 
hours (less 4 min, 3.4 sec), the right ascension of the Shuttle 
orbit varies through 360 deg with respect to that of the number 
1 GPS satellite. The offset time parameter also affects the 
time zero anomalies of the GPS satellites with respect to that 
of the Shuttle, but as previously noted, this variation has 
little effect on the general properties of the results. 

Figures 2.2-26 through 2.2-35 are visibility density 
plots for the orbital visibility studies. The title of each 
figure specifies the number of GPS satellites assumed to be in 
the constellation and the orbital inclination of the Space 
Shuttle. The vertical axes are time on orbit, from 0 to 3 
hours, spanning two complete Shuttle orbits. The horizontal 
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Figure 2.2-26 


Orbital Visibility, 0 Deg Inclination, 
4 Satellites 
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Figure 2.2-27 Orbital Visibility, 0 Deg Inclination, 

6 Satellites 
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Figure 2.2-28 


Orbital Visibility, 30 Deg Inclination, 
4 Satellites 


2-43 




4445444443322221111 1111 1 1 221 122222233333SS555555 
54455444444333311 II 1222 122221 2222222333334444455 
554555544433333322222222322222222111113333444445 
55555555553333332333322232222221 1 1 1 1 1223344445 

555665555544333433333322332222 111 122334445 

5556666 55444 2233233 3333233321 11 1111 2234445 

5556666654322233233333323332221 1 11 11 2233445 

5556666652233233233333343343221 111 111 1223345 

5556666332233233233333344443321 1111 1111223333 

44556333332332332333336554433331 1111 1111 13333 

3444333333333333333356655543333 22 1111 111111 333 

33223333334 3333333566665555433222211 11 11111 123 

22223333334333444566666555544422222111 11111 

122223333344445555666665554444322221 111 111 

1 122223333444555555666555544443322221 111 11 

1 1 122223344445555555554445444433322221 1 111111111 
1 1 1 122222333455555555544444444333322221 111111112 
2111 1222233334555444443444344 433344333332 1111111 
111111 122333334433334433433344444555433332221 111 
1 1111 1 11333333333333433333344456555553332221 11 1 
1 111 1223322233333323333332566666555443322211 1 

1 1 1222333333333323322224466666555544332211 1 

1 122343333333332322 134455666665554433221 1 1 

1 1333333333333321 1123445555666555544332111 

1 333333333333332 222333555555655554433331 1 

1 2333233333323332 1 1223334455556555554433322 

212333332332233222 1 1 123334445555655555533332 

3433333222322221 1 111 12 223 44 44 5 565 55 55 5 554 333 

444333322221221 11 1 122223444455556655555544 

444333322221 111 11 1122223344445555666655554 

444433322221 111 1111 22223 344445555566665555 

544443322221 111 1 122 1222223334444555555555455 

55444442221 111111111 2222322233333333335555555545 
55544433331 1 111 1 1 22232223222233221 23333445555545 
5555543333221112332233223322221 112333344555555 

655555433222 2333332333233321 11133334455556 

66555544321 1223333333332333 1 1122334455556 

666655443 1 233233333333343221 1111 2233445556 

6666644222233233333333332221 1 1111 2233445555 

66666322332332333333324332211 11 1 1 1 1223335555 

5665332333333333333344433321 11111 11 1 1223333455 

54333343334333333345554333321 11111 1111 122333345 
3333334433443344455555533332221 11111111111 133333 
23333344334444555555555444222221 111 111111111 222 
223333444 544456665555544443322221 111 1111 122 

22233444455556666655554444 3332222 1111 11112 

22223444455556666655554444333322221 1 11 1112 

222224444555566665444434443333222221 11111 1111 

12222334455555554433443343333332233321111 111 111 
1 1 1223333555554334333433333333344333331 11111 11 

1111 2333445543333333323333333555544333221 1111 1 

1111 123344433333332 2323322346666555333322 lilt 
1 1 1223433433333323323222566666555543322 1111 
1 1222333333333323321 23456666655554433221 1 1 
1 1 33333333333323212444555666555544332211 1 1 

122222322232222221 ! 1233334445 5SS4555S444333321 1 1 
22222222222121 111111 1233344444554455555543333332 
33322221222111111 1 1 1 *222234444455445555555533443 
33322221111 1 1111222233444455445555555554444 


Figure 2.2-29 Orbital Visibility, 30 Deg Inclination. 

6 Satellites 
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SHUTTLE TIME ON ORBIT {hr) 


THE ANALYTIC SCIENCES CORPORATION 
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0 

r 


OFFSET TIME (hr) 

8 12 16 




T 


20 


24 

n 


X 


T 


T 


T 


333333333221 111 
433333333322222 
443333333322221 

444333333222211 2222212221111111 

444443333222111 222222222221 1 221 
444444332211111 12222223322212111 
4444444221 1 1 112 12222223322331 1 1 1 
444444331111 12212222233333322111 
444443321 1 1 1 1222222223334332221 1 
3444332221 1 122222222334444222221 
3343222221 1 122222223344443322222 
331222222212222222334444333332211 
1111 2222222222333444444433333321 1 1 
1 1222223333444444444333333331 111 
l 11222233334444444443333333321 111 
l 1 1 1 22333344444444433333333222 1111 
L 1 1 1 1 1 23334444444443333333 32222 1111 
L 11 1 122223444443332223333322222322222 
till 1222234443222222223332222 233322222 1 
L 1 1 1222233442222221222332222344433222221 
1 1 1 22233343222221 1 1 222222 233444433322221 
1 222333322222221 11221 12333444443333222 1 
2233332222222 1111211 133334444443332222 
1 1 333322122222 1111 1 42233344444433332221 
122333222122221 1 l 1 22223344444433332222 1 
1222233222122221 1 1 22222444444433333222 1 1 


1 111 111 1 112222244444444 
111 11111 1111111222222333444 
11111121111 111 1 111112222233344 

11122233344 
11112223344 
1111222344 
111 1222244 
111122234 
111122233 
11122233 
11112223 
1112223 
11111 
111 
11 
1 


1 1222222222111221 
332222221222111 1 
332222221 11211 
33322221111 
33332221111 
3333321111 
3333332111 
333333222 

3333322221 122232221 1 122221 

433332221 1 1 1 2222332221 2221 11 

443322211 1 1 1 122222332222221 1 1 
4443221 11 11212222233322332111 
444321 1 1 1 122 1222223333333 1111 
4443211 1 1 222222222333432221 1 1 
443222 1 1 1 2222222223344222221 
4322221 1222332222344433222221 
2222222222233322344433332221 1 1 
1222223322333333444333333321 111 
1 1 222233333333444433333333321 1 1 1 
1 1 122333344444444433333333222 1111 
1 1 1 1233334444444433333333322221 1 1 1 
till 1333444444444333333332222221 1 1 1 
1 1 1 12223444444444 122233332222222 1 2221 1 
1 1 1 122234444443221 1 2223322222233322221 1 
1 1 1 122233444432221 1 12222222233344322222 
1 1 12223333433222211 122222223334443322222 
1 1222333333222221 1 122122233344443332222 
122333323222221 1 1 121 11333344444333 32222 
123333221222221 1 1 1 1 1 1 233334444433333222 
123332221222221 11 122233444444333332222 

222322211122221 22222444444433333322211 

122221221 111 1121 1 1 22222334444433333332221 1 

222211 1 21 1111111 111111 223333444433333333222221 2 

2221 111 11111111111111111 2333344433333333 3322332 
2221111 111 1111 1122333344333333333333333 


1 1 1222233444444333322211 1 
1 1 1 1 1 122333444444333322221 
1111111 3333344444443333333 
11111111 2333344444444433333 
11111111 2233334444444443333 
1111111 1 22223333444444444333 
1 1 1 1 1 1 1 122222322223444444433 
111221111122 222 1 12222224 4444443 

2222224444444 
1122223444444 
1112222444444 
1111222244444 
11 1222244444 
111122234444 
111122223344 
111112223334 
111111222222 
1111 111111 
111111 
11 11 
111 
11 
11 
1 


2-30 Orbital Visibility, 45 Deg Inclination* 
4 Satellites 


2-45 




4445444 44332222 1 1 1 1 1 1 11 1 1 2211 2222223333355555355 
544554444433333111 1 1222 1222222222222333333444555 
554555544433332 1 1 122222232223222221 1 1 12333344455 
$56555555333322 1 3333323332333221 1 1 133344455 

556665555433221 2233333333433431 11113334455 

S56666S54432111 12233334434434321 1112333455 

55666664433 1112 12233334444553322 111 1333355 

556666553222122122333445555442221 1 11 1 133345 

56666554 2222 122222333455653443221 1 1 11 1 133344 

466655453222222222234566664443321 111 11 123344 

456544333222222222245666655433331 111 111 12234 

453433333323322222456666555543322 1111 11 12224 

2222333333333333455666665554443222 4111 11111 

122223333344445555566665555444422221 111 111 

1 1222233334444555555666555544443222211 1 1 1 1 

1 1 122223344445555555554445444433322221 1111111111 
1 1 1 12222234445555555354444444433332222 1111111111 
1111 1222333345555544433343344333334333331 1111111 
1 11 11 1 23333455543333333333333333555443332111 11 1 
111 123334456433333233233223346665544333211 111 
11 122344565433332222222222456666555333321 11 1 
122234554444333222221 1233466666555533321 1 1 
2234554444433222121 1 1 3344566666555333 o l 1 1 
1 13445443443332211 1 1224445566665555333211 
12244434434433221 1 2223445556665555333321 

1233344344344332 1 1 2233355555665555533322 

1 123333333333343 1 1 1 1333344555566555 5443222 

3333333323332232 1 1 1 12233444555666555544332 

443333332223221 1 111 12224444455566665555544 

444333322221 1 1 1 1 1 1 122223444455556666655554 

444433322221 111 1111 2222334444SS5S56666S5S5 

544443222211 11 2221222233334444555555555555 

55444432221 1 1 1 1 12222322233333433334555555555 

555444333111 1 2223323332333332 1 22233355553565 

5555533332 1 1233433333343332 2233335555566 

655554 33221 1 23334434434432 11 11 23334555566 

665544422 1 1 11223334434444431 1 1 1 13333555566 

666544331 1 12 12223344444554322 1 1 1 1333355566 

666543321 l 22 1 2233344455553322 l 111333355566 

66654322222222233344565444222 11 1 1 1 133345566 

66543322222222233345664443321 111 111133334456 

6533332233233223346 6655433332 1111 11111 3334446 

3333333333333333566655553332221 111 11111 1233333 

23333344334444556665555444322221 111 11111111222 

223333444444446666555554444322221 111 1 111 122 

2223344445555666665555444433322221 111 11112 

2222344445555666655555444433332222 1 111 1112 

222224445555666665444333433333322221 1111 111 

122233345555666652333333333333332333221 111 111 

1 1 123334555666533223323323333345544332211 1 1 1 1 

1 1 1 13334456665433222222222334456654433311 1 1 1 

1 1 122344456554333222222222245566655433331 1 11 
1 122244555544333221 221222355666655543331 1 1 1 
1 2233554544333222 12 1 1 1 33446666655553333 1 11 
1234454434433322 1 l 1 1 1234445666655554333 1111 


12444344344333221 122344555666555543333111 

233433333333332 223335555556555554333221 

2333323323332232 1 1 233334455555555555433322 1 1 

333322232223221 1 1111 123344445555555555553333323 
333222212221 11111111 *222234444555455555555533443 
33322221111 1 1 1 1 1222233444455445555555554444 


Figure 2.2-31 Orbital Visibilitjs 45 Deg Inclination, 

6 Satellites 
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333333333221 III 

4333333333 2222 

44333333322221 

44433333322211 

44443333222111 

44444332221111 

44444422211111 


11111111 112222244444444 
11 111111111111222222334444 
11112211111122211112222233444 

11222223444 
1122222444 
112222444 
112222344 
1 1 222334 


2222222211122221 
122223322212221 1 1 
1 22223333223221 1 1 1 
1 1 222233332333 1111 
4444433211 1 1111 112223333343321111 
4444433 1 1 1 1 1 1 1 1 1 1 22 333334433 2 2111 11 222333 

3444322211 1 1121122233344444222211 1222333 

334322221111121122333444443322221 1 122333 

3222222221 1 1221 123334444433332221 122233 

21 1 1222222 1 2221 23334444433333322 11 11 222 

1111 12222223223334444444333333321 11 111 

11111 1 222233334444444443333333321 111 11 

1 1 1 1 1 1 122333344444444433333333222 till 1 

111111111 33334444444443333333222221 I 1 1 
1 1 1 1 1 1 122234444444332223333222222222222 
I 1 1 1 1 1 2222344444 22 2 1 12223322222234322222 1 
11111 222334444322 2111 222222 223334442222 2 
1 1 12223334443222 M 1 1 22 1 1222333444432222 1 
122233334433222 1111111 122333344443322221 


122333334332221 11 11 
113333332332221111 
1122333322322211 11 
1 122233332232221 1 1 
1 1 122222332222221 1 1 

2222222s r 1222111 

33222224 1121 
333222221111 
33332221111 
3333321111 
333333211 
333332221 
333322221 
4333222111 
4432221111 
4442211111 


1 1 1 133334444443222211 
1 122333444444332221 1 1 
1222233444444332221 4 1 1 
12222244444444332221 11 
1 11222234444444332221 111 
111112 233444444443 3 222 1 1 1 
1111111 2333444444443333333 
1 11 11112333344444444433333 
11111111 2233334444444443333 
1111111 122223333444444444333 
1111 11 11222222332223444444433 
2322211 12222211 12222344444443 


1123332221222221 1 
11223333223222211 1 
1 12223333243221 111 
11222333333333111 11 
443321 1111111 12223333343321 1 1 1 
443311 1 1 1 1 21 1222233334432221 1 1 
433221 111121 12223333444222221 
3222221 1122222223344443322222 
22222221222332233444433332221 1 
12222223222333344444333333211 11 
1 1222233323334444443333333321 111 
1 1 12233334444444443333333322211 1 1 
1 I 1 1233334444444433333333222221 1 1 1 
111 12333444444444333333332222221 111 
1 1 1 122344444444421 1 1 22332222222322221 1 
1 1 122224444444331 1 1 122222222233343222 1 1 
11 122233444443321 11 U 21 2222233344332221 1 
11222333444433221 111121 1222333344432222 1 1 
1 2233333443322 11111111 12233334444422221 1 
1 2233333433222 11111 1111 3333444444322221 1 
11333333232222111 1 1 1 123333444444332221 1 

ill 1 .iH!i?WJ}»3ii3iiiil. 

2222232221 1 1 221 111 222234*444433333222 til I 

222221221 11111 111 1 1223334444433333322221 222 

22221111111111 11111111 23333444433333333223332 

2221111 1 1111111122333344333333333333333 




1222234444444 
222223444444 
122223344444 
12223344444 
12222334444 
11222333444 
1 1 122233344 
1 1122233334 
111222233 
11 1111 11 
1 11111 
111 11 
111 
11 
1 


2-32 Orbital Visibility, 60 Deg Inclination, 
4 Satellites 


2-47 




4445444443322221 1111111 11221 1 222222333335S8SS55S 
544564444433331111 1 1 122 1222222222222333333445555 
55555554433332 11111 2222 332233223332 111 2333344 SS5 
55655555433322: 2333333333343331 11233334555 

56665555433221 12333443443444211 1123333555 

56666554442211 122334445445443311 112333555 

56666644432111 1 12223445545553321 112233455 

66666554331 11111 12224455568543322 1 1233445 

666665533221 1111 122344556655442221 1 1223444 

866654442221 121 122234566666444322 1 1 1222444 

556544332221 121 122345666665543332111 1 122344 

544443333222221 1234466666555533321111 122235 

433333333323321 23445666655554433221 1 11 11 223 

23322333333433444555666655544443222111 112 

212222333344445555855565555444432222111 1 12 

1 1 122223344445555555554445444433322221 1111111111 
111 12222244445555555554443444333332222111 111111 
1 1 1 1222333455555354433343333333333333331 111111 
1 1 111233334556664332233233233333465443332111 11 
11111 33344566654332222222223345666644333 11111 
1 1 122344566654432221221 122245566665443321 1 1 
1222345566554432221 11 1 1223355666655443321 1 
1223355565544422211 1 111334466666654433221 
11 334555455444221 1 1 12244456666655443222 

1122 444544544432 11 1 222 344556666554 432221 

1 1223444544544421 1 1223335556666655443221 

1 1 12233344344444311 1 1 133334555666655444221 1 

22233333333343321 111 12334456S56666S544421 1 

34333333222321 11 111 12223444555566665555533 

4443333322221 1 11 1 1 122223444455556666665554 

444433 322221 111 1111 22223344445555566665555 

54444322221 111 1221222233334444555555566555 

5544443221 1 1 1 22223322333333443334655556655 

555443332111 343333333333322 12233455555665 

55553333211 11244434434433321 1223345555666 

6555433222 112344454454443211 222 334555666 

6654442221 1 12234445455443311 122334465666 

6664433211 1122244455555533211 12234455666 

665543311111 112234455565543322 12233445666 

665533221121 1222344556654442221 1 1223444666 

65543222112112224455666444332111 11123344566 

5443332222222222456666554333311 11 11 122344456 

4333333233233224566665554333221 111 1 11 2233443 

23333334333433466666555544322221 111 111111 222 

223333444344456666655554444322221 1 1 1 1 111 122 

222^344445555666665555444433322221 1 1 1 111112 

2222344445555666655555444333332222 1111 1112 

12223444555566666554433343333332222111 1 1 111 

122233455556666643222233233333343443221 11 1 1 

1 12333355566665532222222223334456544322111 11 

11 133344566665542221121222234456685443221 1 1 

1 1223444666655442221121 1222445566654433221 1 
122344556655442221 11 11 122335566666444322 1 1 
122335556554442221 1 1 1 1 1 33446666665443323 1 
113344554544432211 t 1 12344486666655433221 
12244454454443221 2223445566666555433221 

1 22334 43443443321 1 2 233355 566665556433222 

233334333333432 1 11 3333455556655555433222 1 

33333233223321 1 1 122334445555655555543332333 

3333222222221 1 till 122234444555555555555434443 
3332222 1111 111 1 1222233444455445555555554444 


Figure 2.2-33 Orbital Visibility, 60 Deg Inclination, 

6 Satellites 
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333333333221 III 
433333333321 II 
443333333221 
444333332221 
444333322211 
4444332221 II 
444442221111 
444432211111 
444333111111 
343332111111 
323322211111 
222222221 111 


111111111 12222244444444 
1 1 11 1 11 12222222223444444 
11221111 1222222122222344444 


1 22332221 1 22222 1 1 1 222334 444 

1 12333322322221 111 1 122333444 

1 1 13333333332221 111 122233344 

1 133333344332211 111 112233334 

1 1233333444333111111 1 1233334 

1 1 23333444443221 111 11 233333 

11233334444442222111 1 1233333 

1 122333444444332221 11 1 133333 

1122233444444333322211 123333 

2211222221111122233444444333332222 13332 

1111 12222222222233444444333333322 1 1 221 

1111 1 1222233322444444443333333331 111 11 

1 1 1 1 1 1 1 22333344444444433333333222 1111 1 

1111111 1333344444444433333333 22222 1111 1 

1 1 1111 1223344444444422223333222222111221 
1 111122224444444431111222322222233232211 1 
111222234444444321111 122122222333333221 1 1 
12222334444443311 111111 122223333443321 111 
1 22233344444332 llllll 11 12233333443322 1111 


12233334444332111111 
1 1 1333334444332 11111 
1 1 12333333443321 1111 
1 122233333433321 1 1 1 
212222233333332111 1 
232222223322321 l I 
3322222221111 1 
*33222221 111 
3 332221111 


333332111 1 
333332111 
33332221 
3332222 
43222211 
432221 l 1 
4322111 1 
333111 111 
3322111 11 
332211111 


111 33333444433 2 11111 
12233334444332 11111 
22233334444432 llllll 
1 1 2 2233444444 33 2 1 1 1 1 1 1 
l 1 122223444444432 llllll 
1 1 1 1 1 222344444444321 1 112 
111 111 1233444444444321223 
1111111 1333344444444433333 
11111111 223333444444 4443333 
1111111 122223333444444444333 
1 1 111 1 1 12222222222233444444433 
1233322122222211 121123334444443 
1 23333323222221 11111 113333444444 
1233333333322211111 1112333344444 


122333334433221 1111 
1 1 2333334443331 1 # 11 

111 23333444443221 1 1 1 
1 1233333444442221 111 
1 1 1333334444432222 1 1 
322221 llllll 123333444443332222 
2222221 11221233334444433332221 
12222222222223344444433333321 11 
1 1222223333444444444333333321 111 
1 1 1222 33344444444433333333222 1111 
1 1 12333344444444333333333222221 111 
11 1233344444444333223333222222212211 
1 122234444444433221 1 1 1222222223222221 1 1 
1 22223444444 44322 1111111 22222333223321 1 1 


12223344444444221111 1 
222233344444422 11111 
12233334444432211 111 
113333334443321 1111 
123333334433311 1 1 1 
122333333333211 1* 
2222333323321111 
222223322222111 
222222222111 
222211111111 


112333334444 
111233334444 
11233333444 
1233333344 
1223333343 
222333323 
1 223322 
1111111 
llllll 
11111 
111 
11 
11 
1 


2221111 


1 1 2223333333321 111 
1 12233333443331 1 11 1 
1 1 13333334443321111 1 
1 1 2333334444332 11111 
22233344444432211 1 11 
122223344444442211 1 11 
1 122223444444443221 1111 
11 112223444444443322 1111212 
1 1 1 1 1 2233344444433332 2211 322 
1111111 1333344444433333333333 2 


111111111 22333344333333333333333 


Figure 2.2-34 Orbital Visibility, 90 Deg Inclination, 

4 Satellites 
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545544422111 1 12234456666655443222 1133445 

444444332211 1 122244SS6666SSSS333221 123445 
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65444222 1 22335556655443221 1 112334456666 

65443221 1 12334SS666SS53331 1 1 111234456666 

555333111 11123345666665443221 112334SS666 
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554422211 111334456666654433221 12233S556S 
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Figure 2.2-35 Orbital Visibility, 90 Deg Inclination, 

6 Satellites 
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axes are offset time, spanning a complete 24 hrs. The numbers 
in the figures show how many satellites are visible for the 
particular time on orbit and offset time. Zeros, indicating 
no satellites visible, have been replaced by blanks. Visibil- 
ity patterns observed by the Orbiter are given in the vertical 
columns, a different pattern for each offset time. 

The only criterion for visibility used in this study 
was that to be visible, a satellite must have been above 5 deg 
above the Orbiter 's horizon. (At 150 nm of altitude, the 
Orbiter' s horizon is abou" 16.5 deg below the local horizontal 
plane. ) No limitations of the Shuttle GPS receiving antennas 
were assumed, thus removing antenna considerations and Orbiter 
attitude considerations from the study. Actual visibility 
with less than full coverage antennas will often be reduced 
from the values given in the figures, but by adjusting the 
attitude of the Orbiter, if necessary, all the satellites in- 
dicated in the figures may be seen by the Orbiter antennas. 

Tables 2.2-11 and 2.2-12 summarize some statistics 
relating to the orbital visibility study. They show percent- 
ages of the time that a given number of satellites are vis- 
ible. The percentages represent averages over both time on 
orbit and over offset time. They may not apply on any par- 
ticular orbit, but the doubly periodic nature of the visi- 
bility density figures implies that they should be nearly 
correct for any orbit at the given inclination. The time 
percentage visibility figures are generally smoothly varying 
with orbital inclination, except that the 0 deg orbit shows 
some sharp changes that may be real or may be an accident of 
the discrete and limited nature of the simulations. These 
summary tables show that three or more satellites are avail- 
able significantly more often than are four or more, especially 
with the reduced, four satellite constellations. 
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TABLE 2.2-11 

TIME AVERAGE ORBITAL VISIBILITY, 4 SATELLITES 


TIME PERCENTAGES FOR ORBITAL GPS SATELLITE VISIBILITY 

(CONSTELLATION OF 4 PHASE I GPS SATELLITES, 5 deg MASK ANGLE) 
(SHUTTLE ORBITS ARE CIRCULAR AT 150 nm ALTITUDE) 

SHUTTLE 

ORBITAL 

INCLINATION 

(deg) 

NUMBER OF SATELLITES VISIBLE 

■1 

fl 

2 

3 

■ 

5 

6 

3 OR 
MORE 

4 OR 
MORE 

0 

14 

31 

13 

31 

10 



41 

10 

30 

15 

25 

27 

21 

12 



33 

12 

45 

16 

22 

28 

20 

13 



33 

13 

60 

17 

23 

25 

20 

14 



35 

14 

90 

17 

25 

20 

23 

15 



38 

15 


TABLE 2.2-12 

TIME AVERAGE ORBITAL VISIBILITY, 6 SATELLITES 


TIME PERCENTAGES FOR ORBITAL GPS SATELLITE VISIBILITY 

(CONSTELLATION OF 6 PHASE I GPS SATELLITES, 5 deg MASK ANGLE) 
(SHUTTLE ORBITS ARE CIRCULAR AT 150 nm ALTITUDE) 


SHUTTLE 

ORBITAL 

INCLINATION 

(deg) 

NUMBER OF SATELLITES VISIBLE 

0 

1 

2 

3 

■ 

5 

6 

3 OR 
MORE 

4 OR 
MORE 

0 

4 

20 

21 

13 


B 

2 

55 


30 

7 

17 

17 

27 


mm 

4 

59 


45 

9 

16 

17 

22 

14 

15 

6 

58 


60 


16 

mm 

18 

15 

15 

8 

56 

39 

90 


16 

n 

17 

14 

15 

10 

! 56 

40 
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2.3 DEORBIT NAVIGATION 
2.3.1 Introduction 

The performance levels to be expected in using Phase 
I GPS for deorbit navigation are given in this section. It 
first examines the issues of GPS satellite uploading and clock 
quality, the receiver measurement schedule, the receiver navi- 
gation filter, and the "truth model" used to assess filter 
performance. Nominal results are then presented, followed by 
sensitivity studies considering effects on performance due to 
the duration of blackout, to having IMU platform misalignment 
states in the navigation filter, and to adding drag update and 
baro-altimeter measurements. Conclusions are summarized in 
Subsection 2.3.11. 

The deorbit visibility results of Section 2.2.2 are 
used here as input to a navigation performance study of the 
deorbit mission phase. The Phase I GPS satellite constella- 
tion and the Shuttle deorbit trajectory have already been pre- 
sented along with the visibility results. The method used to 
generate performance predictions is one of covariance analy- 
sis. A reasonable navigation filter consistent with flight 
software limitations is postulated, and its statistics are 
simulated. A much more comprehensive "truth model" is then 
statistically simulated to provide an assessment of filter 
performance considering all important error sources, including 
many not acknowledged by the filter itself. Finally, key 
assumptions (including filter design) are varied to assess 
their impact on the results in a sensitivity study. A similar 
approach has been used in the past for studies of the integra- 
ted landing navigation system of Refs. 1 through 4. 
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2.3.2 GPS Satellite Uploading 

Four ground cracking stations track all of the GPS 
satellites, and current best estimates of the satellite ephem- 
erides and clock states are maintained on the ground. As each 
satellite makes its closest approach to Vandenberg AF6 (once 
per day) the ground estimate of ephemeris and clock errors for 
that satellite are uploaded to the satellite. For the next 24 
hours, the satellite then broadcasts this estimate as part of 
its GPS signal. 

Because the Phase I satellites are intentionally 
bunched, they are all uploaded during a five hour period each 
day. Since the satellite ground tracks move from west to east 
and the deorbit trajectory is west of Vandenberg, the period 
of good visibility selected in Section 2.2.2 takes place prior 
to the five hour uploading interval (with one exception) , while 
the period of reduced visibility takes place after the upload- 
ing interval. 

2.3.3 Phase I Satellite Clock Quality 

Depending primarily upon the quality of the satellite 
clocks, the GPS measurements may degrade rapidly with time 
since uploading. Phase I GPS will utilize rubidium clocks, 
while in Phase III more accurate cesium clocks are planned. 

The specified accuracy of the Phase I rubidium clocks, the 
accuracy of a 5-state model used in this study to simulate 
them, and some current data from laboratory tests of the ru- 
bidium clocks are all presented in Fig. 2.3-1 for comparison. 

As the figure shows, the test data is significantly better 
than either the specification or the assumed math model. 
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R -27409 


PHASE I RUBIDIUM CLOCK 
SPECIFICATION vs MODEL vs DATA 



Figure 2.3-1 Phase I Rubidium Clock Model vs 

Specification and Test Data 

From the mathematically modeled clock performance it 
is possible to extract a good approximation to the satellite 
clock contribution to pseudo range measurement error. A graph 
of that contribution for a typical satellite, as a function of 
the time since uploading, is given in Fig. 2.3-2. It should 
be noted that the value of about 130 ft of error at the end of 
one day since uploading is not all due to clock inaccuracy 
accumulated during the day. Part of the error is due to an 
inability of the ground network to accurately calibrate the 
clock frequency at the beginning of the day. This inability, 
itself, is partly related to clock quality, and accounts for 
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0*27409 

SATELLITE CLOCK CONTRIBUTION 
TO RANGE ERROR 



Figure 2.3-2 Satellite Clock Contribution 

to Range Error 

about as much of the error as does the clock error buildup 
since uploading. Also included in Fig. 2.3-2 is a similar 
graph applying to a 3- state model used to represent the Phase 
111 cesium clocks. The current rubidium clock test data is 
closer to the performance of the modeled cesium clocks than it 
is to that of the modeled rubidium clocks. 

With the modeled clock performance, the accuracy of the 
GPS measurements varies considerably with time since uploading, 
and the two segments of satellite visibility postulated earlier, 


^Ground test data circa 1978. 
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coming just before (good visibility) and just after (reduced 
visibility) the uploading interval, represent the extremes of 
such accuracy variation. 

2.3.4 GPS Measurement Schedule 

A GPS receiver that makes sequential pseudo range and 
pseudo range-rate measurements to single GPS satellites is 
assumed in this study. The receiver makes measurements at a 
constant rate (one every 12 sec) by cycling among the set of 
visible satellites. Measurement intervals ranging between one 
sec and one min are discussed in planned receiver designs, but 
for a user with an inertial system of the quality of that planned 
for the Orbiter, the measureme. i: rate makes very little differ- 
ence within these limits, provided that initial transient error 
reduction using GPS (for the first 8 measurements or so) is 
considered on a "number-of-measurements” basis and not on a 
time basis. That is, after the fifth measurement, say, re- 
sults are likely to be very similar whether the five measure- 
ments took five sec or five min. Considered on a time basis, 
a higher measurement rate is always better. 

2.^.5 Deorbit Navigation Filter 

The filter chosen for GPS-aided deorbit navigation is 
an 11-state Kalman filter containing 3 positions states, 3 
velocity states, 3 platform alignment states (all in inertial 
coordinates), and 2 clock states (phase and frequency). Table 
2.3-1 summarizes the initial state standard deviations and the 
process noise densities for the filter. 

Measurement noise values for GPS measurements are 
shown in Table 2.3-2. Note that there are two sets of values, 
one appropriate for use prior to updating and another set with 
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TABLE 2.3-1 

DEORBIT NAVIGATION FILTER INITIAL CONDITIONS 
AND PROCESS NOISES 


T-1138 




AXIS 



DR* 

CR* 

R 

INITIAL 

Position (ft) 

WmjTjjM 

2500 

2500 

RMS 

Velocity (fps) 


3.3 

26.3+ 

ERROR 

Alignment (sec) 

mm 

729 

729 

PROCESS 

Position (ft2/sec) 

— - 

— 

— 

NOISE 

SPECTRAL 

Velocity (fps^/sec) 

0.1076 

0.1076 

0.1076 

DENSITY 

Alignment (sec^/sec) 

4.41 

4.41 

4.41 

INITIAL 

RMS 

Phase ( f t ) 

1200^ 

CLOCK 

ERROR 

Frequency (fps) 

0.25 * 

CLOCK 

PROCESS 

NOISE 

Phase (ft^/sec) 

0.1 

Frequency (fps 2 /sec) 

3 . OxlO” 5 

SPECTRAL 

DENSITY 


♦Referenced to inertial velocity. 

tlnitial downrange position and radial velocity errors are 
-90% correlated 

^Initial clock phrase and frequency errors are +50% 
correlated. 


TABLE 2.3-2 

DEORBIT FILTER MEASUREMENT NOISE 


RMS 

MEASUREMENT 

ERROR 

PRIOR TO 
UPLOADING 

AFTER 

UPLOADING 

Pseudo Range (ft) 

100 

12 

Pseudo Range-Rate (fps) 

0.05 

0.05 
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a much reduced pseudo range error appropriate for use after 
updating. These variations are an attempt to account for the 
unmodeled (in the filter) effect of satellite states on pseudo 
range measurements as discussed in Section 2.3.3. 

The deorbit filter incorporates underweighting as 
discussed in Section 2.3 of Ref. 4. Underweighting reduces 
the possibility of filter divergence due to nonlinear effects 
and also reduces the sensitivity of the filter to a single bad 
measurement. It slows convergence when all is well, however. 

An underweighting factor of 1.2 is used until the root sum of 
the three position variances falls to less than 3281 ft (1000 
meters), at which point the underweighting is turned off. Due 
to the error growth during the blackout zone of deorbit, the 
underweighting is turned back on subsequently, and then, final- 
ly, turned off again if and when the position error falls back 
below the limit above. 

When drag updating and baro-altimeter aids are added 
to GPS aiding, the filter includes an additional state. This 
state represents either drag altitude error or barometric alti- 
tude error, depending on whichever aid is currently being used, 
since the time periods when the two additional aids are being 
used does not overlap. Numerical values appropriate to this 
filter are those discussed in Chapter 5 of Volume 2. 

2.3.6 Deorbit Navigation Truth Model 

The truth model used to evaluate deorbit navigation 
with GPS aiding alone is similar to that described in Section 
2.4.1 of Ref. 4. In particular, the 20-state model for the 
inertial measurement unit (IMU) and the 6-state model for the 
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master timing unit (MTU, also called the Orbiter clock) are 
identical. The gravity model error section of the earlier 
study is not exercised here due to the negligible effect of 
gravity errors in past studies at lower altitudes where such 
errors would be even larger. The satellite states in the truth 
model are changed because a 5- state rubidium clock model is 
used here for Phase 1 GPS instead of tae 3-state cesium model 
used for Phase III GPS in the earlier study, and because there 
are only 6 Phase 1 satellites considered here. Thus, there 
are a total of 78 satellite states (6 satellites x 13 states 
per satellite). Including the 11 filter states, the current 
truth model contains 115 states. 

The data base for the current truth model reflects 
the likely initial values for the filter states and for the 
MTU for deorbit navigation, assuming GPS had not been used for 
the two hours preceding deorbit. The IMU values are identical 
to those assumed in Table 4.2-2 of Volume 2. The GPS satel- 
lite state numerical values are used in the NESA program (de- 
scribed in Section 2.2 of Ref. 4) and are appropriate to Phase 
I GPS. Tabic 2.3-3 summarizes the data base for the deorbit 
navigation truth model. 

When drag updating and baro-altimeter measurements 
are taken along with the GPS measurements, the truth model 
necessarily becomes more complex. The additional models are 
described in Sections 4.3-1 and 4.3-2 of Ref. 4. The truth 
model data base for these additional states includes all drag 
and baro-related error sources as presented in Ref. 4 and used 
in generating the baseline error budgets of that report. Table 
4.3-4 and a portion of Table 4.3-5 of Ref. 4 summarize the 
added data base. 


*When this study was done, it was thought that the MTU would 
serve as the GPS receiver clock. This is no longer the case. 
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TABLE 2.3-3 

DATA BASE FOR DEORBIT NAVIGATION TRUTH MODEL 


i=im 


ERROR SOURCE 

STANDARD DEVIATION 

0. Uncor related GPS Receiver 


Measurement Errors 


Pseudo Range 

12 ft 

Pseudo Range^Rate 

0.05 fps 

1. Initial Filter State Errors 


Position 

(11265, 860, 1038) ft* 

Velocity 

(1.086, 0.942, 12. 468) fps* 

Alignment 

284 sec each of 3 axes 

2 . IMG Errors 


Accelerometer 


Biases 

50 ug each of 3 axes 

Scale Factors 

100 ppm each of 3 axes 

Asymmetries 

100 ppm each of 3 axes 

Nonort hogonal it les 

15 sec each of 3 axes 

Gyro 

Bias Drifts 

0.035 deg/hr each of 3 axes 

Mass Unbalances 

0.025 deg/hr/g each of 5 components 

3. Shuttle MTU Errors 


Initial Phase 

1130 ft f 

Initial Frequency 

0.142 fps 

Aging 

3. 80 x 10~3 ft/sec2 

Random Frequency 

20.2 fps# ^ 

Random Frequency Rate 

1.37 ft /sec** 

Acceleration Sensitivity 

0.1 ft/sec 2 /g 

4. GPS Satellite Errors 

See Text 


•Given in downrange-crossrange- radial coordinates referenced 
to Inertial velocity; significant cross correlations exist 
among initial position and velocity errors 


^Correlations among initial clock errors exist; also, clock 
process noises exist as given in Section 2.4.1 for Ref. 2 

#Clock flicker noise is attenuated by a factor of 100 before 
combining with ordinary frequency error 


2-3.7 Results of GPS Aiding Alone 

Position and velocity error results of a truth moJel 
evaluation of GPS aiding during deorbit, assuming good visibil- 
ity prior to satellite uploading, are presented in Figs. 2.3-3 
and 2.3-4. Similar results for a period of reduced visibility 
after satellite uploading are presented in Figs. 2.3-5 and 
2.3-6. The figures are largely self-explanatory, but because 
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Figure 2.3-3 RMS Position Error versus Time for GPS-Aided 

Deorbit Navigation with Good Visibility Prior 
to Satellite Uploading 



Figure 2.3-4 RMS Velocity Error versus Time for GPS-Aided 

Deorbit Navigation with Good Visibility Prior 
to Satellite Uploading 
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Figure 2.3-5 RMS Position Error versus Time for GPS-Aided 

Deorbit Navigation with Reduced Visibility 
After Satellite Uploading 



Figure 2.3-6 RMS Velocity Error versus Time for GPS-Aided 

Deorbit Navigation ith Reduced Visibility 
After Satellite Uploading 
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of Che scales necessary to present all of the results on a 
linear graph they tend to obscure some details of the tran- 
sient error reductions, both the initial reductions and the 
reductions following the blackout period. Figures 2.3-7 and 
2.3-8 show these initial reductions of position errors on a 
logarithmic scale and on a number-of -measurements basis rather 
than a time basis. It should be pointed out that during the 
initial portion of deorbit with good visibility there are only 
three satellites visible, fewer than are needed for complete 
position determination based on GPS alone. 






NUMBER OF MEASUREMENTS 
AFTER BLACKOUT 


Figure 2.3-7 Initial Error Reduction with GPS < Based on 

Good Visibility Prior to Uploading) 


Tables 2.3- / » thru 2.3-6 contain error budgets prior 
to blackout, after blackout, and at the terminal 20 kft of 
altitude points, both for the good visibility situation and 
the reduced visibility situation. After blackout, the results 
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R'JMOC 



1 2 3 4 5 6 7 8 012348678 

NUMBER OF MEASUREMENTS NUMBER OF MEASUREMENTS 

AFTER INITIALIZATION AFTER BLACKOUT 


Figure 2.3-8 Initial Error Reduction with GPS (Based on 

Reduced Visibility After Uploading) 


TABLE 2.3-4 

GPS-AIDED DEORBIT NAVIGATION ERROR BUDGETS PRIOR TO BLACKOUT 


T4301 


BASED ON 
GOOD 

VISIBILITY 
PRIOR TO 
UPLOADING 


ERROR GROUP 

RMS POSITION 
ERROR (FT) 

DR CR R 

RMS VELOCITY 
ERROR (FPS) 
DR CR R 

INITIAL CONDITIONS 

0 

0 

1 

0.01 

0.01 

0.01 

MEASUREMENT NOISE 

7 

5 

14 

0.05 

0.05 

0.06 

INERTIAL MEASUREMENT UNIT 

1 

1 

1 

0.02 

0.04 

0.04 

ORB ITER CLOCK 

8 

8 

35 

0.12 

0.14 

0.26 

SATELLITE STATES 

133 

100 

300 

0.03 

0.04 

0.42 

RSS 

133 

101 

302 

u . 14 

0.16 

0.50 


BASED ON 

REDUCED 

VISIBILITY 

AFTER 

UPLOADING 


INITIAL CONDITIONS 

539 

110 

82 

0.57 

0.17 

0.13 

MEASUREMENT NOISE 

52 

8 

14 

0.07 

0.10 

0.22 

INERTIAL MEASUREMENT UNIT 

962 

53 

48 

1.21 

0.22 

0.30 

ORB ITER CLOCK 

CiTPf 1 T TP PT.'Prf 

44 2 

1039 

652 

0.40 

0.37 

0.13 

SATELLITE STATES 






MM 

RSS 

1189 

1046 

659 

1.40 

0.47 

0.42 
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TABLE 2.3-5 

GPS -AIDED DEORBIT NAVIGATION ERROR BUDGETS AFTER BLACKOUT 


T-1023 





RMS 

POSITION 

RMS VELOCITY 


ERROR CROUP 


ERROR (FT) 

ERROR (FT) 




DR 

CR 

R 

DR 

CR 

R 


INITIAL CONDITIONS 


2976 

9088 

11915 

14.35 

24.43 

34.86 

BASED ON 
GOOD 

MEASUREMENT NOISE 


58 

42 

83 

0.05 

0.04 

0.17 

VISIBILITY 

INERTIAL MEASUREMENT 

UNIT 

1329 

1950 

2619 

4.46 

5.78 

8.15 

PRIOR TO 
UPLOADING 

ORB ITER CLOCK 


195 

132 

338 

0.19 

0.12 

0.67 


SATELLITE STATES 


410* 

210* 

627* 

0.16 

0.17 

1.62 


RSS 

3292 

9299 

12233 

15.02 

25.11 

35.84 


INITIAL CONDITIONS 


2965 

9109 

11942 

14.59 

24.49 

35.08 

BASED ON 
REDUCED 

MEASUREMENT NOISE 


154 

133 

194 

0.12 

0.11 

0.39 

VISIBILITY 

INERTIAL MEASUREMENT 

UNIT 

1475 

1968 

3654 

4.73 

5.77 

9.24 

AFTER 

UPLOADING 

ORB I TER CLOCK 


1013 

343 

1552 

0.31 

1.25 

2.40 


SATELLITE STATES 


— 

— 

— 

— 

— 

— 


RSS 

3467 

9326 

12586 

15.34 

25.19 

36.36 


•SATELLITE CLOCKS CONSISTENT SflTH RECENT TEST DATA WOULD YIELD POSITION 
ERROR CONTRIBUTIONS 10-20 TIMES SMALLER. 


are about the same in either case, as the errors are dominantly 
due to initial platform misalignments with a lesser amount due 
to the IMU. Because there is very little sensed acceleration 
prior to blackout, the initial platform alignments are not 
observable and thus cannot be estimated during this time in- 
terval from th'j GPS measurements. The IMU errors are not 
estimated by the filter, so their contribution during blackout 
cannot be reduced by GPS measurements prior to blackout. 

Platform alignment errors grow somewhat from their 
initial values up to the end of the blackout period, due to 
gyro drift contributions. After blackout, the filter begins 
to estimate the alignments and to reduce the rms alignment 
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TABLE 2.3-6 

GPS -AIDED DEORBIT NAVIGATION ERROR BUDGETS 



AT 20,000 FT 

ALTITUDE 



T-1024 

ERROR GROUP 

RMS POSITION 
ERROR (FT) 

RMS VELOCITY 
ERROR (FPS) 

DR 

CR 

R 

DR 

CR 

R 



INITIAL CONDITIONS 

1 

1 

0 

0.07 

0.06 

0.06 

GOOD 


MEASUREMENT NOISE 

11 

10 

14 

0.08 

0.07 

0.05 

VISIBILITY 


INERTIAL MEASUREMENT UNIT 

2 

2 

0 

0.07 

0.13 

0.07 

PRIOR TO 
UPLOADING 


ORB I TER CLOCK 

18 

22 

35 

0.18 

0.16 

0.24 



SATELLITE STATES 

314* 

329 

425 

0.05 

0.02 

0.30 



3SS 

314 

330 

427 

0.22 

0.23 

0.39 



INITIAL CONDITIONS 

1267 

317 

1143 

1.44 

0.36 

1.44 



MEASUREMENT NOISE 

265 

71 

171 

0.13 

0.09 

0.32 

REDUCED 

VISIBILITY 


INERTIAL MEASUREMENT UNIT 

2474 

678 

1364 

2.85 

0.68 

2.56 

AFTER 


ORBITER CLOCK 

1617 

1945 

3811 

4.01 

1.04 

3.79 

UPLOADING 


SATELLITE STATES 

— 

— 

• 

— 

— 

— 



RSS 

3226 

2085 

4209 

5.12 

1.30 

4.80 


•SATELLITE CLOCKS CONSISTENT WITH RECENT TEST DATA WOULD YIELD POSITION ERROR 
CONTRIBUTIONS 10-20 TIMES SMALLER 


error , as Table 2.3-7 shows. Sensed acceleration during and 
after blackout is largely in the downraage direction, with a 
lesser amount in the radial direction, and very little n the 
crossrange direction, so alignment observability is best in 
the crossrange direction and worst in the downrange direction, 
which the table confirms. 

2.3.8 Effect of the Duration of Blackout 

Previous Figs. 2.3-3 through 2.3-6 show the error 
growth during the blackout period as a function of time. An 
earlier emergence from blackout would cause errors to peak at 
the lower values plotted for times before 37.3 min. Time and 
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TABLE 2.3-7 

PLATFORM ALIGNMENT ERRORS FOR GPS -AIDED 
DEORBIT NAVIGATION 


T-1141 



GOOD VISIBILITY PBIOB 
TO SATELLITE UPLOADING 

REDUCED VISIBILITY AFTER 
SATELLITE UPLOADING 

EVENT 

RMS MISALIGNMENTS 

(SEC) 

RMS MISALIGNMENTS 

(SEC) 

DR 

CR 

tm 


CR 

R 

INITIAL CONDITIONS 

284 

284 

284 

284 

284 

284 

PRIOR TO FIRST MEASUREMENT 
AFTER BLACKOUT 

290 

291 

292 

292 

293 

294 

EIGHT MEASUREMENTS AFTER 
BLACKOUT 

112 

46 

77 

151 

51 

130 

AT 20 KFT ALTITUDE 

47 

38 

39 

52 

53 

39 


altitude are related in a nonlinear way during blackout, how- 
ever, so it is not directly obvious from the previous figures 
how error growth would be reduced if the altitude (nominally 
145 kft) assumed for the end of blackout were raised. 

The other way in which the blackout interval could be 
reduced, by lowering the altitude (nominally 300 kft) assumed 
for the onset of blackout, would have a less obvious effect on 
the error growth. A separate investigation is necessary to 
examine each variation in the altitude for the onset of black- 
out. Figures 2.3-9 and 2.3-10 show how errors grow during 
blackout as a function of altitude, so that the effect of end- 
ing blackout at a higher altitude is easily assessed, and how 
the errors at the end of blackout (the nominal 145 kft) reduce 
if blackout begins at altitudes lower than 300 kft. A general 
conclusion is that an early emergence from blackout would clear 
ly reduce errors, but that a delay in the onset of blackout 
would be helpful only if the delay were to an altitude below 
240 kft. 
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ft -274M 


EFFECT OF A DELAY IN THE ONSET OF 

ERROR GROWTH DURING BLACKOUT BLACKOUT ON ERRORS AT THE ENO OF BLACKOUT 



300 250 200 150 300 2S0 200 150 

ALTITUDE Ikftl ALTITUDE AT ONSET OF BLACKOUT Ikftl 


Figure 2.3-9 Sensitivity of Aided Navigation Results to 

the Duration of Blackout - - Position Errors 


ft-27412 

EFFECT OF A DELAY IN THE ONSET OF 

ERROR GROWTH DURING BLACKOUT BLACKOUT ON ERRORS AT THE END OF 8LACKOUT 



300 250 200 150 300 250 200 ISO 

ALTITUDE Ikftl ALTITUDE AT ONSET OF BLACKOUT ikftl 


Sensitivity of Aided Navigation Results to 
the Duration of Blackout -- Velocity Errors 
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2.3.9 Effect of Misalignment States in the Filter 

A possible simplification of the deorbit navigation 
filter would be to omit the three platform alignment states. 
These states cannot possibly be estimated until after black- 
out, because before blackout there is inadequate sensed 
acceleration to have platform misalignment cause any effect, 
and during blackout there are no measurements. The misalign- 
ment states are useful, however, in quickly reducing the nav- 
igation errors after blackout, as Fig. 2.3-11 shows. The 
reason for this behavior is that with alignment states pres- 
ent, the filter correctly accounts for the cross correlations 
in position and velocity error components that build up during 
blackout, making subsequent identification of errors an easier 
process . 



NUMBER OF MEASUREMENTS 
AFTER BLACKOUT 



AFTER BLACKOUT 


Figure 2.3-11 Effect of Misalignment States in the Filter 

(Based on Good Visibility Prior to Uploading) 
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2.3.10 Effect of Adding Drag Updates and a Baro- Altimeter 

One way to control the error buildup during blackout 
and to assist in recovery after blackout in a reduced visibil- 
ity situation is to add drag updating and baro- altimeter meas- 
urements. Both of these aids are useful in reducing errors in 
the plane of motion, especially in the radial direction. 

Drag updating begins around 282 kft of altitude, just 
after blackout begins, and continues past the nominal end of 
blackout down to 100 kft of altitude. One update is assumed 
each 3.6 seconds in this analysis. Baro-altimeter measure- 
ments do not begin until 31 kft of altitude is reached, and 
continue to the 20 kft point at which this study stops. 
Baro-altitude measurements are also assumed at 3.6 second 
intervals. Because these additional aids are most useful 
under reduced visibility conditions, such conditions are 
assumed to exist. Reduced visibility conditions here follow 
satellite uploading, so that questions of satellite clock qual- 
ity are obviated. 

Figure 2.3-12 shows the error growth during blackout 
when dreg updating is used, while Fig. 2.1-18 shows the init- 
ial error reduction following the# blackout period. Table 2.3-8 
gives error budgets after blackout and at 20 kft of altitude 
when drag updates and baro-altimeter measurements are used. 

Note especially the reductions in radial position and velocity 
errors when the additional aids are used in conjunction witn 
even reduced visibility GPS measurements. 


*After the init.'al series of Shuttle test fights, the air 
data system will be sufficientlv well calibrated to pernu ' 
operation at higher altitudes, possibly up to 85-100 kft. 
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HIM 

ERROR GROWTH DURING BLACKOUT 



ALTITUDE (left) ALTITUDE Iktt) 

Figure 2.3-12 GPS + Drag Update + Baro-Altimeter Aided 

Deorbit Navigation (Error Growth During 
Blackout) 

2.3.11 Deorbit Navigation Conclusions 

An important aspect of deorbit navigation concerns 
visibility coverage for the deorbit mission with Phase 1 GPS. 

It is shown in section 2.2.2 that good visibility is obtained 
only about five hours each day, with reduced visibility obtained 
for a 12-16 hour interval each day. For the remaining 8-12 
hours each day, satellite visibility during deorbit is poor to 
nonexistent. This study shows that GPS aiding alone produces 
good deorbit navigation results during periods of good visibil- 
ity. while acceptable results with reduced visibility may be 
obtained by adding drag updating and baro-altimeter aids . 
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TABLE 2.3-8 

GPS + DRAG + BARO-ALTIMETER AIDED NAVIGATION ERROR BUDGETS 
(BASED ON REDUCED VISIBILITY AFTER UPLOADING) 


2zim 



ronnn GDrtTTD 

RMS POSITION 
ERROR (FT) 

RMS VELOCITY 
ERROR (FPS) 



DR 

CR 

R 

DR 

CR 

R 


INITIAL CONDITIONS 

1366 

7897 

1661 

D 

21.98 

6.72 


GPS MEASUREMENT NOISE 

83 

180 

110 

I ' 

0.22 

0.18 

AFTER 

BLACKOUT 

INERTIAL MEASUREMENT UNIT 

621 

1890 

874 


5.57 

2.78 

ORBITER CLOCK 

124 

279 

1477 

IBB 

2.96 

1.13 

DRAG UPDATE STATES 

343 

459 

1153 

1.08 

1.24 

2.50 


BARO-ALTIMETER STATES 

— 

— 

— 

— 

— 

— 


RSS 

1547 

8140 

2654 

4.56 

22.91 

7.78 


INITIAL CONDITIONS 

358 

64 

293 


0.06 

0.28 


GPS MEASUREMENT NOISE 

126 

19 

56 

0.15 

0.08 

0.18 


INERTIAL MEASUREMENT UNIT 

1249 

188 

245 

1.65 

0.18 

1.15 

AT 

ORBITER CLOCK 

2144 

382 

401 

0.32 

0.15 

0.63 

20 , 000 FT 
ALTITUDE 

DRAG UPDATE STATES 

76 

11 

29 

0.09 

0.01 

0.07 

BARO-ALTIMETER STATES 

678 

55 

640 

0.72 

0.06 

0.77 


RSS 

2601 

434 

840 

1.88 

0.26 

1.55 


Provided the time between measurements is on the order 
of one min or less, a sequential receiver cycling among vis- 
ible GPS satellites is adequate. The navigation results are 
insensitive to the actual time between measurements, if tran- 
sient error reduction is considered on a "number-of-measure- 
ments" basis rather than on a time basis. 

Error buildup during the blackout period is severe 
with GPS aiding alone . It may be reduced either by shortening 
the blackout interval or by adding drag updating. 
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Mnw 

INITIAL ERROR REDUCTION AFTER BLACKOUT 
(based on reduced visibility after uploading) 



MEASUREMENT NUMBER AFTER BLACKOUT 


UNDERWEIGHTING 
REMAINS ON 


Figure 2.3-13 GPS + Drag Update + Baro-Altimeter Aided 

Deorbit Navigation (Post Blackout Error 
Reduction ) 


With good satellite visibility prior to satellite 
ephemeris uploading, and assuming the modeled 5-state rubidium 
satellite clock accuracy, "steady-state" errors are dominated 
by satellite ephemeris errors. With reduced visibility after 
uploading, "steady-state" errors are much larger and are dom- 
inated by the Orbiter clock and IMU. 
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Platform alignment states in the deorbit navigation 
filter are of limited utility in the early stages of deorbit, 
but help in rapid transient error reduction following blackout. 
The addition of drag updates and baro-altimeter measurements 
to the GPS measurements substantially reduces errors at low 
altitudes, especially in the radial direction. 


2.4 ORBITAL NAVIGATION 
2.4.1 Introduction 


Orbital navigation performance with Phase I GPS is 
treated in this section in a manner similar to that used for 
deorbit navigation in the previous section. There are key 
differences. No IMU information is utilized on orbit, as 
sensed acceleration levels are generally below the IMU thresh- 
old. Instead, low-level venting thrusts are the dynamical 
error drivers and must be modeled. The satellite constel- 
lation and orbital trajectory are discussed along with the 
visibility results in previous Section 2.2.3. In this section, 
the orbital navigation filter in the GPS receiver is discussed, 
as is the truth model used to evaluate the filter performance. 
The nominal results are presented, followed by a set of re- 
sults for a situation slightly modified to assess error build- 
up during periods of no satellite visibility. Finally, con- 
clusions are summarized. 

2.4.2 Orbital Navigation Filter 

There are fundamental differences between orbital and 
deorbit navigation that ca.l for changes in the navigation 
filter. The filter for orbital navigation remains an 11-state 
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Kalman filter with underweighting, but the three platform align- 
ment states used during deorbit are not useful on orbit. Vehicle 
dynamics on orbit, besides those due to the accurately predict- 
able earth gravitational field, are primarily due to venting 
effects, and these are below the threshold of the IMU. In 
place of the three platform alignments states, the orbital 
filter has three thrust acceleration states, modeled as first- 
order Markov processes in downrange-crossrange-radial coordi- 
nates. This approach is similar to that used in Section 3.3 
of Ref. 21. 


Discussions with JSC (Ref. 22) related that the Shuttle 
venting details are not yet firm, but that rms values in the 0.5 
to 1.0 lb of thrust range (1 lb of thrust is 5pg's of acceler- 
ation on a 200,000 lb vehicle) are probable. Venting thrusts 
will probably be dominately along the Orbiter body forward-aft 
and "vertical" axes with a smaller amount in the port-starboard 
direction. Expecting that the forward-aft body axis will be 
near the orbital plane more often that not, and using some 
conservation in filter design, the rms thrust values used in 
the filter are those of Table 2.4-1. Venting thrusts may be 
highly correlated in Orbiter body coordinates, but they will 
be less so in the orbital coordinates of downrange-crossrange- 
radial. The first-order Markov process time constants assumed 
for thrust in Table 2.4-1 reflect that fact. 

,2.4.3 Orbital Navigation Truth Model 

The truth model for GPS-aided orbital navigation is 
similar to that used for deorbit navigation with a few excep- 
tions. Due to the very low level of sensed acceleration on 
orbit, the IMU states are not important. Hence, the 20 IMU 
states are removed from the truth model. For the same reason, 
one of the MTU states, an acceleration sensitive clock state, 
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TABLE 2.4-1 

THRUST PROCESS PARAMETER VALUES USED IN 
THE ORBITAL NAVIGATION FILTER 


T-1152 


THRUST 

PROCESS 

PARAMETER 


AXIS 


DR - 

CR 

1 R 

Time Constant (sec) 

600 

300 

600 

RMS Initial Thrust 
(lb) 

1.5 

0.3 

1.5 

Thrust Process Noise 
Intensity (lb/ /sec) 

0.0866 

0.0408 

0.0866 


is removed, reducing Che MTU model to five states. And of 
course, the 11 filter states in the truth model correspond to 
those just discussed, incorporating the three thrust states in 
place of the three platform alignment states used for deorbit 
studies. 


The data base for the orbital navigation truth model, 
where the two models have corresponding states, is entirely 
the same as that used for deorbit and summarized in Table 2.3-3. 
The truth model values for the thrust states are shown in Table 
2.4-2. Thr thrust state time constants used in the filter are 
kept, but the initial rms state values and the corresponding 
process noises are reduced by 1/3. 

2.4.4 Orbital Navigation Results 

Position and velocity error results of a truth model 
evaluation of GPS aiding on orbit, assuming a recent satellite 
uploading, are presented in Figs. 2.4-1 and 2.4-2. As the 
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TABLE 2.4-2 

THRUST PROCESS PARAMETER VALUES USED IN 
THE ORBITAL NAVIGATION TRUTH MODEL 


T-1156 


THRUST 

DDAprcc 

AXIS 


rnUttoo 

PARAMETER 

DR 

CR 

R 


Time Constant (sec) 

600 

300 

600 


RMS Initial Thrust 
(lb) 

1.0 

0.33 

1.0 



NO MEASUREMENTS 



SEQUENTIAL RECEIVER 
20 ac MEASUREMENT OKIE 


1210 



0 <8 70 30 40 S8COns090 

E LAPSE 0 TIME Mart 

* REFEREHCED TO IRERTIAL VELOCITY 


Figure 2.4-1 Orbital Navigation: RMS Position Error 

Versus Time 
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Figure 2.4-2 Orbital Navigation: RMS Velocity Error 

Versus Time 

figures show, initial errors are reduced during the first 10 
min segment of fading satellite visibility. The errors grow 
during the next 20 min segment of no visibility, and then they 
fall again as visibility improves. An error budget at 89 min 
is presented in Table 2.4-3. After more than one orbit, er- 
rors prior to a segment of no visibility would be smaller than 
those at 10 min of Figs. 2.4-1 and 2.4-2. In the figures, the 
effect of the large initial errors was not yet overcome before 
the 20 min segment of no visibility was encountered. 

2.4.5 Results for a Modified Measurement Schedule 


To see how large orbital navigation errors would be- 
come during a typical period of no satellite visibility after 
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TABLE 2.4-3 

ORBITAL NAVIGATION ERROR BUDGET 
(89 MINUTES) 


ERROR SOURCE 
GROUP 

POSITION ERROR 
(FT) 

1 ICSI 1 

CLOCK ERROR 

DR 

CR 

R 

DR 

CR 

R 

PHASE 

(FT) 

FREQUENCY 

(FPS) 

INITIAL CONDITIONS 

0 


n 

0 

0 

0 

0 

0 

MEASUREMENT NOISE 

21 

16 

El 

.04 

.02 

.04 

21 

.04 

THRUST PROCESS 

7 

B 


.03 

.02 

.03 

8 

.02 

CLOCK 

27 

Q 


.09 

.03 

.10 

23 

.08 

SATELLITE ERRORS 

59 

D 

D 

.02 

.02 

.04 

62 

.02 

RSS 

66 

45 

45 

.10 

.04 

.12 

69 

.09 

VECTOR RSS 

93 

0.16 

i 


first stabilizing during a period of good visibility, the or- 
bital navigation results were redone after adjusting the pre- 
vious measurement schedule so it simulates such a situation. 
Figure 2.4-3 shows the adjusted visibility pattern and meas- 
urement schedule shifted from that in Section 2.2 "end -around” 
by 50 min. 

Position and velocity error results of a truth model 
evaluation of GPS aiding for this visibility situation on or- 
bit, assuming, as before, recent satellite uploadings, are 
presented in Figs. 2.4-4 and 2.4-5. The figures show that 
position errors are reduced to reasonably stable low values 
after about 10 min, while velocity errors require about 20 
min. Note from Fig. 2.4-3 that 4 satellites do not become 
visible until about 10 min of elapsed time, largely explaining 
the initial error behavior. 
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Figure 2.4-3 Revised Orbital Satellite Visibility and 

Measurement Schedule 

Downrange position error and radial velocity error grow 
most rapidly during periods of no visibility. In the results 
previously reported, rms downrange position error grew to 1200 
ft and rms radial velocity error grew to 1.82 fps during the 
20 min period of no visibility. Here, as Figs. 2.4-4 and 2.4-5 
show, rms downrange position error grows to 870 ft while rms 
radial velocity error grows to 1.34 fps during the 20 min pe- 
riod of no visibility that occurs after GPS aiding has had time 
to overcome the effect of the large initial navigation errors. 
The present results are expected to be typical of error growth 
during periods of no visibility on orbit using Phase I GPS . 

2.4.6 Orbital Navigation Conclusions 

GPS appears particularly well suited to Shuttle orbital 
navigation. Even Phase I GPS with its limited coverage should 
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A-MMt 



ELAPSED TIME (min) 


Figure 2.4-4 Orbital Navigation: RMS Position Error 

Versus Time 

provide good satellite visibility over some segment of almost 
every orbit, whatever the orbital inclination. During periods 
of good visibility subsequent to a satellite uploading inter- 
val , rms orbital navigation errors should reduce to less than 
100 ft in position and 0.2 fps in velocity . After 20 min of 
no visibility following a period of good visibility, rms nav- 
igation errors should grow into the 600-900 ft range in posi- 
tion and be less than 1.5 fps in velocity . 

An additional conclusion evident in the detailed exam- 
ination of the truth model results is that while the filter 
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Figure 2.4-5 Orbital Navigation: RMS Velocity Error 

Versus Time 

does not do well in estimating the thrust states, it uses the 
thrust states to compensate for correlated frequency noise in 
the Orbiter MTU, making for a slight improvement in velocity 
estimation . 
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SINGLE SATELLITE GPS ORBITAL NAVIGATION 


3 . 1 INTRODUCTION 

This chapter presents results of a study of orbital 
Space Shuttle navigation using a single GPS satellite. The 
study meets several practical objectives. It demonstrates the 
minimum capability of GPS for Shuttle navigation in a worst 
case situation. While the current schedule for GPS satellite 
launches calls for several satellites to be available in the 
time frame projected for GPS user equipment implementation 
aboard the Shuttle, schedule changes or satellite failures 
could conceivably result in only one or two operational satel- 
lites during some Shuttle flights with a GPS receiver. This 
study provides a valuable preliminary examination necessary 
for guiding decisions on the use of GPS for orbital navigation 
with an incomplete satellite system. Even if several satel- 
lites are always available, this study still is useful for the 
insights it provides. The results from a single satellite 
study bring out fundamental relationships in the satellite 
navigation system, thus leading to a better understanding of 
the complex interactions of the multisatellite navigation 
system. 


3 . 2 METHODOLOGY 

In all of the results reported here, the single GPS 
satellite was assumed to be in a typical GPS circular orbit of 
63 deg inclination at a radius of approximately 4.164 earth 
radii. The Shuttle was also assumed to be in a circular orbit, 
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always at 150 nni altitude (1.044 earth radii), sometimes at a 
38 deg inclination typical of some early missions, and some* 
times at a 63 deg inclination, matching that of the GPS sat- 
ellite. to consider extreme Shuttle-GPS geometry. Right 
ascension and mean anomaly of the Shuttle orbit relative to 
the GPS satellite orbit were varied. 

A preparatory satellite visibility study indicated 
that typical visibility periods of 45 minutes and gaps of 55 
minutes could be expected. Variations in initial mean anom- 
alies made little difference when 12 hour periods were con- 
sidered. Variations in right ascensions , however, could cause 
gaps as large as 104 minutes and visibility periods as short 
as 6 minutes in extreme cases. This preliminary analysis 
helped establish the reasonable Shuttle-GPS satellite orbital 
relationships for the navigation performance study. 

An 8-state suboptimal GPS navigation filter and a 14- 
state truth model are used for error analysis. The states of 
the filter and of the truth model are presented in Table 3.2-1 
Detailed development of the models and state variables is con- 
tained in Refs. 3 and 4, and in Chapter 2. A 14-state truth 
model is sufficient for the present analysis because only one 
satellite is considered; some small error contributing factors 
are not modeled. 

Many simplifications are possible in this single sat- 
ellite study since simulation of a particular mission is not 
being attempted. Circular orbits are assumed with no maneu- 
vering on orbit, and no gravity errors per se are included. 
Satellite ephemeris and clock errors are also absent, as they 
are minor contributors (Ref. 4) and are reasonably difficult 
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TABLE 3.2-1 

FILTER AND TRUTH MODEL STATES 


Shuttle Filter 

Truth Model 

Position 3 

Velocity 3 

Clock 

Phase 1 

Frequency 1 

1 

Position 3 
Velocity 3 
Thrust 3 
Clock 

Phase 1 
Frequency 1 
Aging 1 
Random Frequency 1 
Random Frequency Rate 1 

Total 8 

Total 14 

Position and velocity are in ECI (earth centered 
inertial) coordinates in both the Filter and the 
Truth Model. 


to model over 12 hour periods, especially when uploading is 
★ 

considered. 

The truth model contains three thrust states to 
account for venting and other unmodeled acceleration effects. 
The navigation filter has velocity white noise added to pre- 
vent it from being overly optimistic, since it has no thrust 
state error contributions. 


’''Satellite ephemeris and clock errors are minor for the single 
satellite study, but would be significant contributors for a 
multiple satellite study where greater navigation accuracy 
was possible. 
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3 . 3 PROCEDURE 

Three situations of varying Shuttle-to-GPS-satellite 
geometry are used in producing the navigation resuts reported 
here. These three situations are termed Cases 1, 2, and 3. 

As is mentioned above, in all cases the GPS satellite is in a 
circular, 63 deg inclination orbit at the appropriate radius 
to yield a 12 hr orbital period. The Shuttle orbit is also 
circular in all cases, at an altitude of 150 nm, so that it 
has a period of 90 min. In each case, a 12 hr span of GPS- 
aided navigation is simulated, with measurements at 2 min in- 
tervals whenever the satellite is visible. (The satellite is 
assumed visibile if it is at least 10 deg above the Shuttle 
horizon. ) The initial conditions in each case have the Shuttle 
and the GPS satellite at their respective ascending nodes cross 
ing the earth’s equatorial plane. In Cases 1 and 2, the Shuttl 
is in a 38 deg inclination orbit. In Case 1, the Shuttle right 
ascension is initialized to be equal to the initial right ascen 
sion of the GPS satellite, so that the Shuttle starts directly 
beneath the satellite. In Case 2, the Shuttle right ascension 
is moved by 60 deg, so it does not start directly beneath the 
GPS satellite. Case 3 is selected to make the orbit planes 
for the Shuttle and the GPS satellite coincide, so that cross- 
range information is not present in the GPS measurements. 

Thus, in Case 3, the Shuttle orbit has a 63 deg inclination, 
and. as in Case 1, the Shuttle starts directly beneath the 
satellite. Table 3.3-1 summarizes the geometry for these 
three cases. 

Data bases used for the navigation filter and for the 
associated truth model are summarized in Tables 3.3-2 and 3.3-3 
These values have all been used in prior GPS-aided Shuttle 
navigation studies, and the tables give references to their 
prior usage. 
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TABLE 3.3-1 

SHUTTLE -GPS SATELLITE GEOMETRY COMPARISON BY CASE 


ORBIT 

SPACECRAFT 

CASE 

PARAMETER 

1 

2 

3 

Inclination 

Shuttle 

38 

38 

63 

(deg) 

Satellite 

63 

63 

63 

Initial Right 
Ascension 
(deg) 

Shuttle 

1 

i o 

-60 

0 

Satellite 

o 

0 

0 


TABLE- 3.3-2 
FILTER DATA BASE 


ERROR SOURCE 

STANDARD 

DEVIATION 

DATA SOURCE 

Measurement Noise 
Pseudo Range 
Pseudo Range-Rate 

Clock Process Noise 
Phase 
Frequency 

Velocity White Noise 

12 ft 
0.05 fps 

0.1 ft 

3.0xl0’ 5 fps 
2.89xl0‘ 3 fps 

Ref. 2 

RMS ERRORS FOLLOWING INITIALIZATION (no measurements) 

Type 

Radial 

Downrange 

Crossrange 

Data 

Source 

Position (ft) 
Velocity (fps) 

2500 

26.3 

25000 

3.3 

2500 

3.3 

l 

Ref. 2 

i 

RMS Clock Errors 
Phase 1200 ft 

Frequency 0.25 fps 
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TABLE 3.3-3 
TRUTH MODEL DATA BASE 


ERROR SOURCE 


STANDARD DEVIATION 


DATA 

SOURCE 


Measurement Noise 
Pseudo Range 
Pseudo Range-Rate 

Clock Process Errors 
Initial Phase 
Initial Frequency 
Aging 

Random Frequency 
Random Frequency-Rate 


12 ft 
0.05 fps 

1130 ft 
0.142 fps 
3.86xl0‘ 6 ft/sec 2 
20.2 fps 
1.37 ft/sec 2 


Ref. 4 


Refs. 3 
and 4 and 
Chapter 2 


Thrust 

Axial Parameter 

Values 


Radial 

Down Range 

Cross Range 

Time Constant 
( sec ) 

600 

600 

300 

RMS Initial 
Thrust (lb) 

1.0 

1.0 

.33 

Process Noise 
( lb/vsec ) 

.0578 

.0578 

- 

.0272 

RMS Initial 
Conditions 

i i 

1 

1 

| 


i 

Position (ft) ; 

| 1040 

; ii3oo 

866 

Velocity (fps) 

12.5 

t 

1.09 

.094 


Chapter 2 


Ref . 4 
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3.4 RESULTS AND ANALYSIS 

Large initial errors are assumed at the start of ob- 
servation for all cases under consideration. Truth model rms 
position errors are assumed to be approximately 1000 ft radial, 
11.000 ft downrange, and 900 ft crossrange at initialization. 
Figure 3.4-1 displays error development for all three cases 
during the initial phase of the simulations. In all cases, 
the position errors have stabilized within 200 min (i.e., by 
the third visibility period). For Cases 1 and 3 where the 
first visibility period is short (18 min), rms downrange posi- 
tion errors grow to about 20,000 ft during the first visibility 
gap. The longer (36 min) initial visibility in Case 2 prevents 
rms downrange error from reaching 10,000 ft in the first gap. 

Steady state stabilized position errors for the three 
cases are presented in Fig. 3.4-2. RMS downrange position 
errors predominate and have maximum values, near 4,000 ft. at 
the end of the nonvisibility propagation phase. RMS radial 
errors are less than 900 ft and actually begin to decline after 
ibout 45 min (one-half an orbit) of nonvisibility. Radial and 
downrange errors fall below 500 ft rms in all cases at some 
point within the satellite visibility periods. Crossrange rms 
position error remains consistently below 500 ft for the first 
and second cases, with oscillations at a period of one h?lf 

tSt 

the orbital period. In Case 3, where the orbital geometry is 
such that the Shuttle orbit is always in plane with the single 
satellite orbit, no crossrange information is obtained. Cross- 
range rirs error therefore demonstrates a slow but steady growth 
due to thrusting, in addition to its oscillation. 


' f Crossrange error has a period equal to the orbital period. 
The graphs demonstrate a frequency doubling effect because 
they concern rms statistics on the error and not actual 
error values. 
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(b) Case 2 (c) Case 3 

Figure 3.4-1 Initial Phase RMS Position Errors 
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too SSO too 8S0 700 

Time imtm 


500 SSO 


600 650 

TIME tmm) 


(b) Case 2 


(c) Case 3 


Figure 3.4-2 


"Steady State" RMS Position Errors 
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Initial phase and steady state rms velocity errors 
for Case 1 appear in Fig. 3.4-3. The shape of these graphs 
when compared with Case 1 position errors in Figs. 3. 4- la and 

3.4- 2a illustrates important relationships between velocity 
and position error growth. Radial velocity and downrange posi- 
tion error curves display almost identical development. This 
correspondence is also observable w’hen comparing the downrange 
velocity with radial position error curves. It can be noted 
that the major rms radial velocity errors peak at about 4.5 
fps and all velocity errors drop to below 0.5 fps during part 
of the visibility segment. 

Clock rms phase and frequency errors for Case 1 are 
presented in Figs. 3.4-4 and 3.4-5, respectively. The phase 
errors increase slowly to the 1300 ft rms range when the sat- 
ellite is not visible and decrease to about 300 ft rms during 
visibility periods. Frequency error graphs demonstrate a quick 
rise to 0.34 fps rms when visibility is lost and lower values 
of 0.2 fps rms during visibility. Some of the abruptness of 
changes in frequency error as indicated by the truth model may 
be caused by the differences between the simplistic Shuttle 
filter two state clock and the truth model five state clock. 

An important consideration of the single satellite 
navigation results is the source of the errors. Error budgets 
for Case 1 at a time near maximum error, and at a time during 
good visibility, identify the important error sources. Table 

3.4- 1 presents these error budgets. As can be seen, initial 
conditions contribute little to velocity or position errors at 
this stage, about 1C hours after initiation of the GPS filter. 
Measurement errors are also a minor consideration during both 
visibility and nonvisibility time periods. At 590 min. when 
propagated errors in position and velocity are large, the 
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a) Initial Phase b) Steady State 


Figure 3.4-3 RMS Velocity F.rrors (Case 1) 

predominant influence is thrust. When the satellite is visible 
(610 min), thrust errors remain a major consideration, but clock 
errors become a dominant error contributor. An improvement in 
clock accuracy would cause a major improvement in position and 
velocity errors only during visibility periods when accuracy 
is already fairly good. 
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a) Initial Phase b) Steady State 


Figure 3.4-4 Clock Phase Errors (Case 1) 
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a) Initial Phase b) Steady State 


Figure 3.4-5 Clock Frequency Errors (Case 1) 
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3.5 CONCLUSIONS 

The key to accurate orbital navigation with a single 
GPS satellite is controlling error growth vhen the satellite 
is not visible . To some extent, this error growth can be con- 
trolled by limiting the duration of the visibility gaps them- 
selves. through constraints on the vehicle orbit. A second 
method for controlling error growth is to reduce the size of 
Shuttle vehicle thrust errors due to venting, to gravitational 
modeling errors, and to other unmodeled or incorrectly modeled 
acceleration sources. Any other system improvement -- in GPS 
receiver clock quality or in measurement accuracy -- would do 
very little to reduce peak navigation errors during visibility 
gaps. 


Vhen a single GPS satellite is visible, navigation 
errors are reduced to levels comparable to TDRSS navigation 
errors as reported in Ref. 23 . Visibility gaps with the two 
satellite TDRSS are relatively small, and no error growth 
characteristics comparable to those reported here are evident 
in the TDRSS results. However, it is not clear from the ref- 
erence that the assumed level of unmodeled vehicle acceleration 
is as high in the TDRSS simulation results as is assumed in 
this GPS study. 

in summary, a venting Shuttle vehicle may be navigated 
on orbit with only a single GPS satellite to useable accuracy 
provided that: (1) the orbit plane of the Shuttle does not 

coincide with that of the GPS satellite, and (2) exceptionally 
long visibility gaps, possible for some orbital geometries. 
are avoided. 
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U . POST BLACKOUT NAVIGATION WITH PHASE III GPS 


A major benefit of carrying GPS user equipment aboard 
the Space Shuttle is the timely provision of accurate naviga- 
tion data following reentry blackout. Chapter 2 of this re- 
port presents navigation accuracy results for Shuttle deorbit 
with the Phase I GPS constellation of six satellites that pro- 
vides coverage for deorbit for about a six hour period each 
day. This chapter considers navigation following reentry 
blackout with the Phase III GPS constellation of 24 satellites. 
Comparable accuracy would be anticipated if the Phase III GPS 
constellation were reduced to as few as 18 satellites. 

There are three additional changes from the scenario 
assumed for Chapter 2. First , a new Shuttle GPS antenna cover- 
age pattern (provided by Rockwell International for the Space 
Shuttle GPS panel) is used here. Second, the effect on navi- 
gation accuracy of incorporating a ground GPS transmitter to 
augment the satellite constellation is evaluated. (Chapter 5 
treats ground transmitters in considerably more detail.) Third , 
this study assesses the effect on navigation accuracy of three 
user clock (oscillator) types: the standard quartz clock, a 

quartz clock compensated for acceleration sensitivity, and a 
rubidium clock. 

Section 4.1 documents a variety of background ele- 
ments concerning deorbit navigation, principally those that 
are treated here in the same way as they were in past studies. 
The deorbit trajectory, the initial conditions on the naviga- 
tion system, the quality of the 1MU, the duration of reentry 
blackout, and the details of the navigation filter are all 
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discussed. Section 4.2 treats the issue of visibility, in- 
corporating the new antenna patterns, the Phase 111 GPS con- 
stellation, and the matter of measurement scheduling. Section 
4.3 deals with the three user clock types included in the study 
and Section 4.4 presents the navigation results and conclusions 


4 . 1 BACKGROUND 

Many of the elements used in this study of deorbit 
navigation were introduced and documented in Chapter 2 and in 
earlier work. Capsule descriptions of these elements are given 
here, and references to the earlier studies are provided. 

The overall methodology of evaluating aided naviga- 
tion errors is discussed in Chapter 4 of Ref. 3. In capsule, 
a covariance simulation of a suboptimal navigation filter is 
performed using a reference trajectory and an assumed meas- 
urement schedule. The filter gains are stored and used as 
input to a more comprehensive "truth model" covariance analy- 
sis program that produces the navigation accuracy results. 

The Shuttle trajectory used here is taken from a JSC- 
supplied trajectory tape and is an OFT-1 trajectory with a 
landing at Edwards AFB runway 17. It is discussed in Chapter 
2. The same conservative assumptions about the duration of 
reentry blackout are made here, blackout being assumed from 
300 kft to 145 kft of altitude, a period of nearly 15 min. 
During this segment, no GPS measurements are assumed possible. 

The navigation filter is an eight-state filter incor- 
porating three position error states, three velocity error 
states, one clock phase error state, and one clock frequency 
error state. This filter is the minimum necessary to process 
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GPS pseudo range and pseudo range-rate measurements. The cur- 
rent Shuttle GPS Panel baseline navigation filter adds either 
three platform misalignments or three acceleration error states. 

The IMU error model data base used in the truth model 
is the same as was used in the previous deorbit navigation 
studies, and corresponds to the IMU specifications for a pe- 
riod of time greater than 15 hrs since the last calibration. 

GPS measurement noises in both the navigation filter and in 
the truth model analysis program are 12 ft for pseudo range 
and 0.05 fps for pseudo range-rate. 

Since the Phase III GPS constellation offers general- 
ly excellent orbital visibility that results in excellent or- 
bital navigation, the initial navigation errors in position, 
velocity and the clock upon entry into reentry blackout are 
assumed to be negligibly small. The initial IMU platform mis- 
alignments, however, are not small. As in past studies, these 
misalignments, observable on orbit only by comparison with 
star tracker outputs, are assumed to be 284 seconds of arc 
(lo) about each of the three, space stable platform axes. 

The GPS satellite ephemeris and clock errors treated 
by the NESA program (Ref. 3) in some previous studies are not 
included in the truth model analysis of this study. Phase III 
GPS satellites are assumed to carry cesium clocks resulting in 
uploaded ephemeris errors that produce pseudo range errors of 
five feet or less (lo) to the user. These errors are negli- 
gible in the context of the current study. 


*A pseudo range-rate accuracy of 0.05 fps assumes a doppler 
count of approximately one second. 
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4.2 VISIBILITY 

There are three major issues affecting GPS satellite 
visibility from the shuttle during deorbit: 1) the receiving 

antenna gain pattern, 2) the GPS satellite constellation, and 
3) the time relationship between the deorbit trajectory and 
the satellite constellation. 

At the time this analysis was performed, the baseline 
the antenna configuration consisted of two antennas, one top- 
side with its maximum gain or symmetry axis pointing "straight 
up" in Shuttle body coordinates, the other a bottom antenna 
with its symmetry axis pointing straight down. The bottom an- 
tenna has a conical gain pattern; that is, antenna gain is a 
function only of the angle between the direction to the signal 
source and the direction of the antenna symmetry axis. For 
visibility purposes, such an antenna is treated as it if can 
receive all signals originating from directions within a cir- 
cular cone about the symmetry axis. The topside antenna is 
more complex. It has an elliptical visibility pattern in 
Shuttle coordinates. It is assumed to receive all signals 
originating inside of an elliptical cone about its symmetry 
axis . 


Table 4.2-1 contains the limits of visibility for 
these two antennas in terms of half-cone angles for two differ- 
ing sig-nal strength criteria. The primary zone, within which 
antenna gain will always exceed -1 dBiC (decibels ; ith respect 
to an isotropic antenna for right hand circularly polarized 
signals), is the zone in which usable signal reception is con- 
fidently assumed. In the larger secondary zone, usable signals 
might well be received, but the confidence is lower. The pri- 
mary zone is used in all of the results in this studv. 
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TABLE A. 2-1 

BASELINE ANTENNA VISIBILITY LIMITS 
(HALF-CONE ANGLES IN DEGREES) 


ANTENNA 

PRIMARY ZONE 
-1 dBiC POINT 

SECONDARY ZONE 
-4 dBiC POINT 


FORE-AFT 

LATERAL 

FORE -AFT 

LATERAL 

Topside 

00 

o 

o 

63° 

89° 

75° 

Bottom 

83° 

83° 

0 

r*- 

O0 

87° 


(See Figure 5.2-2, showing a slight variant of these 
specifications, for a better visualization). 

The Phase III GPS satellite constellation of 24 satel- 
lites has eight satellites uniformly spaced 45 deg apart in 
each of three orbital planes inclined at 63 deg to the equator 
and separated by 120 deg in right ascension. The satellites 
are in circular orbits at about four earth radii so they com- 
plete one orbit in just under 12 hrs and repeat their ground 
tracks once every sidereal day. Table 4.2-2 gives a set of 
right ascension and anomaly values for the 24 satellites at 
one instant. 

Missing from Table 4.2-2 is the "satellite epoch," 
the absolute time at which given orbital elements are valid. 
While there are baseline values for this epoch, they mean 
little to the present study. Rather, what is important is the 
time relationship between the satellite constellation and the 
deorbit trajectory. As part of the preliminary visibility 
studies, this time relationship was varied in a search for the 


*Recent changes put the inclination for Phase 111 at 55 deg 
and reduced the number of satellites to 18. Only minor dif- 
ferences in visibility, mostly improvement for this study, 
will result from the change in inclination. 
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TABLE *.2-2 

PHASE III GPS SATELLITE CONSTELLATION RIGHT ASCENSIONS 

AND INITIAL ANOMALIES 


ORBIT PLANE 1 
RIGHT ASCENSION 
0 deg 

ORBIT PLANE 2 
RIGHT ASCENSION 
120 deg 

ORBIT PLANE 3 
RIGHT ASCENSION 
240 deg 

SATELLITE 

NUMBER 

ANOMALY 

■ - - 

SATELLITE 

NUMBER 

ANOMALY 

SATELLITE 

NUMBER 

ANOMALY 

1 

0 

9 

15 

17 

30 

2 

45 

10 

60 

18 

75 

3 

90 

11 

105 

19 

120 

4 

135 

12 

150 

20 

165 

5 

180 

13 

195 

21 

210 

6 

1 225 

14 

240 

22 

255 

7 

270 

15 

285 

23 

300 

8 

315 

16 

330 

24 

345 


poorest visibility following reentry blackout. Visibility 
with the 24 satellites of the Phase III constellation is gen- 
erally excellent, and navigation accuracy given excellent vis- 
ibility is quite adequate for any Orbiter vehicle (as opposed 
to payload) navigation purpose. 

Figure 4.2-1 illustrates the worst case deorbit visi- 
bility encountered. As the key to the figures indicates, the 
thick line implies that a satellite is visible to the topside 
antenna (within the -1 dBiC gain pattern limit) whi'e the thin 
line implies visibility to the bottomside antenna. Following 
the indicated blackout segment of the deorbit trajectory, only 
three Satellites. Numbers 3, 4. and 15, are visible to the 
Space Shuttle. A little more than two minutes after the black 
out, several other satellites become visibile, due to an atti- 
tude change of the Orbiter. 

An important fact illustrated by Fig. 4.2-1 is that 
almost all of the satellites visible after blackout (with the 
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• PHASE III OPS CONSTELLATION. SITUATION SELECTED FOR POOR POST-BLACKOUT VISIBILITY 

• THICK SOLID LINE FOR TOPSIDE ANTENNA. THIN SOLID LINE FOR BOTTOM ANTENNA 
(BASELINE ANTENNAS) 



Figure 4.2-1 Worst Case Deorbit Visibility 

Following Blackout 


exception in this case of Satellite 23) are visible for some 
time segment prior to blackout . This means that if the GPS 
receiver/processor has storage space for the precision ephem- 
eris data on all 24 satellites (or a large enough subset of 
them), then precision GPS measurements may be taken after 
blackout without spending the 30 seconds on each satellite 
necessary to acquire the ephemeris data (see Chapter 6). This 
pre-acquisition of ephemeris data will require special consid- 
eration in the satellite selection algorithm, but will reduce 
the critical time-to-f irst fix following blackout by 60 sec for 
a two channel receiver, assuming four sat.ell.tes are required 
for the fix. 
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From the visibility results shown in Fig. 4.2-1. a 
measurement schedule was prepared that takes into account the 
time to acquire and track the signal from each satellite, omit- 
ting the 30 sec ephemeris data acquisition time. Figure 4.2-2 
shows this measurement schedule in terms of satellite availa- 
bility. Measurements were scheduled at 6 sec intervals to the 
available satellites. 




CO 

i 

z 


P « 

_j 

ui 30 

t- 
« 

W 21 -j 
22 -1 


GT H 


* 1 1 1 ' — l 1 1 1 

36 39 40 41 42 43 44 

MISSION TIME (mint 

Figure 4.2-2 Measurement Schedule After Blackout 

(Worst Case Visibility) 


To evaluate performance sensitivity to Shuttle antenna 
coverage, satellite avai 1 abi 1 i t les for a reduced antenna cover- 
age situation are also used to generate performance results, 
figure 4.2-3 shoves these reduced coverage availabilities, which 
correspond to antennas with conical viewing limits of 55 cleg 
half angles. With reduced coverage, only two satellites are 
available immediately following blackout. 
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Figure 4.2-3 Reduced Coverage Satellite Availability 

After Blackout (Including Ground 
Transmitter) 

Because of the low coverage with reduced visibility, 
the possibility of using a ground GPS transmitter is considered 
in some cases. Figure 4.2-3 also shows when the ground trans- 
mitter (assumed located under the Orbiter trajectory at San 
Simeon) would be available. No antenna limitations on ground 
transmitter visibility are assumed, due to the potentially 
high power of the ground transmitter; only line of sight lim- 
itations from the transmitter to the Orbiter are taken. 

4.3 CLOCK QUALITY 

Three different clock types are included in this study 
Figure 4.3-1 illustrates the quality of two of them -- the ordi 
nary (uncompensated) quartz clock and the rubidium clock -- in 
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Figure 4.3-1 Quartz and Rubidium Oscillator Quality 

an unaccelerated environment. Clock modeling is discussed in 
Chapter 4 of Ref. 3. The same five-state clock model (with 
different data bases) is used to simulate both the quartz and 
the rubidium oscillators. 

An additional source of clock error arises from accel- 
eration sensitivity. Ref. 6 discusses such sensitivity. A 
quartz crystal has an axis of sensitivity to non-gravi tat ional 

acceleration, and the frequency error produced bv the oscilla- 

q 

tor is typically 1 part in 10 per g of acceleration along the 
sensitive axis. Accelerations normal to the sensitive axis 
produce little or no effect. 

Rubidium clocks depend for long term (10 sec and up > 
stability on the atomic rubidium standard, but they use a quartz 
crystal oscillator for short term (1 sec and less' stability. 
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The rubidium cell is put into an outer loop about the quartz 
crystal, and the time constant of this loop (from 1 to 10 sec) 
determines the intermediate term stabiliity. The rubidium cell 
is not subject to appreciable acceleration sensitivity, being 
at least two orders of magnitude less sensitive than is the 
quartz crystal. Tracking stability of the GPS receiver depends 
upon short term stability of the oscillator, so quartz and 
rubidium clocks are equal in this respect. But the ability of 
a sequential receiver to "tie together" measurements at differ- 
ent times depends on the long term stability of the oscillator, 
so rubidium oscillators provide an improvement over quartz 
oscillators in this respect. In particular, a rubidium oscil- 
lator can maintain a good time reference through the blackout 
period of 15 min, so that only three satellites are necessary 
to obtain a good "fix" after blackout. 

Much of the error that a quartz crystal clock would 
acquire during the blackout period would arise from the accel- 
eration sensitivity. As Ref. 24 points out, this sensitivity 
may be calibrated and accounted for by mounting an accelerom- 
eter with the crystal and building a compensation control sys- 
tem about the oscillator. Figure 4.3-2 diagrams such a control 
system. Alternately, for the Space Shuttle, acceleration com- 
pensation could be accomplished in the GPS receiver/processor 
software using acceleration and attitude inputs from the IMU’s. 

The quality of such compensation is potentially quite 
good. Figure 4.3-3 shows experimental results reported in Ref. 
24, while Figure 4.3-4 indicates the range of frequencies over 
which the compensation is valid (the primary limitation being 
the accelerometer and not the math model). For navigation pur- 
poses, as previously remarked, it is the long term (low fre- 
quency) stability that is important. Because of this potential 
for improving the quartz clock, the third type of clock con- 
sidered here is an acceleration compensated quartz oscillator. 
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ACCELERATION 



A33>M 




OSCILLATOR 

PACKAGE 


I I 


Figure 4.3-2 Acceleration Compensated Quartz Oscillator 





Figure 4.3-3 Acceleration Compensation Test Result (Ref. 24) 
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Figure A. 3-4 Acceleration Compensation Frequency 

Dependence (Ref. 2A) 

q 

A frequency error of one part in 10 is equivalent to 
a pseudo range-rate error of about 1 ft/sec (since the speed 

Q 

of light is about 10 ft/sec). Since acceleration levels on 
the Orbiter during deorbit are limited to three g's or less, 
and since compensation is good to 2%, the remaining acceler- 
ation sensitivity of a compensated quartz or rubidium oscilla- 
tor for navigation purposes may be ignored. 


A. 4 POST BLACKOUT NAVIGATION RESULTS 

The combination filter covariance/truth model covari- 
ance analysis programs were run on several different post-black- 
out scenarios varying the post-blackout GPS transmitter availa- 
bility and the user clock type. Fortunately, the results may 
be summarized without having to look in detail at the case con- 
ditions and output covariances versus time for each scenario. 

A series of tables will suffice. 

Table A. 4-1 shows the clock errors after blackout. 

These errors change very little in the time period following 


4-13 


THE ANALYTIC SCIENCES CORPORATION 


TABLE 4.4-1 

POST- BLACKOUT CLOCK ERRORS* 


RMS CLOCK 
ERROR 

QUARTZ 

COMPENSATED 

QUARTZ 

Rb 

PHASE (FT) 

900 

125 

8 


(40) 

(30) 


FREQUENCY (FT/SEC) 

1.3 

0.25 

0.005 


(0.25) 

(0.20) 



blackout until the fourth emitter (satellite or ground trans- 
mitter) is acquired, whereupon the filter obtains the informa- 
tion necessary to correct the relatively large errors of the 
quartz oscillators, and the values in parentheses are obtained. 
This table shows that the dominant portion of the quartz clock 
errors accumulated during the blackout are due to acceleration 
sensitivity of the clock , for the only difference between the 
quartz clock and the compensated quartz clock is the accelera- 
tion sensitivity. 

Table 4.4-2 is an error budget for Shuttle position 
and velocity errors following the blackout period before any 
GPS measurements are taken. The budget assumes that an ordi- 
nary quartz oscillator is used in taking GPS measurements prior 
to the blackout period. If a compensated quartz or rubidium 
oscillator were used prior to blackout, the budget line items 
for the clock and for the clock acceleration term would change 
to smaller values, but the totals would be unaffected, since 
clock errors do not contribute significantly to post-blackout 
navigation error until after GPS measurements are processed by 
the filter. 
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TABLE 4.4-2 

POST- BLACKOUT ERROR BUDGET PRIOR TO MEASUREMENTS 



RMS ERRORS i 

ERROR 

GROUP 

POSITION (FT) 


VELOCITY (FPS) 



DR 

CR 

R 

DR 

CR 

R 

k 

Measurement Noise 

366 

252 

670 

0.33 

0.12 

1.28 

+ 

Platform Alignment 

3018 

9073 

12030 

14.60 

24.58 

35.28 

IMU 

1305 

1933 

2600 

4.41 

5.74 

8.11 

Clock" 

103 

104 

319 

0.11 

0.07 

0.59 

k 

Clock Acceleration 

10 

10 

30 

0.01 

0.01 

0.06 

RSS 

3310 

9280 

12330 

15.26 

25.24 

36.19 

Vector RSS 


15790 



46.69 



*These errors have propagated from pre-blackout conditions. 

tPlatform alignment error is initial value. Subsequent growth during 
blackout is included in IMU error group. 

In fact, the only significant contributors to post- 
blackout navigation errors (before GPS measurements are taken) 
are the initial IMU platform alignments going into blackout and 
the IMU component (gyro and accelerometer) errors accumulated 
during blackout . Gyro errors during blackout increase the IMU 
platform misalignment errors, but the budgeted navigation errors 
due to platform misalignments refer to those errors due only 
to the initial values of the misalignments, with the errors 
due to subsequent growth being lumped into the IMU entry. 

Tables 4.4-3 and 4.4-4 show how navigation errors de- 
crease after the acquisition of independent emitters following 
blackout. In this general scenario, the navigation error is 
not a strong function of time, per se, or of the number of 
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TABLE 4.4-3 

POST-BLACKOUT NAVIGATION ERRORS AFTER 
ACQUISITION OF 1ST AND 2ND EMITTERS 


RMS VECTOR 
ERROR 

BEFORE 

MEASUREMENT 

AFTER 

MEASUREMENTS TO 
1 EMITTER 

AFTER 

MEASUREMENTS TO 
2 EMITTERS 

Position (ft) 
Velocity (ft/sec) 

16,000 

47 

13,000 

40 

5,000 

14 


TABLE 4.4-4 

POST- BLACKOUT NAVIGATION ERRORS AFTER 
ACQUISITION OF 3RD AND 4TH EMITTERS 


CLOCK 

TYPE 

RMS VECTOR 
ERROR 

AFTER 

MEASUREMENTS TO 

AFTER 

MEASUREMENTS TO 


Position (ft) 

1500 

60 

QUARTZ 

Velocity (ft/sec) 

2.0 

1.5 

COMPENSATED 

Position (ft) 

180 

40 

QUARTZ 

Velocity (ft/sec) 

1.3 

1.3 


Position (ft) 

30 

25 

Rb 

Velocity (ft/sec) 

1.0 

1.0 


measurements (several measurements to one emitter being pos- 
sible), but is primarily a function of the number of different 
emitters to which measurements have been made . The results 
from which Table 4.4-3 was compiled show, by implication, that 
there is no substantial variation in navigation error with user 
clock type when only one or two emitters have been acquired . 

In Table 4.4-4, on the other hand, the clock type is explicit, 
and variations in navigation error are noted. The clock type 
is especially important when three emitters have been acquired . 
This is to be expected. Without accurate time information 
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(the '"normal" case) it takes four GPS satellites to determine 
a user fix. With accurate time information, three satellites 
suffice . 


The only time (in this limited study) when an accurate 
user clock is important is when only three emitters are avail- 
able following blackout. In the same context, the only time a 
ground GPS transmitter is useful is when three satellites are 
available and an uncompensated quartz oscillator is used, or 
when two satellites are available and a compensated quartz or 
rubidium oscillator is used . 

Table A. 4-5 is an error budget following acquisition 
of the third GPS emitter. It shows the dominance of the clock 
errors when an uncompensated quartz clock is used. Table 4.4-6 
is an error budget following acquisition of the fourth emitter. 
After the fourth emitter has been acquired, position errors 
are reduced to small values, but velocity errors are not cor- 
respondingly small. The dominant source of the velocity er- 
rors is platform misalignment. (The 1MU entry is also mostly 
platform misalignment growth during blackout.) The eight- 
state baseline navigation filter does not recognize platform 
alignment errors and thus can do little to mitigate their 
effects. A significant reduction in post-blackout velocity 
errors would require the addition of three platform misalign- 
ment states to the navigation filter / Additionally, a com- 
pensated quartz or rubidium clock would help. Measurement 
noise contributions to velocity errors will fall as the number 
of measurements increases, while the other error contributions 
will not. 


*These misalignment states have subsequently been incorporated 
into the baseline filter design. 
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TABLE 4.4-5 

ERROR BUDGET FOLLOWING ACQUISITION OF THE 3RD EMITTER 
(UNCOMPENSATED QUARTZ OSCILLATOR) 


ERROR GROUP 



RMS 

ERRORS 



POSITION (FT) 

VELOCITY (FPS) 

DR 

CR 

R 

DR 

CR 

R 

Measurement Noise 

25 

29 

39 

0.41 

0.55 

0.25 

Platform Alignment 

3 

2 

1 

0.40 

0.32 

0.06 

IMU 

1 

0 

0 

0.13 

0.10 

0.03 

Clock 

83 

100 

139 

0.15 

0.15 

0.23 

Clock Acceleration 

633 

769 

1065 

0.80 

0.76 

1.41 

RSS 

639 

766 

1074 

1.00 

1.01 

1.45 

Vector Rss 


1470 



2.04 



TABLE 4.4-6 

ERROR BUDGET FOLLOWING ACQUISITION OF THE 4TH EMITTER 
(UNCOMPENSATED QUARTZ OSCILLATOR) 


ERROR GROUP 



RMS 

ERRORS 



POSITION 

(FT) 

VELOCITY (FPS) 

DR 

CR 

R 

DR 

CR 

R 

Measurement Noise 


11 

29 

0.60 

0.11 

0.34 

Platform Alignment 

8 

3 

15 

0.84 

0.39 

0.25 

IMU 

2 

1 

0 

0.24 

0.11 

0.08 

Clock 

3 

1 

6 

0.16 

0.12 

0.24 

Clock Acceleration 

9 

3 

44 

0.49 

0.37 

0.54 

RSS 

22 

12 

53 

1.19 

0.58 

0.73 

Vector RSS 


58 



1.51 



4-18 























THE ANALYTIC SCIENCES CORPORATION 


5. GPS GROUND TRANSMITTERS FOR SPACE SHUTTLE SUPPORT 


During ascent and entry of the Space Shuttle, GPS sat- 
ellite visibility may be inadequate to support accurate naviga- 
tion. Depending upon the type of clc~k used in the GPS receiver, 
either three satellites (for an atomic clock) or four satellites 
(for a quartz crystal clock) must be in view of the Orbiter 
antennas to accurately navigate during these dynamic flight 
phases. Before the Phase III GPS constellation reaches initial 
operational capability (18 satellites) in 1986* there will at 
most times be fewer than three satellites visible to the Orbiter. 
Even after the Phase III constellation is complete, the external 
fuel tank will block the downward GPS antennas on the Orbiter 
during ascent. Since the Orbiter flies "upside down" during 
ascent, the fuel tank may severely impact satellite visibility. 

This chapter examines some of the issues involved in 
using GPS ground transmitters to augment the satellite constell- 
ation and support Shuttle navigation during ascent and entry. 
Section 5.1 briefly examines some of the issues involved, while 
Sections 5.2 and 5.3 present results for the entry and ascent 
phases, respectively. Section 5.4 contains a summary of con- 
clusions regarding the use of ground transmitters. 


5.1 GROUND TRANSMITTER ISSUES 

A wide range of design, performance, and support issues 
would need resolution to establish the feasibility of GPS ground 


^^Restructuring of the GPS program may result in a final Phase 
111 constellation of 18 satellites, to be operational in 1987. 
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transmitters as an aid to the Space Shuttle during ascent and 
entry. This report addresses a subset of those issues regard- 
ing the use of the satellite signals, a subset that comes gener- 
ally under the headings of signal availability, signal usability, 
and navigation accuracy. In particular, the following questions 
are examined: 

• How many ground transmitters would be 
required, and where would they best be 
placed? 

• What demands would ground transmitters 
place on the antenna management and track- 
ing loop functions of the Orbiter GPS 
receiver processor assembly (R/PA)? 

• What Orbiter navigation accuracy could 

be expected to result from the use of 
ground transmitters? 

The first question, that of number and placement of 
ground transmitters, depends, in turn, on a range of further 
issues, including: 

• fixed versus mobile transmitters 

• the range (variety) of ascent or entry 

trajectories to be accommodated 

• the type and quality of the R/PA clock 

• the GPS antenna coverage of the Orbiter 

• the interim GPS satellite constellation(s ) 

• the attitude dynamics of the Orbiter 

The number of transmitters required might also be expected to 

be related to the navigation accuracy demanded of the system, 

but the actual relationship here is a simple one: depending 

on the type of clock in the R/PA, if three or four signals are 
available, navigation accuracy will be good: if fewer than three 
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or four signals are available, then (GPS-based) navigation 
accuracy will be less than adequate unless the navigation system 
Kalman filter can take advantage of changing transmitter- to- 
Orbiter geometry. So the question of number and placement of 
ground transmitters generally is resolved on matters of signal 
availability (visibility) alone. 

The second question, that of the demands a ground 
transmitter signal places on the receiver, is partly a matter 
of visibility -- Is the signal available to one antenna long 
enough to be acquired and tracked? -- is partly a matter of 
signal dynamics -- Can the tracking loops follow the signal 
dynamics, especially the geometrical range acceleration arising 
f-om high Orbiter velocities relative to nearby ground transmit- 
ters? -- and is partly a matter of signal strength relationships 
Can the ground transmitter signal strength be selected so as 
to be usable over the large relative variation in range to the 
Orbiter and yet not interfere with other ground transmitters 
or with the satellite signals? These questions may be answered 
with visibility charts, range dynamics computations , and range 
ratio computations. 

The third question, that of expected accuracy, is not 
addressable with just the visibility analysis tools useful for 
the other issues considered here. The source of the GPS signals 
has far less impact on Shuttle accuracy than does the number 
of such signals, so previous results for post-blackout entry 
accuracy will give a good idea of the accuracy to be expected 
in this flight phase if GPS ground transmitters were substituted 
for satellites. Only limited performance analysis has been 
done for the ascent phase, so no comparable accuracy estimates 
are presented here for that phase. 
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5 . 2 ENTRY 

In investigating the use of ground transmitters to 
support entry, it would be desirable to have a full range of 
entry trajectories available so that variations in signal 
availability could be considered. Unfortunately, only one 
entry trajectory tape, for a landing at Edwards Air Force Base 
(EAFB) (the tape used in all previous deorbit studies and dis- 
cussed in Section 2.2.2), was available for this study. The 
problem of signal availability is determined by the relative 
locations of the ground transmitters to the trajectory, so by 
assuming numerous ground transmitter sites, one trajectory can 
be made to show a wide variety of situations. Figure 5.2-1 is 



Figure 5.2-1 Deorbit Ground Transmitter Locations 
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a map showing the ground track of the trajectory utilized for 
this study and the locations of the 14 ground transmitter 
sites simulated. The numbers labeling the tic marks along the 
trajectory ground track are minutes to touchdown. 

Due to the mountainous terrain in the wide vicinity 
of the landing site, high elevation locations were often 
chosen for the ground transmitters to avoid the question of 
terrain blockage of visibility. Some of the exact locations 
assumed might not be practical, but equivalent, practical lo- 
cations are almost surely available. In all cases in this 
entry study it is assumed that the Orbiter is visible to the 
ground transmitter if the Orbiter is over five degrees above 
the local horizon at the transmitter. This does not, however, 
necessarily imply signal availability. In addition to the 
Orbiter being visible to the ground transmitter, the ground 
transmitter must be within the beamwidth of the Orbiter GPS 
antennas in order for the signal to be received. 

Figure 5.2-2 illustrates the antenna coverage assumed 
for the Orbiter for the entry study (Ref. 25). This is approx- 
imately the same coverage assumed for the deorbit analysis in 
Chapter 4. (The ascent study, done later, uses different an- 
tenna patterns resulting from model tests at JSC. Test re- 
sults were not available when the entry study was done.) For 
ground transmitter visibility during entry, it is the bottom 
antenna that is of primary importance. Note that the bottom 
antenna pattern assumed is an axisymmetrical "conical” pattern, 
while the top antenna pattern is an asymmetrical "elliptical" 
pattern . 


The maximum range of visibility from a ground trans- 
mitter to the Orbiter is about 240 nm. and Fig 5.2-1 illus- 
trates the point of first contact for three of the ground 
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I ! R40I2K 


TOP ANTENNA AXIS TOP ANTENNA AXIS 



BOTTOM ANTENNA AXIS BOTTOM ANTENNA AXIS 

i : 


TOP ANTENNA "ELLIPTICAL" 

BOTTOM ANTENNA "CONICAL" 

Figure 5.2-2 Antenna Coverage Assumed for Entry 

(-1 dBi Gain Limits) 

transmitters. This point occurs following the (assumed) exit 
of down-looking L-band communications blackout. (GPS signals 
as well as TACAN signals are in the L-band.) Only the final 
12 min or so of the trajectory prior to landing need be con- 
sidered for ground transmitter aiding. 

Figure 5.2-3 summarizes the visibility situation for 
all 14 ground transmitter sites. As the key describes, the 
narrow lines indicate periods when the Orbiter is below the 
five deg elevation minimum to be visible from the ground trans- 
mitter. During the wide portions of the lines, the Orbiter is 
visible to the ground transmitter, but the ground transmitter 
is not always visible to the Orbiter GPS antennas. During the 
unshaded portions of the wide lines, the ground transmitter is 
not visible to either Orbiter GPS antenna. During the solid 
portions of the wide lines, the ground transmitter is visible 
to the bottom antenna, and during the cross hatched portions 
of the line, it is visible to the top antenna. 
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MINUTES TO TOUCHDOWN 




KEY 


BELOW S <Nf 

ELEVATION 

f— | NOT VISIBLE TO 
* — EITHER ANTENNA 
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BOTTOM ANTENNA 


EE 3 


VISIBLE TO 
TOP ANTENNA 


-1 dB* ANTENNA PATTERN 
LIMITS ASSUMED FOR 
VISIBILITY 


Figure 5.2-3 Entry Visibility of Ground Transmitters 


In preparation for landing, the Orbiter undergoes two 
sharp attitude changes, noted in Fig. 5.2-3, that produce numer 
ous simultaneous changes in the ground transmitter visibilities 
In general, these attitude changes arc not totally predictable 
as to timing or amount, but the figure illustrates the type of 
unpredictable visibility situation that must be dealt with 
during a typical entry. 

Alt of the 14 ground transmitters would not be needed 
to support entry, as only three or four signals are necessary 
for accurate navigation. If no signals were available from 
the GPS satellite constellation, then three to five ground 
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transmitters would be necessary for the specific trajectory 
studied, depending on the receiver clock and on the period of 
coverage desired. Generally speaking, a period of about two 
minutes of simultaneous availability of three signals as early 
in the trajectory as possible would provide most of the poten- 
tial benefit to be derived from ground transmitters (assuming 
an atomic clock in the R/PA). The sites at Cone Peak, Edwards, 
and Magic Mountain would provide such signal availability. 

If a range of entry trajectories were to be consid- 
ered, and fixed ground transmitters were to be used, then, of 
course, more than three would be required. If the entry tra- 
jectory were known sufficiently in advance, the use of mobile 
ground transmitters could be considered. There are two ways 
to use mobile transmitters. One involves locating the mobile 
transmitter and synchronizing its clock by reference to the 
GPS satellite system. During GPS Phase II, this might take 
the better part of a day to accomplish, as the six (or fewer) 
satellites are not visible to ground locations except for a 
period of a few hours a day. A second way to set up mobile 
transmitters would be to locate them at previously surveyed 
sites and depend on atomic time standard accuracy to maintain 
clock synchronization to each other (and to the satellite sys- 
tem, if necessary). In this mode, mobile transmitters might 
be set up in a period of a few hours. In either case, mobile 
transmitters might not be able to support contingency entries. 

A second vantage point on ground transmitter visibil- 
ity is provided in Table 5.2-1. This table ignores the effect 
of the Orbiter antennas in limiting visibility of the ground 
transmitters and assumes the transmitters are visible if the 
Orbiter is above five degrees elevation when viewed from the 
transmitter. The duration of the visibility period is noted 
for each of the 14 ground transmitter sites, as is the peak 
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TABLE 5.2-1 

GROUND TRANSMITTER VISIBILITY RANGES 
(ANTENNA EFFECTS NOT INCLUDED) 




| GROUND TRANSMITTER 

VISIBILITY 

DURATION 

(MIN) 

PEAK 

ELEVATION 

(DEC) 

AZIMUTH 

VARIATION 

(DEG) 

NO. 

location 

1 

LICK OBSERVATORY 

3.2 

M.O 

64 

2 

VISALIA 

6.0 

11.5 

103 

3 

MATUKANGO PEAK 

4.4 

7.0 

43 

4 

CONE PEAK 

5.2 

20.8 

124 

5 

ALDER PEAK 

5.5 

25.2 

127 

6 

CJUNA LAKE 

5.0 

11.9 

49 

7 

BAKEKSPIELD 

7.4 

28.5 

156 

8 

EDUARDS 

8.4 

27. A 

96 

9 

BARSTOW 

5.8 

9.0 

13 

10 

vandenberc 

6.5 

30.1 

170 

11 

TWENTYN1NE PALMS 

1.3 

5.2 

6 

12 

MAGIC MOUNTAIN 

6.9 

19. b 

77 

13 

RIVERSIDE 

4.6 

7.8 

38 

14 

PENDLETON 

2.9 

6.0 

27 


elevation of the Orbiter, and its variation in azimuth as viewed 
from the transmitter. The elevation and azimuth variation 
values are related to navigation accuracy. In particular, 
accuracy in the vertical direction depends on the elevation 
angles, higher angles producing better accuracy. 

Variations in azimuth during the viewing period can 
improve overall accuracy. Since the ground transmitter ob- 
servations are assumed to be processed by a Kalman filter that 
has access to IMU-sensed velocity information, the cumulative 
effect of two measurements in different directions is much 
stronger than two measurements in a single direction. That 
is, one ground transmitter viewed from two different direc- 
tions during the trajectory has about the same effect as two 
independent ground transmitters in different directions. This 
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aspect of navigation accuracy with ground transmitters was not 
studied explicitly, but it might be concluded that even fewer 
than three ground transmitters could support navigation during 
entry. Fewer than three transmitters would not support ac- 
curate navigation at the earliest opportunity, however, for to 
be effective with less than three transmitters requires wait- 
ing until substantial changes in the Orbiter- to- transmitter 
geometry have occurred. 

Accuracy results from the Chapter 4 study on this 
same trajectory are appropriate for the first establishment of 
accurate navigation following entry blackout. Table 5.2-2 (a 
repeat of Table 4.4-3 and 4.4-4) summarizes these results. 
Following the blackout, typical navigation errors would be 
16,000 ft. and 47 ft/sec. Regardless of the type of clock used 
in the R/PA, these errors will be reduced when the first two 
(independent) GPS signals are processed, as the first part of 
the table shows. The results after processing the third and 
fourth independent GPS signals depend on the type of clock in 
the R/PA, as the second part of the table shows. The accel- 
eration-compensated quartz clock refers to a quartz crystal 
oscillator which is compensated by special added hardware (or 

TABLE 5.2-2 

POST-BLACKOUT NAVIGATION ERRORS 
FOLLOWING ACQUISITION OF 
FIRST AND SECOND EMITTERS 


RMS VECTOR 
ERROR 

BEFORE 

MEASUREMENT 

AFTER 

MEASUREMENTS TO 
1 EMITTER 

AFTER 

MEASUREMENTS TO 
2 EMITTERS 

Position (ft) 

16,000 

13,000 

5.000 

Velocity (ft/sec) 

47 

40 

14 
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TABLE 5.2-2 (Cont.) 
(THIRD AND FOURTH EMITTERS) 


CLOCK 

TYPE 

RMS VECTOR 
ERROR 

AFTER 

MEASUREMENTS TO 
3 EMITTERS 

AFTER 

MEASUREMENTS TO 
4 EMITTERS 


POSITION (FT) 

1500 

60 

QUARTZ 

VELOCITY (FT/SEC) 

2.0 

1.5 


POSITION (FT) 

180 

40 

QUARTZ+ 

VELOCITY (FT/SEC) 

1.3 

1.3 


POSITION (FT) 

30 

25 

Rb 

VELOCITY (FT/SEC) 

1.0 

1.0 


through R/PA software) so as to reduce its acceleration sensi- 
tivity by about an order of magnitude. The ’'Rb" heading refers 
to a rubidium time standard. Such an atomic clock has a negli- 
gible acceleration sensitivity for navigation purposes. This 
table illustrates that accurate navigation will be achieved 
when the third or fourth independent signal is processed, de- 
pending on the type of clock in the R/PA. 

One more result of the entry ground transmitter study 
is given in Table 5.2-3, where the ratios of the minimum and 
maximum ranges from the Orbiter to the various ground transmit- 
ters during the periods of potential visibility are expressed 
in terms of an equivalent power ratio (20 log^Q rmax/rmin). 

The various C/A codes used for GPS signals can interfere with 
each other if they do not have the proper signal power rela- 
tionships. The peak cross correlation of one C/A code to 
another is -21.6 dB (Ref. 26), which means that if one signal 
is close to 21.6 dB stronger than another, it will interfere 
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TABLE 5.2-3 

GROUND TRANSMITTER SLANT RANGE VARIATIONS 




GROUND TRANSMITTER 

SLANT RANGE 
ABOVE 5 DEC ELEVATION 

NO. 

LOCATION 

MAX 

(MI) 

MIN 

(MI) 

POWER 

RATIO 

(dB) 

1 

LICK OBSERVATORY 

230 

180 

2.1 

2 

VISALIA 

230 

97 

7.5 

3 

MATURANGO PEAK 

186 

82 

7.1 

4 

CONE PEAK 

242 

65 

11.4 

5 

ALDER PEAK 

242 

55 

12.9 

6 

CHINA LAKE 

210 

52 

12.1 

7 

BAKERSFIELD 

232 

35 

16.4 

8 

EDWARDS 

212 

2 


9 

BARSTOV 

196 

51 

11.7 

10 

vandenberg 

. 240 

42 

15.1 

11 

TWENTYN1NE PALMS 

152 

124 

1.8 

12 

MAGIC MOUNTAIN 

220 

48 

13.2 

13 

RIVERSIDE 

192 

84 

7.2 

14 1 

PENDLETON 

185 

124 

3.5 


PEAK C/A CODE CROSS CORRELATION (WITH DOPPLER SHIFT) = -21.6 dB 


with the acquisition process for the other code. As long as 
the power ratio of the minimum and maximum ranges for a single 
transmitter is safely below this figure, it would be possible 
to select a noninterfering transmitter power level that is 
constant. When the power ratio climbs to within a few dB of 
21.6 dB, however, such a power level setting is not possible, 
and either: 


• the transmitter power level must be dy- 
namically adjusted in response to range 
to the Orbiter during entry. 

• the transmitter power level may be kept 
constant but will be inadequate for 
acquisition at maximum range. 
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• the transmitter power level may be kept 
constant but will interfere with other 
signals at minimum range. 

As the table shows, only the transmitter located at the land- 
ing site poses any interference problems. 

* 

One final issue of signal usability was not investi- 
gated as part of the entry study. It concerns the dynamics of 
the signals from the ground transmitters and the stress they 
place on the receiver tracking loops. This stress is consid- 
erable and would substantially affect the receiver design, as 
is shown in the next Section, where it is considered for the 
ascent phase. 


5 . 3 ASCENT 

Figure 5.3-1 (Ref. 27) illustrates the range of allow- 
able launch azimuths from the Kennedy Space Center. Three 
ascent trajectories (Ref. 28) were available for this study 
that reasonably span the range of those allowable. Two are 
termed northward launches, one to a nearly limiting 55 deg 
orbital inclination, and another to a more normal 40.5 deg 
orbital inclination (sometimes referred to as 40 deg incli- 
nation). A third trajectory is eastward, to a 28 deg orbital 
inclination. This eastward ascent exhibits about the poorest 
ground transmitter visibility of all launches due to the scar- 
city of locations for ground transmitters in this direction. 

A total of 17 ground transmitter sites were selected 
for study although only 14 were used at any one time. Figure 
5.3-2 maps the three ascent trajectory ground tracks and shows 
the ground transmitter locations, although it omits the names 
of the ground transmitters. The length of the trajectory ground 
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Figure 5.3-1 Allowable Orbit Inclinations and Launch 

Azimuths from Kennedy Space Center 








Figure 5.3-2 Map of Ascent Trajectory Ground Tracks And 

GPS Ground Transmitter Study Locations 
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tracks is indicative of how much of the trajectory was avail- 
able for study. In the 40.5 deg inclination case, only that 
portion of the trajectory prior to Orbiter main engine cutoff 
(MECO) was available. Post-MECO data was available for the 
other two trajectories. 

To correlate the locations of the ground transmitters 
with the names of the ground transmitter sites given in the 
sequel, it should be noted that the transmitters are always 
listed from north-most to south-most, and Bermuda is an 
obvious reference point included in all lists. Bermuda is just 
slightly north of Savanna Beach, Georgia, which appears at an 
equal latitude in the figure. Cape Kennedy and Patrick AFB 
are so close to one another as to be indistinguishable on the 
map in the figure, where two arrows point to the same location. 

There are no terrain visibility restrictions for any 
of these ground transmitter sites, and so the visibility cri- 
terion used in the entry study was slightly relaxed here -- 
during ascent the Orbiter must be above three degrees above 
the horizon at the ground transmitter to be visible to the 
ground transmitter. Of course, the ground transmitter must 
also be within the receiving antenna beamwidth of the Orbiter 
for signal reception. Figure 5.3-3 illustrates the visibility 
situation for ground transmitters during ascent. Note that the 
antenna pattern has also been altered from that used in the 
entry study. Now an 80 deg half-cone angle, conical antenna 
is assumed as nominal (Ref. 29), while a full hemispherical 
coverage, 90 deg half-cone angle antenna is used as an alter- 
nate to study the sensitivity of the visibility results to 
antenna coverage. 

Ground transmitters were studied for entry primarily 
to assist during GPS Phase II when no satellites are available 
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Figure 5.3-3 GPS Ground Transmitter Visibility Geometry 

During Ascent 

for significant portions of each day. For the ascent study, 
the external fuel tank will block satellite visibility to the 
downward antenna even during GPS Phase III when satellites are 
constantly available. But satellites are sometimes visible to 
the upward antenna during ascent, so the contribution of the 
satellites to the ascent signal availability must be examined. 
The visibility criteria for satellite signals are illustrated 
in Fig. 5.3-4. Because satellite signals have limited strength 
t^e horizon mask angles are more conservative for satellites 
than for ground transmitters. 

Figures 5.3-5 through 5.3-7 summarize the ground trans- 
mitter visibility situations for the three ascent trajectories. 
The key to these figures is almost the same as that for Fig. 
5.2-3. except that three deg elevation is now the visibility 
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Figure 5.3-4 GPS Satellite Visibility Geometry During Ascent 
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Figure 5.3-6 Ground Transmitter Visibility For KSC Ascent 

(40 Deg Inclination Orbit) 
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Figure 5.3-7 Ground Transmitter Visibility For KSC Ascent 
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criteria, and the roles of the top and bottom antennas have 
been interchanged, since the Orbiter flies "upside down" dur- 
ing ascent. The solid portions of the wide bars thus imply 
visibility to the top antenna. The cross hatched portions of 
the wide bars imply visibility to the bottom antenna without 
regard for blockage by the external fuel tank. Thus, these 
cross hatched segments should be regarded as periods of no 
visibility, since the tank does block the bottom antenna. 

Figure 5.3-6 shows the excellent ground transmitter 
coverage for ascent to a 40 deg inclination orbit. Figures 
5.3-5 and 5.3-7 show that coverage reduces for extreme 
launches, especially for the eastward launch. Generally speak- 
ing, to support accurate navigation when it is needed during 
ascent would require three signals available during roughly 
the five to seven min period following liftoff. The Montauk 
Point, Cape Hatteras, Bermuda, and Abaco Island sites would 
probably cover all launch azimuths voith marginal capabilit} 1 
(assuming no satellite visibility). 

The eastward launch suffers the poorest ground trans- 
mitter coverage. Figure 5.3-8 shows how this coverage would 
dramatically improve if the receiving antenna beamwidth could 
be enlarged to full hemispherical coverage. 

Figures 5.3-9 and 5.3-10 show what GPS satellite visi- 
bility might be available to the 40 and 28 deg inclination 
ascents with the 24 satellite, Phase III constellation. Unlike 
the case with ground transmitters, satellite visibility varies 
depending on the time of day, being roughly similar in each 
eight hour period. (It takes the earth eight hours to rotate 
the 120 deg separation between the three GPS satellite orbital 
pianes.) Thus, the figures give two time variables, one for 
time since liftoff in minutes, and one for time of day (called 
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Figure 5-3-8 Ground Transmitter Visibility For KSC Ascent 

(28 Deg Inclination Orbit) 

(90 Deg Half-Cone Antennas) 

offset time) in hours. The figures show the number of satel- 
lites visible to the (80 deg half-cone) top antenna during 
ascent, with a blank space indicating no satellites visible. 

Figure 5.3-9 shows that for a northward launch to a 
40 deg orbital inclination no satellites are ever visible after 
about four min after launch. Satellite visibility for a more 
northward launch would be even worse; thus, the 55 deg results 
are not shown. Clearly, no satellite visibility can be expected 
in the critical five to seven minute period following launch. 
Figure 5.3-10 shows an improving situation for an eastward 
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launch. Still, there are significant portions of each day 
when no satellites may be seen during the critical period, 
even though one or two are often available. 

Travelling at high speed close to a ground transmitter 
site produces a large range acceleration that puts severe de- 
mands on the receiver tracking loops. Figure 5.3-11 shows a 
time history of range acceleration for the ascent trajectory 
to the 55 deg orbital inclination and the Wallops Island ground 
transmitter site. The curve peaks out at over 13 g's and spends 
most of the period of visibility higher than three g’s. The 
receiver tracking loops will have to accomodate three g's to 
account for Orbiter acceleration levels, but it would strongly 
affect receiver design to demand higher levels. 


A-A6731 



Figure 5.3-11 Example Plot of Range Acceleration During 

Ascent 
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5.4 CONCLUSIONS OF THE GROUND TRANSMITTER STUDY 


Several conclusions summarize the results of the two 
ground transmitter studies for ascent and entry. They are: 


Satellite Visibility 

During entry, satellite visibility will be ade- 
quate once the full Phase 111 GPS constellation is estab- 
lished. Substantial periods of no visibility will exist 
during Phase II. During ascent, northward launches can ex- 
pect no satellites visible to the baseline antenna system 
even after the Phase III constellation is in place. Eastward 
launches may occasionally have one or two satellites in view. 


Ground Transmitter Visibility 

Three ground tiansmitters would be adequate to 
support any particular entry trajectory, with four or five 
being necessary to support a range of trajectories to one 
landing site. An additional ground transmitter would be 
required if a quartz crystal clock were used in the R/PA 
instead of an atomic clock. Four ground transmitters 
would provide adequate support for ascents from Kennedy 
Space Center. 


Ground Transmitter Power Levels 


In general , it would not be a problem to select 
constant power levels for ground transmitters that would 
not interfere with other GPS signal sources for the Space 
Shuttle during either ascent or entry. Some care would 
have to be used with a transmitter at the landing site for 
entry. Signal interference to other (aircraft) users of 
GPS is a possibility. 


Ground Transmitter Signal Dynamics 

Signal dynamics for ground transmitter signals 
during either ascent or entry are severe. Receivers de- 
signed for the baseline three g's of signal acceleration 
would generally be unable to use ground transmitter signals. 
Design of receiver tracking loops not aided by IMU informa- 
tion and capable of tracking at the lOg's and higher levels 
necessary is very difficult and would have a major program 
impact . 
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Alternatives to Ground Transmitters 


During entry, there appears to be no alternative 
to ground transmitters other than to ’’give up" on GPS 
navigation until the Phase 111 constellation is suffi- 
ciently established to provide adequate satellite vis- 
ibility, probably in 1986 or 1987. During ascent, there 
are likewise no alternatives until Phase III is established. 
Once the satellites are in place, alternatives to ground 
transmitters include widening the topside antenna field 
view, mounting a GPS antenna on the external fuel tank, 
or changing the attitude profile of the Orbiter so the 
satellites become visible to the unshaded topside antenna. 
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6. GPS RECEIVER SUBSYSTEM ANALYSIS 


The objective of GPS is navigation, but its methods 
are those of communication, and there are several GPS user 
equipment communication issues that impact navigation capabil- 
ity. The user equipment communication issues may be broken 
into two general areas: signal acquisition/reacquisition and 

signal tracking. Each of these areas is treated in this chapter. 


6.1 GPS SIGNAL ACQUISITION/REACQUISITION 

This section discusses the issues involved in acquiring/ 
reacquiring a GPS signal and the time required to establish a 
navigation fix. The discussion includes an explanation of the 
need for an indirect method of acquiring the P-code, as well 
as the functional dependence of time- to-acquire/reacquire on 
signal- to-noise density, C/N Q , and on residual line-of-sight 
dynamics. The nature and complexity of the GPS acquisition 
process results from the GPS signal structure. The relevant 
properties of the signal structure are discussed in the next 
section, followed by a discussion of the actual acquisition/ 
reacquisition process including navigation fix times. 

6.1.1 GPS Signal Structure 

The Global Positioning System (GPS) Navigational Sig- 
nal is a composite waveform consisting of a Protected (P) Sig- 
nal and a Clear/Acquisition (C/A) Signal transmitted in phase 
quadrature. The C/A-signal is designed for commercial and 
other users not requiring the highest accuracy and desiring to 
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use low cost equipment. From a military standpoint, its pri- 
mary function is in providing a method of acquiring the pro- 
tected code. The protected code is designed for both the 
highest navigation accuracy and for high anti- jam protection. 
The C/A-signal, by design, is easy to acquire but is highly 
susceptible to jamming. On the other hand, the P-signal, 

JL 

having been designed to resist jamming” , is nearly impossible 
to acquire directly without a priori range information; hence, 
the need for an indirect acquisition procedure. 

Each signal transmitted by a given satellite (P and 
C/A) consists of a carrier wave modulated by both 50 Hz data 
and a bi-phase digital ranging code. The data consists of 
space vehicle ephemerides, system time-of-day, C/A to P hand- 
over information^, etc. The bi-phase digital codes are de- 
rived from a class of digital codes known as pseudo-random 
noise sequences (PRN sequences). The spread spectrum nature 
of the codes makes them ideal for highly accurate ranging while 
providing high anti- jam (A/J) protection. 

The PRN sequence used for the C/A ranging code is 

composed of -1000 nsec (300m) chips with a code repetition 

period (i.e., code length) of 1023 chips (1 msec). Whereas, 

the PRN sequence for the P-code is composed of -100 nsec (30m) 

12 

chips with a code repetition period of 6x10 chips (7 oays). 
The narrower chip width gives the P-code a factor of ten im- 
provement in ranging accuracy over the C/A code. The longer 
code period gives the P-code added protection against jamming 
while denying its use to unauthorized users. The sequences 


*The P-code has a seven day period which, while making direct 
acquisition quite difficult, is the very feature which gives 
the P-code its anti-jam resistance. 

tThe Handover information enables acquisition of the P-signal 
with a minimum of search time. 


6-2 



THE ANALYTIC SCIENCES CORPORATION 


are similar in that they have very sharp autocorrelation func- 
tions: i.e. , if a PRN sequence is correlated with a replica of 
itself that differs by more than a chip, then the resulting 
signal has essentially zero power. Thus, for a GPS receiver 
to detect a GPS RF signal, it must correlate the incoming sig- 
nal with a replica of the code delayed by less than a chip. 
Otherwise, the effective signal- to-noise ratio, SNR, is = -30dB 
for the P-code (-20 dB for the C/A-code). This compares with 
a SNR as great as 40 dB for perfect correlation. Therefore, a 
transmitted signal modulated by a PRN sequence is inherently 
nonobservable to a receiver that does not know the repetition 
time of the PRN sequence. 

6.1.2 GPS Signal Acquisition and Reacquisition 

The term acquisition denotes the synchronization of 
carrier and code phases of the user GPS receiver with the sat- 
ellite transmitted P-code. The acquisition procedure is there- 
fore the first step in establishing a navigation fix. The 
user must demodulate the GPS data after closed loop tracking 
is achieved. This demodulated data is then decoded and com- 
bined with pseudo range and pseudo delta range measurements in 
an appropriate algorithm before computing a navigation fix. 

The time required to perform this entire operation is called 
the time-to-first-fix , TTFF, and is composed of the acquisi- 
tion time, the data demodulation time, tracking time, and the 
navigation fix time. 

The structure of the GPS signal is such that two dis- 
tinct methods of acquiring it are available to the user: direct 

and normal. For the direct method , a priori information is 
utilized to define a time- frequency uncertainty region for the 
P-signal. This region is then searched to determine the sig- 
nal's code position and carrier frequency so that tracking may 
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be initiated. Because the P-signal is designed to be secure, 
direct acquisition may require unsatisfactorily long acquisi- 
tion times unless accurate a priori information is available 
to the user. 

For the normal method , a priori information is em- 
ployed to define a time- frequency uncertainty region for the 
C/A- signal . This region is then searched to determine the 
code position and carrier frequency so that tracking can begin. 
The tracking operation (for the C/A-signal) and associated 
data demodulation continues until the handover word (which is 
present in the data signal) is recovered -- handover words are 
present in the GPS data every six seconds. This word provides 
the information required to acquire the P-signal directly in a 
reasonable length of time. A brief, direct acquisition of the 
P-signal is the final step in the normal method. When tl &. a 
priori information available to the user is not extremely 
accurate, the r rmal method will provide substantially better 
acquisition times than the direct method. 

There are several factors which affect the acquisi- 
tion time of a receiver. An obvious factor is the quality of 
the a priori information about range delay and doppler shift 
of the GPS signal. This factor is determined by uncertainties 
in the position and velocity of the user relative to the sat- 
ellite, and by uncertainties in the time and frequency of the 
user's clock. Other factors influencing acquisition time in- 
clude the received SNR, the desired probability of a correct 
acquisition, the search pattern used to check each cell in the 
time- frequency uncertainty region, and the structure of the 
acquisition receiver. 

The a priori range-velocity uncertainty defines a code 
posi tion- frequency uncertainty region which must be searched 
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by the receiver". The acquisition procedure is initiated by 
testing a code position and a doppler frequency. This test is 
accomplished by replicating the received signal as it would 
exist under the assumption that the test parameters are cor- 
rect. In the simplest case, this signal is then correlated 
with the actual received signal, bandpass filtered, and then 
signal power detected. If the assumed doppler frequency is 
incorrect, the mixed signal will not pass the bandpass filter. 
Similarly, if the assumed code position is incorrect by more 
than a chip, the SNR will be -20dB for the C/A-code (-30 dB 
for the P-code). If the test parameters are close to the cor- 
rect values, the SNR will be as great as AO dB. If the cross- 
correlation of the received signal and the locally-generated 
signal estimate (with the assumed code position and doppler 
frequency) do not pass a predetermined energy threshold, those 
signal parameters are rejected, and the local signal estimate 
is sequentially stepped through 'the uncertanty region until 
the energy exceeds the threshold. When the threshold is ex- 
ceeded, operation is switched to the tracking mode. 

Assuming a low cost GPS set in which all tests must 
be performed serially', Fig. 6.1-1 depicts the time to acquire 
a C/A-signal for a 300 meter range uncertainty as a function 
of doppler uncertainty and C/N Q , i.e., GPS signal strength at 
the receiver input. The results assume probability of detec- 
tion, P^=0.9, and probability of false alarm, Pp=10~^. Each 
curve within the family of 5 curves is composed of three regions 
designated by a solid curve, a dashed curve, and a dotted curve. 


*lf previous receiver operation has not provided GPS system 
time, then an additional range uncertainty equal to the time 
error must be included. This is typically small compared 
with the range uncertainty. 

tThis assumes there are not multiple channels available which 
could be employed to perform testing in parallel. 
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C/N 0 (dB-Hz) 


Figure 6.1-1 Time to Acquire/Reacquire the C/A-Code as a 

Function of C/N , ar.r* Position and Velocity 
Uncertainties 


The dashed and dotted curves employ slight variations of the 
acquisition procedure defined above in order to optimize acqui- 
sition performance under weak signal conditions . For a 600 
meter uncertainty, the times indicated in the table double, 


*The solid curve assumes only a single doppler frequency is 
searched. As C/N decreas a sequence of independent samples 
at .he output of the defector are required as the basis for 
detection, (dashed curve). The number of samples available 
depends on how long the doppler accuracy will keep the test 
code within a half chip of the received code. A further de- 
crease in C/N requires multiple doppler regions to be searched 
in addition to multiple code position, (dotted curve). 
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etc. It should be noted that, in the worst case, no a priori 
range information is required for C/A acquisition; i.e., since 
the signal is periodic, the time uncertainty is limited to the 
code period, 1023 msec. Thus, with no a priori range informa- 
tion a maximum 1023 chips must be tested (and typically codes 
are tested every 1/2 chip) for a total number of 20^6 tests. 

In Fig. 6.1-1 a scale has been added for the total acquisition 
time (right hand side) for one C/A channel as a function of 
C/N q assuming no a priori range information. 

Figure 6.1-2 is analogous to Fig. 6.1-1; it depicts 
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Figure 6.1-2 Time to Acquire/Reacquire the P-Code as a 

Function of C/X , and Position and Velocity 
Uncertainties 
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the time to acquire a P-signal for a 300 meter range uncer- 
tainty as a function of doppler uncertainty and C/N 0 . Again, 
the times double for a 600 meter uncertainty. Since the P- 
code has a 7-day period, there is no corresponding right hand 
scale depicting acquisition time if no a priori range infor- 
mation is known; i.e., the time would be prohibitive. 

These acquisition times can be combined with the data 
demodulation time, tracking time, and the navigation fix time 
to provide the time-to-first-fix (both normal and direct). 

The data demodulation time is the time required to demodulate 
the entire data block (1500 bits) in order to recover the 
necessary satellite data for utilization in the navigation 
algorithm. Since the data is transmitted at a 50 Hz rate, 
this operation requires exactly 30 seconds following bit syn- 
chronization (<1 sec). The tracking time governs the accuracy 
of pseudo range and pseudo delta range measurements. Its value 
depends on the receiver bandwidths; values of 2 to 6 seconds 
are reasonable. Finally, the navigation fix time, which de- 
pends upon the user's computation capabilities and the fix 
algorithm, will typically require approximately 1 second. 

Table 6.1-1 lists the times associated with the indi- 
vidual factors that make up the time-to-first-fix for the 
P-signal. Several points should be noted: 

• If the range-doppler uncertainty is 
small enough to allow P-code acquisition 
in less than -8 seconds (see Fig. 6.1-2), 
then direct acquisition is faster than 
normal acquisition 

• The time required to demodulate the data 
block (30 seconds) can make up the major 
portion of time-to-first-fix, TTFF 
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TIME 6.1-1 

TIME-TO-FIRST-FIX FOR THE P-SIGNAL 


ACQUISITION PHASE 

NORMAL ACQUISITION 

DIRECT ACQUISITION 

Acquire Carrier and C/A-Code 
Phase (Depends on Position and 
Velocity Uncertainties) 

t c/a, acq 

(Fig. 6.1-1) 

— 

Obtain Bit Synch and Frame Synch 

< 1.0 Sec 

— 

Demodulate and Decode Handover 
Word 

6.0 Sec 

— 

Acquire P-Code 

< 1.0 Sec 

T P, ACQ 
(Fig. 6.1-2) 

Demodulate Data Block 
(1500 Bits @ 50 Hz) 

30 Sec 

30 Sec 

Derive Pseudo Range and Pseudo 
Delta Range Measurements 
(Depends on Receiver Bandwidths) 

T 

Trackingf 

T 

TrackingT 

Sum of Previous Times 

T SUB 

T SUB 

N=4 in Low Cost Sequential Set 
(N=l in 4 Channel Set) 

W 

t sub x n 

Perform Navigation Fix (Depends 
on User Computer and Fix 
Algorithm) 

* 1 Sec 

*v 1 Sec 


+T-. . . Governs the Accuracy of the Correlation Process and Therefore 

1 * the Measurement Accuracy (Values of 2 Sec to 6 Sec are 
Reasonable) . 


• A low cost sequential set (1 channel 
set) affects TTFF by introducing a 
factor of four increase in TTFF* rel 
ative to a four channel set. 


*Less the time required to perform the navigation fix. 
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An example rnry help to clarify the discussion. Assume the 
signal- to-noise density (C/N Q ) is 30 dB-Hz, the total velocity 
uncertainty is ±150m/sec, the total range uncertainty is 3000 
meters (~2 miles), and the user has a single channel receiver. 
Then, assuming also 4 seconds to derive pseudo range and pseudo 
delta range, both normal and direct acquisition would take ~169 
seconds; 120 seconds of which are required to demodulate the data 
block. The direct approach would be faster if C/N c increased or 
if doppier uncertainty and/or range uncertainty decreased. 

The term reacquisition denotes the resynchronization 
of carrier and code phases of the user GPS receiver with the 
satellite transmitted P-code that must be accomplished whenever 
the tracking loops lose lock. Because closed loop tracking 
has already been achieved, demodulated GPS data is available 
to the user attempting reacquisition. Since GPS data changes 
slowly, this available data may be employed in the navigation 
algorithm once closed loop tracking has been re-established. 
Consequently, when the outage is brief, the normal reacqui- 
sition-f ix-time would be the acquisition time of the P-signal 
(Fig. 6.1-2) plus ~4 sec to derive pseudo range and pseudo 
delta range measurements, plus the navigation fix time. (This 
value is multiplied by four if a sequential single channel set 
is used.) Reacquisition-fix-time is naturally shorter than 
the time-to-first-fix since the necessity of ^ata demodulation 
has been eliminated and since the position and velocity uncer- 
tainties are small. When outage is sufficiently long, faster 
reacquisition can be accomplished by first locking-up the C/A- 
signal and then reacquiring the P-code by the handover word; 
this is referred to as the indirect reacquisition method. 

Table 6.1-2 lists the times associated with the factors that 
make up the reacquisition- fix- time for the P-signal for both 
the normal and indirect reacquisition methods. This table 
differs from Table 6.1-2 primarily in there being no time re- 
quired to demodulate the data block. 
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TABLE 6.1-2 

REACQUISITION- FIX-TIME FOR THE P-SIGNAL 


REACQUISITION PHASE 

NORMAL 

REACQUISITION 

INDIRECT 

REACQUISITION 

Acquire Carrier and C/A-Code Phase 
(Depends on Position and Velocity 
Uncertainties) 

• • • 

T C/A, ACQ 
(Fig. 6.1-1) 

Obtain Bit Synch and Frame Synch 


< 1 Sec 

Demodulate and Decode Handover Word 

— 

6 Sec 

Acquire P-Code 

T P, ACQ 
(Fig. 6.1-2) 

< 1 Sec 

Derive Pseudo Range and Pseudo Delta 
Range Measurements (Depends on Receiver 
Bandwidths) 

T 

Tracking! 

T 

Tracking! 

Sum of Previous Times 

T SUB 

t sub 

N=4 in Low Cost Sequential Set 
(N=l in 4 Channel Set) 

t sub xS 

t sub xs 

Perform Navigation Fix (Depends on 
; User Computer and Fix Algorithm) 

* 1 Sec 

^ 1 Sec 


TT t . . Governs the Accuracy of the Correlation Process ana Therefore 
* the Measurement Accuracy (Values of 2 to 6 Seconds are 
Reasonable) 

6.1.3 Acquisition/Reacquisition Summary 

The complexity of the acquisition/reacquisition proc- 
ess results from the structure of the GPS signal; namely, un- 
less both the pseudo-random code modulating the GPS signal as 
well as the time when the code repeats are known to a receiver, 
the GPS signal is inherently non-observable; i.e., buried in 
the noise. To account for the GPS signal structure, there are 
two methods available for performing both acquisition and re- 
acquisition: the normal and indirect method for reacquisition. 
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The time required to first establish a navigation fix 
is referred to as time-to-first-fix , TTFF. TTFF for the 
P-signal is tabulated in Table 6.1*1 for both the normal and 
direct methods of acquisition. The normal method requires 
acquisition of the C/A signal; the acquisition time is shown 
in Fig. 6.1-1 in terms of C/N Q and position and velocity un- 
certainties. The direct method of acquisition requires acqui- 
sition of the P-signal; this acquisition time is shown in Fig. 
6 . 1 - 2 . 


When the GPS receiver loses lock, resynchronization 
of the carrier and code phases with the transmitter P-signal 
must be accomplished before the next navigation fix can be 
obtained. The time required to obtain the navigation fix is 
referred to as reacquisition- fix- time . Reacquisition- fix- time 
for the P-signal is tabulated in Table 6.1-2 for both the nor- 
mal and indirect methods of reacquisition. The normal method 
requires reacquisition of the P-signal; this time is given in 
Fig. 6.1-2 as a function of C/K T 0 and position and velocity 
uncertainties. The indirect method of reacquisition requires 
reacquisition of the C/A-signal; this time is given in Fig. 
6 . 1 - 1 . 


The method used to acquire/reacquire a GPS signal 
will depend on the range and velocity uncertainties as well as 
on the signal- to-noise density, C/N o - lr particular, if uncer- 
tainties and C/K q are such that the P-code can be acquired 
within 6-8 seconds, then the direct acquisition process will 
be faster than the normal method . Similarly, if the P-code 
can be reacquired within 6-8 seconds, then the normal reacqui- 
sition process will be faster than the indirect method . 
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6.2 GPS SIGNAL TRACKING ANALYSIS 

The analysis reported in this chapter was done when 
the phased approach to putting GPS aboard the Shuttle was being 
considered. The most attractive candidate for a receiver for 
the experimental phase is called GPSPAC, under development by 
Magnavox for use in satellites. While the phased development 
approach and the experiment plan have been abandoned, the re- 
ceiver analysis results of this chapter are presented because 
they are still relevant. In fact, it is possible that a deriv- 
ative of the GPSPAC could still be a candidate for the Shuttle 
GPS receiver. 

The baseline navigation system for the Space Shuttle 
will require measurements from a variety of navigation aids in 
order to maintain accurate position and velocity estimates, 
including navigation updates from a ground tracking system 
during orbital operations. One of the greatest potential bene- 
fits of GPS to the Space Shuttle is the possibility of obtain- 
ing accurate, autonomous navigation across the entire flight 
regime from a single navigation aid. Unfortunately, none of 
the GPS user equipment under development is versatile enough 
to fulfill all potential Space Shuttle navigation requirements 
without either software or hardware modifications. However, 
the GPSPAC receiver being developed by Magnavox is sufficient- 
ly attractive in terms of both its operational flight envelope 
and its potential availability to warrant a more detailed eval- 
uation. The results of a preliminary analysis of the GPSPAC 
receiver design are presented in this chapter with emphasis on 
the changes required to meet the Space Shuttle navigation 
requirements . 

The five modules present in most user equipment de- 
signs are identified in Fig. 6.2-1. The heart of the equipment 
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R- 25364 


AUnilARV SENSORS 
(IMU. ALTIMETER) 


PRIME 

POWER 


MAN /MACHINE 
INTERFACE 

Figure 6.2-1 Major Elements of GPS User Equipment 

is the receiver subsystem which tracks the GPS signal and gen- 
erates range and range- rate measurements. In addition to the 
measurements, the receiver subsystem outputs system data trans- 
mitted by the satellites. In a high dynamics environment, the 
receiver tracking loops are generally rate-aided with processed 
IMU outputs generated by the data processor. The data proces- 
sr ” also incorporates the GPS measurements into the navigation 
solution. In most user equipment designs, this is accomplished 
with a Kalman filter which is also capable of processing meas- 
urements from auxiliary sensors such as an altimeter or IMU. 

The GPSPAC receiver is a two channel, space-rated set 
designed for a low dynamics environment. Measurements from four 
different satellites are obtained by sequencing between the sat- 
ellites at well-defined intervals. The performance specifica- 
tions require the receiver subsystem to maintain signal lock 
at orbital velocities and in the presence of accelerations as 
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large as lb m/sec . Rate-aiding would not be possible without 
a software modification. The general processor uses a Kalman 
filter, out IMU data could not be incorporated into the navi- 
gation solution without a software me dification. 

It is apparent that the GPSPAC receiver subsystem as 
presently designed could operate satisfactorily during orbital 
operations for the Space Shuttle. The remainder of this chap- 
ter investigates potential receiver performance during re-entry. 

6.2.1 Receiver Subsystem 

Much of the nature and complexity of GPS receiver 
structures derives from the nature of the GPS signal. General- 
ly, the satellite signals consist of RF carrier waves bi-phase 
modulated (±90°) with pseudo-random codes and 50 baud (bit-per- 
second) data. The carrier is tracked for user- to-satellite 
relative velocity information. The codes are tracked to pro- 
vide user- to-satellite pseudo ranges (i.e., actual range plus 
user clock bias error) . The 50 baud data has information con- 
cerning satellite ephemerides, etc. 

Although there can be substantial variation in the 
design of an appropriate receiver, some essential elements are 
common to most proposed designs. Noncoherent delay-locked 
loops are used for locking to the pseudo-random code trans- 
mitted by each satellite, and Costas carrier tracking loops 
are used for locking to the doppler-shifted suppressed carrier. 
The pseudo ranges obtained from the code tracking provide the 
positional fix, while the range-rate estimates obtained from 
the Costas loop provide the velocity fix. The carrier phase 
reference obtained from the Costas loop is also used in the 
data demodulation, which is accomplished using 20 msec inte- 
grate-and-dump circuits. 
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The proposed GPSPAC receiver (Ref. 30) utilizes a 
single IF (Intermediate Frequency) and VCO (Voltage Controlled 
Oscillator) which is common to both the carrier and code 
tracking loops. Figure 6.2-2 provides a block diagram of the 
receiver subsystem. The common IF design is accomplished by 
time multiplexing the two tracking loops. In particular, the 
pseudo-noise (PN) demodulator (1 in the figure) is time shared 
between two distinct modes. In mode 1, the PN demodulator 
mixes an on-time version of the locally generated P-code with 
the incoming IF signal. This, in effect, provides a recon- 
structed carrier which is utilized by the Costas loop for 
phase tracking . In mode 2, the PN demodulator alternately 
mixes the input signal with an early and late version of the 
P-code. This results in a code phasing error signal which is 
used by the delay- locked loop for code phase alignment. The 
common IF design eliminates dedicated code channel hardware . 
However, as will be discussed, dynamic errors can result while 
the PN demodulator is in mode 2 and the Costas loop is coasting . 

The VCO (3 in the figure) is used in a Costas- type 
phase-locked loop to track the carrier signal phase. The code 
loop clock is also derived (A in the figure) from the same 
carrier tracking VCO. The common VCO design eliminates the 
effect of vehicle dynamics on the code loop ; i . e . , when the 
Costas loop is tracking the carrier (mode 1) the code clock is 
inherently locked in frequency to the received code. However, 
it remains to align the code phase (mode 2) to derive the range 
measurement . 


*ln actuality, the carrier is still suppressed by the data; 
however, the data is removed by subsequent inphase/quadra- 
ture multiplication (2 in the figure). 
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The receiver mode of operation is first to establish 
carrier tracking while the PN demodulator is in mode 1. Having 
locked to the received carrier signal, the PN demodulator then 
switches to mode 2 during which time code phase alignment is 
accomplished. (The Costas loop coasts via a prediction of the 
carrier signal during this period.) Once aligned, the PN de- 
modulator may switch back to Costas loop operation (mode 1) to 
remove any carrier tracking errors (due to vehicle dynamics) 
which may have occurred while in mode 2. 

A model for the Costas loop is shown in the upper 
half of Fig. 6.2-3. The loop provides the coherent phase ref- 
erence for data demodulation and range-rate determination. 
Because of its improved acquisition performance compared with 
that of a third-order loop, the Costas loop employs a second- 
order loop (35 Hz bandwidth, 0.707 damping coefficient). 
However, use of a second-order loop can result in sizeable 
phase tracking errors during acceleration unless some form of 
rate-aiding is employed . (See Section 6.2.2) 

A model for the code loop is shown in the lower half 
of Figure 6.2-3. The loop provides code phase alignment for 
range determination. The code loop clock is derived from the 
carrier tracking VCO by dividing the VCO output (consisting of 
some IF frequency plus carrier lcop doppler) by the code/ 
carrier ratio (154) to obtain the doppler-shifted code fre- 
quency. This signal is then mired with a reference to. produce 
the 10.23 MHz P-code clock which in turn is used during mode 2. 

During mode 2 a code phase error signal is derived by 
alternately modulating the input signal with an early and late 
version of the P-code. This technique is referred to as x- 
dithering and has the major advantage over the early-late 



IF ft CARRIER LOOP DOPPLER 



Kigure 6.2-3 Simplified Receiver Block Diagram 



tmE ANALvTIL, SCIENCES CORPORATION 


delay-locked loop approach that it employs only one correlation 
channel. Thus, this implementation minimizes hardware , while 
eliminating the problems associated with two channel gain im- 
balances. However, the time sharing of the early and late 
channels does result in an approximate 3 dB degradation in the 
code loop signal-to-noise ratio. In so much as the system 
need not consider ultra low signal-to-noise operation (i.e., 
high anti-jam operation), this seems to be a reasonable trade- 
off. The code phase error signal consists of an amplitude 
modulation of *■»!€ carrier occurring at the dithering rate (25 
Hz). The amplitude difference of the two components (early 
and late) is proportional to the code phase error; whereas, 
the sign of the twc components (relative to the r -dither ref- 
erence signal) contains the polarity of the code phase error. 

The code phase error signal is detected by the code 
loop during mode 2. The baseband envelope detected signal is 
passed through an up/down integrator which is controlled by 
the t - dither reference. The polarity of the integrator output 
indicates whether a positive or negative phase adjustment to 
the code loop clock is required. The magnitude of the inte- 
grator output is input to a VFC (Voltage to Frequency Converter). 
The VFC output provides a phase correction to the code loop 
clock. This correction is applied at a fixed rate; the r**te 
defines the code loop bandwidth. 

6.2.2 Error Analysis 

The essent'al elements which make up the GPS receiver 
are a Costas carrier tracking loop used for locking to the 
doppler shifted carrier, and a noncoherent delay-locked loop 
for locking to the pseudo random code transmitted by each sat- 
ellite. The performance of each loop can be specified in terms 
of rms phase and delay jitter, respectively. In turn, each is 
a function of two components, input noise and user dynamics. 
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The phase- iocked loop requires a higher signal-to- 
noise density ratio, S/N Q , than the delay-locked loop to main- 
tain lock. Thus, one need only investigate the performance of 
the phase- iocked loop to determine the minimum S/N Q of the 
receiver. (This assumes that a phase- locked loop is needed to 
provide a coherent reference for data demodulation. Otherwise, 
if data demodulation is not required, the AFC frequency acqui- 
sition function of the GPSPAC receiver -- which can operate in 
a lower S/N Q than the Costas loop -- will dictate the minimum 
S/N 0 . ) 


With regard to user dynanics, the receiver configura- 
tion has eliminated the effect of vehicle dynamics on code 
loop performance by employing a common VCO design. Thus, user 
dynamics impact the performance of the GPS receiver by affect- 
ing the Costas loop. Coupled with the previous observation on 
minimum S/N Q , this suggests that preliminary analysis need 
only be concerned with the performance characteristics of the 
Costas loop. A complete, rigorous investigation into the per- 
formance of the receiver should be performed at a later time. 

The total allowable rms phase error of the Costas loop 
depends on the receiver function under consideration. A pre- 
defined range-rate accuracy will define a maximum rms phase er- 
ror; this allowable error is referred to as the track threshold . 
However, independent of the track threshold, the loss-of-lock 
threshold defines an inherent upper bound on the rms phase er- 
ror; loss-of-lock threshold is reached at an rms phase error 
of = 0.25 radians. (The Costas loop implementation tracks 
a double frequency carrier which would then have a 0.5 radian 
rms error; 0.5 radian is considered the loss-of-lock threshold 
for a phase- locked loop.) 
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The total rms error of the Costas loop is composed of 

two terms. The first, a . , is a result of input noise, and 

*n 

the second, o A , is due to user dynamics. The total rms error 


is given by 


o 






(6.2-1) 


As is shown in the following section, for a fixed S/N c narrow 
ir.g the loop bandwidth will tend to impro' e tracking perform- 
ance by decreasing the effective noise. However, the useable 
bandwidth is bounded from below by the user's dynamic motion; 
this is discussed in Section 6.2-2. 


Oscillator instability errors have also been investi- 
gated, and, for the nominal Costas loop bandwidth, the rms 
phase error was found to be on the order of 1 mrad. With re- 
gard to vibration- induced fluctuations, these need not be a 
serious problem in the system under consideration. Note, in 
particular, that the Costas loop noise bandwidth is rather 
wide (B^ = 35 Hz) and phase fluctuations below the loop cutoff 
are tracked out by the loop. If necessary, it w’ould be a rela- 
tively simple matter to provide damping and isolation to a 
small oscillator for mechanical vibration frequencies above, 
say, 10 Hz. 

Noise Error - The Costas loop error variance is of 

the form 
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where is closed-loop single-sided bandwidth, S/N Q is the 
equivalent signal-to-noise density assuming all correlation 
losses and other factors, and is defined to be the squaring 
loss associated with a squaring type phase-locked loop (such 
as a Costas loop). The squaring loss would be equal to unity 
for the common phase- locked loop. However the GPS data modu- 
lation results in a suppressed carrier which can only be tracked 
by a squaring- loop method. 

The squaring loss is defined by 

s l " f 1 + k l r “ini (6.2-3) 


where is a function of the type of low-pass filter employed 
(see Fig. 6.2-3) and B^ p is the 3 dB single-sided bandwidth of 
the low-pass filter. is equal to 5/6 assuming the low-pass 
filters are third-order Butterworth. Thus, the noise variance 
is given by 


2 _ 


N 




5“ B L 


1 + 


N o B LP 
3 “ “ 2 “ 


(6.2-4) 


The value for B u> is constrained by the presence of 
the 50 baud d;.ta. Given that 


B Lp = 200 Hz 
B l = 35 Hz 


and 


o. unlock - 0.25 radians 
0 
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eq. 6.2-4 can be solved for the minimum equivalent signal -to- 
noise density in the absence of other disturbances. This is 
roughly (S/N 0 ) m ^ n = 28 dB-Hz; Fig. 6.2-4 plots o^ versus S/N 0 . 

Thus, the minimum signal-to-noise density at the receiver in- 
put, S/N rt (allowing for such factors as cable losses) 

o input ^ 

is somewhere in the vicinity of 30 dB-Hz. In the presence of 
dynamic errors, larger values of S/N 0 are required to maintain 
lock; this is discussed in the following section. 



Figure 6.2-4 RMS Phase Error Versus Signal- to-Noise Density 


Dynamic Errors - The phase error for a phase-locked 
loop ( in the frequency domain) is given by the expression, 


‘''Given a 35 Hz filter bandwidth, these values can be converted 
to the more familiar signal-to-noise ratio by subtracting 
15.5 dB . 
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S0 i (s) 

e e (s> 1 s ” K’RJTTs! 


(6.2-5) 


where F(s) is the carrier loop filter (Fig. 6.1-2), is the 
phase detector gain, K Q is the VCO gain constant, and 0^(s) is 
the Laplace transform of the input phase. 

If the input is a step change in frequency of magni- 
tude Aw (equivalent to a step change in velocity), then 0.(s) 

2 1 

is a phase ramp; i.e., 0^?) = Aw/s . The steady-state error 
is given by the final value theorem of Laplace transforms. 


lim y(t) = lim s Y(s) (6.2-6) 

t-+“> s-*0 


Thus , 


lim © e (t) = lim 
t-»00 S"*0 


Aw _ Aw 

s+ W (s > ' K v 


(6.2-7) 


Where K v is the dc loop gain and is given by the expression 

K v = K o K d F(0) (6.2-8) 

The resulting error is the steady state velocity error, o . 

*v 

For an active filter in a second-order loop, F(0) = A, where A 
is the dc gain of an operational amplifier. 

When the error is due to a doppler shift (and not to 
an actual difference between transmitter and receiver frequen- 
cies). o is given by 
% 
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(6.2-9) 


where r is the relative range-rate along the line-of-sight and 
c is the speed of light. The loop gain, K v , has not been spe- 
cified; however, it can be shown that a reasonable choice for 
the loop gain is on the order of 10^ to 10^. 

Assuming a value of K y = 5 * 10^, Fig. 6.2-5 shows 

the linear relationship between o ar.^ r. (For different 

*v 

values of K v the slope scales inversely with K v >) The maximum 
velocity uncertainty before the loop unlocks is obviously a 
function of S/N 0 (see Eqs. 6.2-1, 6.2-2). This relationship 
is shown in Fig. 6.2-6. For S/N 0 less than -28 dB-Hz the loop 
unlocks due to noise alone; i.e., independent of velocity un- 
certainty. The steady state analysis would indicate that as 
the signal level increases larger velocity errors can be tol- 
erated -- up to a maximum of 373 m/sec for the high S/N Q case. 
(This result is conditioned on there being no other dynamics 
involved and on a loop gain of 5 x 10^). Since velocity er- 
rors are expected to be on the order of 10 m/sec, Fig. 6.2-6 
implies that the velocity uncertainty has negligible impact on 

defining (S/N ) . . However, it should be noted that if such 
o min 

a large velocity step were applied, the loop would actually 
skip cycles for a while following the input transient, before 
locking up again. While the carrier loop is cycle slipping, 
velocity measurements are meaningless. The maximum velocity 
step which can be tolerated without cycle slipping is propor- 
tional to the carrier loop bandwidth; for a bandwidth of 35 
Hz, the velocity step must be less than -7 m/sec. 
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VELOCITY UNCERTAINTY (m/secl 


Figure 6.2-5 RMS Phase Error Due to Velocity Uncertainty 



Figure 6.2-6 Maximum Velocity Uncertainty Before 

Loop Unlocks Versus S/N q 


6-2 






THE ANALYTIC SCIENCES CORPORATION 


For sufficiently large gain the steady state velocity 
error can be made quite small. However, a second-order loop, 
will produce a steady state phase error due to acceleration 
even for infinite gain . In particular, for infinite gain the 
phase error (defined by Eq. 6.2-5) in a second-order loop can 
be written in the form 

s 2 0 . (s) 

BA s) = - 7 (6.2-10) 

e s^ + 2 £w n s + w * 

n n 

where u> n is the natural frequency of the loop and £ is the 
damping coefficient. If the input frequency is changing linear- 
ly at a rate Aui rad/sec^ due to an acceleration, then the in- 
put phase is 6^(s) = Aw/s . This leads to the steady state 
acceleration error , , 

o. = lim 6 (t) = lim — 7 (6.2-11) 

™a t-** s-»0 s" + 2£w s + w_ 

* n n 

* 

_ Aw 

2 

% 

The natural frequency, w n , can be related to the brndwidth of 
the loop by the expression 

w R = 2B l / ^( + Z^) t 0.2- 12) 


where the damping coefficient is given as £ = 0.707. Thus, 
substituting th^s expression into Eq. (6.2-11) 
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(6.2-13) 


where r is relative acceleration along the line-of-sight . 


Assuming a bandwidth of 35 Hz, Fig. 6.2-7 shows the 
linear relationship between o. and r. (The slope scales in- 
2 

versely with B^ ) . Again, the maximum allowable acceleration 
uncertainty before the loop unlocks is a function of S/N Q . 

This relationship is shown in Fig. 6.2-8. Notice that even 
for infinite S/N q the maximum acceleration uncertainty must be 
less than 32.5 m/sec (assuming no additional phase error due 
to velocity errors). The acceleration tolerance of the Costas 
loop can be increased beyond those values indicated by Fig. 
6.2-8 by either: 


• Rate-aiding the receiver 

• Increasing the loop bandwidth, 

• Replacing the second-order loop with a 

third-order loop. 

The phase error due to acceleration decreases with 

2 

increasing B^ since is proportional to 1/B^ . The disad- 
vantage of an increased B^ is primarily a decrease in noise 
tolerance. A B^ of 43 Hz seems to be a reasonable compromise. 
This would provide marginally acceptable loss-of-lock perform- 
ance in a 5 g environment. The advantage of a third-order 
loop is that it has zero steady state phase error due to un- 
compensated acceleration errors. However, the pull-in behav- 
ior when the VCO is not initally tracking the signal is less 
stable than that of a second-order loop. This disadvantage 
could potentially be offset by the AFC acquisition circuitry 
in the receiver (Fig. 6.2-2). 
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Figure 6.2-7 RMS Phase Error Due to Acceleration Uncertainty 


R-27633 



Maximum Acceleration Uncertainty 
Before Loop Unlocks Versus S/N Q 
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Figure 6.2-8 
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6.2.3 Summary 

The proposed GPSPAC receiver utilizes a single IF and 
VCO which is common to both the carrier and code tracking loops 
The common IF design is accomplished by time- sharing the two 
tracking loops. Further hardware is eliminated by generating 
the code phasing error signal using t - dithering. The advan- 
tages of this implementation are: 

• The common VCO design eliminates vehicle 
dynamics from the loop as well as mini- 
mizing hardware 

• t -dithering minimizes hardware while 
eliminating the problems associated with 
two channel gain imbalances. 

The disadvantages are: 

• The common IF design can result in errors 
due to vehicle dynamics during the time 
interval dedicated to the code loop 

• There is a 3 dB degradation in the code 
loop signal-to-noise due to x-dithering. 

The Shuttle performance specifications require sub- 
system operation in a 5 g environment. The dynamic error analy 
sis presented in Section 6.2.2 indicates the acceleration tol- 
erance of the GPSPAC in this environment is, at best, marginal. 
The acceleration tolerance can be increased by either: 

• Rate-aiding the receiver 

• Increasing the loop bandwidth 

• Replacing the second-order loop with a 

third-order loop. 


6-31 



THE ANALYTIC SCIENCES CORPORATION 


Increasing the loop bandwidth to A3 Hz would keep the 
Costas loop in lock for acceleration uncertainties as great as 
5 g while raising the minimum usable signal-to-noise by only 1 
dB. A third-order loop has zero steady states phase error' due 
to acceleration; however, its pull-in behavior is not as stable 
as that of a second-order loop. Further analysis is required 
to determine the most desirable alternative. 
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7. 


SUMMARY AND CONCLUSIONS 


This report addresses a number of system implementa- 
tion and performance issues associated with the possible use 
of GPS as a navigation aid for the Space Shuttle. While each 
of the preceding seven chapters has a separate conclusions 
section, the more important results and conclusions of this 
study are summarized here for convenience. 


7.1 NAVIGATION WITH PHASE I AND II GPS 

Depending upon the final schedules for both the Or- 
oiter and for the GPS system, there may be a period of several 
years when GPS user equipment is aboard the Shuttle but the 
fully operational, Phase III GPS constellation of satellites 
is not yet in place. The reduced coverage provided by the 
Phase II GPS constellation of six (or perhaps even fewer) sat- 
ellites places some serious limitations on the ability of the 
Shuttle to navigate with GPS. The visibility (satellite sig- 
nal availability) studies presented in this report are an im- 
portant prelude to the navigation accuracy studies. Visibility 
patterns depend upon receiving antenna coverage, as well as on 
masking by the earth, and vary considerably with Shuttle flight 
phase and, for the incomplete Phase I and II constellations, 
with "time of day." 

Figure 7.1-1 shows a typical visibility product for 
the reentry phase, assuming a six satellite constellation. 

The vertical axis is time during the flight phase, while the 
horizontal axis is essentially "time-of-day . " Since three or 
four satellites are necessary for reasonable navigation during 


7-1 




THE ANALYTIC SCIENCES CORPORATION 


e 

1 

X 

u 

O 


25h 


5 20 

at 

0 

S 

IM 

u 
Z 
in 

1 30 1— 


2 

S 35 


40 r - 


OFFSET TIME (hours) 

8 12 16 


3 - 
5 - 

10 - 

15 - 


L> 

< 


451- 

TOUCHDOWN* 


T 


T 


T 


20 


T 


R- 26U1 

24 


»*•*••*•< 


33334655555655*4433 

33344455554555444433* • • . 

33344A455 54 55 54444 3 33 

3334444555455544443333* . ** 
3334444455445544443333* * * * 
• **4444455445544444 3333* •* 

* * • 334444 544554444433333 * • 

* * * 3344444445444444333333* 

* *33334444444444445533333* 
• *33334443444444555553333* 
• *333334434444445555553333 

* *333334434433455555554333 
• *333334334433455555555443 

* *333334333333455655555443 
*•3333333333 
**33333333*3 
*•33333*33** 


> *33333* 
>*33333* 


>3333. 

3333 •« 

333 * * * 

• *•*••• 333 * * 

333 * 

••*•••• * 33 * • 






»3« 

»3« 





• 33* 

• 33* 
*33. 
*33. 

• 3* . 
*3*. 

3* . 


** 56665555444 3. 
•*566665554443. 


3466665555443 
3466665555443 
3*5666655544 3 3 ry 


•34566665554433. 

35*5 , gOTH 5 4 T §3 f 

344455665555433* • 
• 44455665555433* • 
*444556665554333* 
*344555665554333* 

• 344555665554333* 
• 334555665555333* 
*334555565555333* 
•333555565555433* 
*333555565555433* 
*333555565555433* 
*333555565555433* 
• 333455565555533* 
*333455566555533* 

• 333455566555433* 
• 333455566555433* 

• 333455556555433 * 
• 333455556555433. 
. 3334SSSS6SS5433. 
.33345555655543. • 
.33345555655543. • 
.33345555555543. • 
.33345555555543. • 


SJL 


KEY: 

BLANK 

3-6 


NO SATELLITES 

1 OR 2 SATELLITES 

INDICATED NUMBER OF SATELLITES 




A- 




Figure 7.1-1 


Satellite Visibility During the 
Deorbit Phase as a Function of 
Offset Time 


a dynamic flight phase, the figure illustrates that the early 
GPS constellations will be capable of supporting reentry only 
during a short portion of each day . 



7-2 



THE ANALYTIC SCIENCES CORPORATION 


On-orbit visibility patterns demonstrate a periodicity 
roughly comparable to the orbital period. Figure 7.1-2 shows 
a statistical summary for a particular orbit (at one time 
slated for mission OFT- 5) with a six satellite constellation. 
Three or more satellites (which will provide high accuracy 
orbital navigation) are available more than half of the time. 
Possible changes to the Phase II constellation to reduce it to 
four satellites would reduce the three-satellite visibility to 
about one third of the time. 



Figure 7.1-2 Percentage Visibility for Typical Orbit 

as Function of Offset Time 

Given one of the good visibility situations possible 
during reentry, the simulation results in Chapter 2 indicate 
that reasonable navigation accuracy is achievable . Figure 
7.1-3 shows a typical result, with errors growing to about 
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Figure 7.1-3 RMS Position Error versus Time for GPS- 

Aided Deorbit Navigation with Good 
Visibility Prior to Satellite Uploading 

13000 ft (lo) during the reentry blackout period during which 
no GPS signals are assumed available, but being well controlled 
at other times. Various navigation filter mechanizations were 
evaluated in this study, with an eleven state filter (incorpo- 
rating 1MU platform misalignments along with the usual vehicle 
vector position and velocity and the clock phase and frequency 
states) proving adequate . Drag update measurements (part of 
the current Shuttle baseline) could control the error buildup 
during blackout . 
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Orbital navigation performance, on the other hand* is 
generally quite good at all times, even with only a six satel- 
lite constellation . Figure 7.1*4 shows a typical result. 

Here too, errors grow during periods when no satellite signals 
are available, but they grew much more slowly and typically 
peak at less than 1000 ft (lo). Low level thrusting on orbit 
due to venting of liquids and to attitude jet usage is the 
dominant contributor to this error growth. During periods of 
good visibility (three or more satellites in view) orbital 
navigation errors are in the area of 100 ft (lo) in position 
and 0.2 ft/sec (lo) in velocity. 


•* m.j 



elapsed time (mm) 


Figure 7.1-4 Orbital Navigation: RMS Position 

Error Versus Time 
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Surprisingly, even a single GPS satellite could pro* 
vide useful navigation information to the Space Shuttle on 
orbit . Figure 7.1-5 illustrates a typical situation. During 
periods when the satellite is visible, denoted by the solid 
bar beneath the time scale, position errors are reduced to the 
300-600 ft (lo) level. During periods of no visibility, down- 
range position errors climb to A500 ft. The situation when 
the Shuttle orbit is coplanar with the GPS satellite orbit is a 
singular one, not typical, wherein crossrange errors would not 
be controlled. Velocity errors in the typical case range from 
about 0.5 ft/sec during periods of visibility to 4.5 ft/sec 
when the satellite is out of view. 


I 

t 


5 f 



Figure 7.1-5 Typical Single Satellite Navigation Errors 
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7.2 POST- BLACKOUT NAVIGATION WITH PHASE III GPS 

The post-blackout period during Shuttle reentry is 
one of the most severe for GPS user equipment and may drive 
many of the ultimate specifications for the antenna/preamp/ 
receiver/processor equipments . Even with a 24 satellite. 

Phase III constellation, only three GPS satellites may be vis- 
ible to the Orbiter GPS antennas following blackout. Further, 
it is undesirable to wait 30 or 60 extra seconds while a one 
or two channel sequential receiver demodulates navigation data 
from four satellites. The performance results presented in 
Chapter 4 verified that an atomic time standard will substi- 
tute for one satellite, in terms of visibility, while processor 
storage for navigation data from all 24 satellites will usually 
solve the data gathering time problem . Table 7.2-1 summarizes 
typical post-blackout navigation performance after each of the 
first four GPS measurements (to different satellites) are proc- 
essed by a sequential (single channel) receiver, and parameter- 
izes on the clock type. Beyond the conclusions already noted, 
it was also discovered that the post-blackout velocity errors 
cannot be further reduced below the 1.0 ft/sec level unless 
IMU platform misalignments are modeled in the navigation filter, 


7.3 GPS GROUND TRANSMITTERS FOR SPACE SHUTTLE SUPPORT 

Before the full GPS constellation of satellites is 
available, there will rarely be enough coverage to provide GPS 
navigation during reentry. Even with the full constellation, 
antenna blockage by the main fuel tank will hamper GPS use 
during ascent. Ground GPS transmitters to augment the satel- 
lite constellation would be a way of filling these coverage 
gaps. Several transmitters would be required to support either 
ascent or reentry, use of transmitters might interfere with the 
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TABLE 7.2-1 

POST-BLACKOUT NAVIGATION PERFORMANCE 



RMS VECTOR 
ERROR 

AFTER PROCESSING GPS MEASUREMENTS TO 

CLOCK 

TYPE 

1 

SATELLITE 

2 

SATELLITES 

3 

SATELLITES 

4 

SATELLITES 


Position 

(ft) 

13000 

5000 

1500 

60 

Quartz 

Velocity 

(ft/sec) 

40 

.14 

2 

1.5 


Position 

(ft) 

13000 

■n 

30 

25 

Rubidium 

Velocity 

(ft/sec) 

40 

14 

1.0 

1.0 


satellite signals to the Orbiter or to other GPS users, would 
impose severe design restraints in the Orbiter receiver due to 
the high signal dynamics associated with nearby signal sources. 
Other alternatives (such as varying the Orbiter attitude pro- 
file) may be available to solve the ascent problem, and the 
reentry problem is a temporary one. For all these reasons, 
ground transmitters do not appear desirable for Shuttle support . 


7.4 GPS RECEIVER SUBSYSTEM ANALYSIS 

GPS signal acquisition and reacquisition processes 
may be time consuming for a sequential receiver with only one 
or two channels. The time required to establish a navigation 
fix is referred to as the time-to-first-fix (TTFF), and may 
depend on the users £ priori knowledge of position and time as 
well as on the signal-to-noise density ratio. Figure 7.4-1 
summarizes one important component of TTFF, the time required 
to acquire the C/A code. 
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Figure 7.4-1 Time to Acquire/Reacquire the C/A-Code as a 

Function of C/N , and Position and Velocity 
Uncertainties 


To the C/A acquisition time must be added about 10 
sec to process the C/A signal sufficiently to acquire the pre- 
cise P-code and to track the P-code for a navigation measure- 
ment. An additional 30 sec is required if the navigation data 
message superimposed on the codes must be interpreted. It 
must be remembered that four such measurements are typically 
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required for a navigation fix* but that the receiver may have 
more than one channel working at a time. 

The receiver signal tracking process was addressed by 
considering a specific design, that of the GPSPAC receiver 
currently being designed for space application. That receiver 
utilizes a single IF and VCO common to both the carrier and 
code tracking loops, time sharing the two loops. This design 
minimizes hardware and eliminates vehicle dynamics from the 
loop, but it allows vehicle dynamics to upset the carrier^Toop 
while the code loop is occupying the common IF stage, and it 
results in a 3 dB degradation in the code loop signal* to-noise. 
The dynamic analysis indicates that the loops are, at best, 
marginal for the Orbiter acceleration environment. They may 
be improved by rate-aiding the receiver (by using 1MU outputs), 
by increasing the loop bandwidth, or by going to a third-order 
loop instead of the second-order loop currently envisioned. 
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APPENDIX A 

RELATIVISTIC EFFECTS ON THE SPACE SHUTTLE GPS CLOCK 


This chapter describes the relativistic effects on 
the rubidium time standard being considered for incorporation 
into the Space Shuttle GPS user equipment set. The relativ- 
istic effects are reviewed and then numerically evaluated for 
ascent, reentry blackout, and elliptical orbit flight segments. 
While the relativistic effects are generally small enough to 
ignore for purposes of navigating the Orbiter, under certain 
conditions they are large enough to be the dominant source of 
clock error if uncompensated. At these times, seme adjustment 
should be made in the operation of the GPS-aided navigation 
algorithm. The simplest approach is to properly compensate 
for the relativistic effects at all times in the state propa- 
gation routines. It should be noted that the effects treated 
here are with reference to a ground-based observer at the north 
pcle and are not adequate to compensate the Shuttle receiver 
clock relative to the various GPS satellite clocks. 


A . 1 RELATIVISTIC EFFECTS ON CLOCKS 

An atomic time standard is a "proper clock" to an 
observer that travels with the clock. For that observer, the 
clock runs at a constant frequency within its statistical limits 
of stability. To an outside observer, however, there are two 
relativistic effects that produce deterministic deviations 
from constancy in the clock frequency. One effect is known as 
"time dilatation" or as "second-order Doppler effect" and is 
related to the velocity of the clock with respect to the ob- 
server. A moving clock appears to run more slowly than a 
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clock at rest with respect to the observer. The second effect 
is known as "gravitational redshift" and is related to the 
difference in gravitational potential between the clock and 
the observer. A clock at a lower gravitational potential than 
the observer ("downhill" from the observer) appears to run 
more slowly than a clock at the observer's potential. 

The normalized frequency offset observed for the clock 
due to the above effects is given approximately by (Ref. 31) 


Af _ AU 

r o"? 


V 2 

5 ? 


(A. 1-1) 


In this expression AU is the difference in gravitational po- 
tential between the clock and the observer, v is the magnitude 
of the velocity of the clock relative to the observer, and c 
is the velocity of light (in vacuo). 

One interesting result cf Eq. A. 1-1 is that all clocks 
fixed on the geoid appear to run at the same rate when viewed 
from the north pole. On the geoid, the gravitational poten- 
tial changes with latitude, but so does the velocity imparted 
by the earth's rotation, and the effects on clock frequency 
cancel . 


For purposes of evaluating the relativistic effects 
on a GPS user clock it is useful to scale the normalized fre- 
quency offset into units of velocity (pseudo range-rate) by 
multiplying Eq. A. 1-1 by the speed of light. The integral of 
the scaled, normalized frequency offset is the "accumulated 
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phase offset" to be denoted by <1> and measured in units of dis- 
tance (pseudo range). * Thus: 


do _ cAf _ AU v 2 

cTt ~ f c Tc. 

o 


(A. 1-2) 


The derivative, dp/dt , is referred to as the "phase rate offset. 


Since a GPS user clock is synchronized in phase and 
frequency to GPS system time through the GPS measurements, 
nominal differences between the user clock and system time are 
not important. It is only changes in the clock behavior rela- 
tive to GPS system time (during periods when GPS satellite 
signal availability is inadequate to maintain clock synchro- 
nization) that are important. For the Space Shuttle there are 
three different periods in which significant relativistic ef- 
fects could potentially accrue. 

One such period is the segment from launch to main 
engine cutoff (MECO). Because the external tank will obstruct 
the bottomside GPS antenna, and because the Shuttle flies upside 
down during this period, no satellite signals may be available 
to either GPS antenna. The Orbiter gravitational potential 
and Orbiter velocity change markedly during launch, of course. 

A second such period may occur during reentry when 
radio blackout will deny GPS signals to the antennas. Again, 
Orbiter gravitational potential and Orbiter velocity are chang- 
ing rapidly at this time. 


*Since the speed of light is roughly one foot per nanosecond, 
a one foot phase offset corresponds to a one nanosecond time 
offset , and a phase rate offset of one foot per second cor- 

9 

responds to a frequencv offset of one part in 10 . 
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The third period potentially impacted by relativistic 
effects on the GPS clock is a segment of an elliptic orbit 
with little or no satellite visibility. Since GPS user equip- 
ment is expected to be aboard the Shuttle in mid-1981 when 
only six GPS satellites will be in place, there will be sig- 
nificant visibility gaps on orbit. A complete visibility gap 
of one-half an orbit (about 45 minutes) will be possible. For 
a circular orbit, both the gravitational potential and the 
Orbiter velocity are constant, so there is no relativistic 
fluctuation in the clod frequency. Elliptic orbits do, how- 
ever, produce relativistic fluctuations. 


A. 2 SPACE SHUTTLE ASCENT AND REENTRY BLACKOUT 

To evaluate rhe relativistic effects on the clock 
during ascent and reentry blackout, it is necessary to have, 
in addition to Eq. A. 1-2, an expression for the gravitational 
potential in the vicinity of the earth and relevant trajectory 
characteristics of the orbiter during these periods. The gravi- 
tational potential is given approximately by 

U = - £ , (A. 2-1) 

where m is the earth's gravitational constant, and r is the 
radial distance from the center of the earth to the point in 
question . 


Calculation shows that the gravitational redshift of 
the clock is of little importance during ascent and reentry 
blackout. Table A. 2-1 summarizes these results. 

An OFT-1 simulated trajectory (Ref. 32) was used for 
the ascent segment. Total elapsed time from launch to MECO 




THE ANALYTIC SCIENCES CORPORATION 


TABLE A. 2-1 

GRAVITATIONAL REDSHIFT DURING ASCENT AND REENTRY BLACKOUT 


PERIOD 

EVENT 

RADIAL 
DISTANCE (r) 
(ft) 

PHASE RATE 
OFFSET (d<j>/dt) 
(ft/sec ) 


LAUNCH 

2. 09099* 10 7 

0.0000* 

ASCENT 

MECO 

2 . 12903*10 7 

0.0122 

REENTRY 

ENTER 

2. 12097*10 7 

0.0000* 

BLACKOUT 

EXIT 

2.10480*10 7 

-0.0052 


^Assumed point of clock calibration. 

for this trajectory is 527 sec. The peak redshift phase rate 
offset is 0.0125 ft/sec, occurring at about 370 sec after launch, 
when the trajectory is at its maximum altitude for this segment. 
Altitude decreases somewhat before MECO, leading to the smaller 
offset value reported in Table A. 2-1. The accumulated phase 
error from launch to MECO due to the redshift is less than 
5 ft . 


The standard deorbit trajectory first discussed in 
Section 2.2.2 was also used for the reentry blackout calcula- 
tions. It is not certain just when the blackout will occur, 
so a conservative assumption was made that blackout will begin 
at 300 kft of altitude and end at 145 kft of altitude. A pe- 
riod of 890 sec elapses between these two altitudes. The peak 
redshift occurs at exit from blackout and is noted in Table 
a. 2-1. The accumulated phase offset during blackout due to 
the redshift would be less than 3 ft. 

The second-order Doppler shift effects during ascent 
and reentry blackout are not so small. Figure A. 2-1 plots Lhe 
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Time imc) 


Figure A. 2-1 Velocity versus Time for Ascent and 

Reentry Blackout 

inertial velocity magnitude during these two flight segments. 
Because exit from blackout is uncertain, various terminal alti 
tudes are noted on this and subsequent figures. (The initial 
altitude for blackout is still taken as 300 kft). Figure A. 2- 
shows the resulting clock frequency offset due to the second- 
order Doppler shift, and Fig. A. 2-3 plots the integrated clock 
phase offset. For both these latter figures the sign of the 
offsets during ascent is negative. 

For purposes of comparison, the statistical phase 
error is plotted in Fig. A. 2-4 for a typical rubidium clock 
calibrated to GPS system time aboard the Orbiter. Only one 
hour of elapsed time is shown in Fig. A. 2-4, as longer time 


ATiOUOCX Pham OffiCI #1) CiOCy PMAtt AATC Offt€T tH/mct 


THE ANALYTIC SCIENCES CORPORATION 



Figure A. 2 
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Figure A. 2-4 Statistical Performance of a Rubidium 

Oscillator Calibrated to GPS Time 
Aboard the Orbiter 

periods are generally not of interest for Orbiter navigation 
purposes. Over this short time period the statistical stabil- 
ity of the clock is much better than Fig. A. 2-4 would indicate. 
The (la) clock errors shown are almost totally the result of 
errors in the initial calibration of the clock to GPS system 
time. Excellent satellite visibility and geometry, use of 
both the Lj and L 0 satellite signals for ionospheric delay 
compensation, and use of Phase III GPS satellites should allow 
a user clock calibration on orbit to about 4 ft in phase and 
about 0.0046 ft/sec in frequency (lo). These values, along 
with a five state stochastic model of typical rubidium oscil- 
lator performance, were used to generate Fig. A. 2-4. 

Clearly the relativistic clock effects during ascent 
and reentry blackout are larger than the expected errors due 
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to clock instability and calibration error. The subject of 
what to do about the relativistic effects is postponed until 
the effects on an elliptic orbit are examined. 


A. 3 RELATIVISTIC EFFECTS ASSOCIATED WITH ELLIPTIC ORBITS 


The gravitational redshift is not important compared 
tc the second-order Doppler shift during ascent and reentry 
blackout. On orbit, however, the two effects are of equal 
importance . 


The square of the total velocity of a spacecraft in 
an elliptical earth orbit is given approximately by 


v 



(A. 3-1) 


where "a" is the semimajor axis of the orbit. Since the clock 
is to be calibrated in orbit, it is only variations in the 
square of the velocity that produce second-order Doppler shift 
effects. These variations are inversely proportional to the 
radial distance of the orbiter, as Eq. A. 3-1 states. The 
variations in the gravitational redshift are also inversely 
proportional to the radial distance, as Eqs. A. 1-2 and A. 2-1 
jointly state. Inserting Eqs. A. 2-1 and A. 3-1 into expression 
A. 1-2 for the total phase rate offset, it is apparent that the 
total effect is the sum of equal amounts of redshift and recond- 
order Doppler shift, and the resultant effect is given by 


do _ 2 m / 1 1 v 

dt c v r Q " r(t)' 


(A. 3-2) 


In this expression, r Q is the radial distance to the spacecraft 
at the time of clock calibration, and the explicit dependence 
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of the radial distance to the spacecraft upon time has been 
noted . 


It would be convenient to have a closed form expres- 
sion for the dependence of the phase rate offset upon time. 
Equation A. 3-2 does not provide such an expression because 
there is no closed form expression for the radial distance to 
a body in an elliptical orbit as a function of time. Equation 
A. 3-2 is the "next best" expression, involving only the radial 
distance. A similarly useful expression is derived for the 
phase offset in the remainder of this section. 

Clearly, the phase offset is the integral of Eq. A. 3-2 
t 


dt 

r( t ) 


(A. 3-3) 



which may be rewritten as 


_ 2jj 1 t 

f r(t) dr 1 

c l r 
( ° 

J rvdr./dt) i 

O ' 


(A. 3-4) 


For a body in an elliptical orbit of semimajor axis 
"a" and eccentricity "e" . 


dr _ 
dt " 


4 



at 1 -e 

y 

r 



(A. 3-5) 


which expresses dr 'dt in terms of r for use in Eq. A. 3-4. The 
upper t + ' sign holds for that portion of the orbit from perigee 
to apogee, while the lower t-> sign holds for the remaining 
half. 
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The indicated integration in Eq. A. 3-4 may be performed 
analytically by use of Eq. A. 3-5. The result is 


0>( t) 




(A. 3-6) 


This result is valid only over segments of an orbit wholly 
contained between apogee and perigee or vice versa. 

To see how large 0(t) may grow in one half of an orbit, 
consider that the spacecraft is at apogee at time zero. That 
is 


r Q = a(l + e) 


Then 


®(t) 


2jj 1 t 
c | a(l+e) 



(A. 3-7) 


(A. 3-8) 


One half of an orbital period, T/2, later, the spacecraft is 
at perigee where 


r(T/2 ) = a(l-e) 


(A. 3-9) 


and 

T = 2na\P - 


(A. 3-10) 


Thus . at perigee 


<t> (T/2 ) 


2 ny au e 

c 1+e 


(A. 3-11 ) 
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af(T/2) 



(A. 3-12) 


For a numerical example, consider a deorbit trajectory 
with an apogee at 400 run altitude and a perigee at the earth's 
surface. For this orbit (an extreme elliptical orbit for the 
Space Shuttle) 

a = 2.214086 x 10 7 ft, 
e = 5.488599 x 10~~ , 

T = 5517 sec (92 min). 


and 

♦(T/2) = -185.6 ft, 

^(T/2) = -0.1423 ft/sec. 

Again, comparison with Fig. A. 2-4 shows that the rela- 
tivistic effects are larger by far than the likely clock errors 
due to instability. 


A. 4 CLOCK COMPENSATION* FOR RELATIVISTIC EFFECTS 

Comparing the relativistic clock effects during ascent, 
reentry blackout, or elliptical orbit with the statistical 
clock performance figures, it is clear that the inherent sta- 
bility of the clock will not be achieved without compensation 
for the relativistic effects. For purposes of Orbiter naviga- 
tion, however, the relativistic clock offsets should be toler- 
able without compensation, even though limited satellite visi- 
bility and poor geometry can multiply clock errors by a factor 
of 10 or more in producing navigation position and velocity 
errors (Chapter 4). 
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If, as is a virtual certainty, a Kalman (or Kalman- 
like) filter is used to process GPS measurements, then some 
accounting for the relativistic effects will be necessary for 
proper filter operation . The clock phase and frequency vari- 
ances of the navigation filter will have to be sized to accom- 
modate any uncompensated relativistic effects. To do otherwise 
would be to produce unrealistically small clock variance esti- 
mates at certain times, leading to improper filter operation 
and probably to the rejection of good measurement data at cri- 
tical times. It is simpler to incorporate deterministic com- 
pensation for the relativistic effects than to (improperly) 
adjust the filter clock variances on special occasions. 


The equations for compensating the clock for the rela- 
tivistic effects are easily derived. Using <? to represent the 
clock phase (measured in units of distance) as before, and 
using f here to represent the clock phase rate (measured in 
units of velocity) then, clearly 


dt 


f 


(A. 4-1) 


df _ p dr v dv 

dt ' 2 dt ‘ c dt 

cr 


(A. 4-2) 


In discrete form 

f, = f + Ar - ^ Av , (A. 4-3) 

1 ° c 
cr 

<1>1 = i> 0 + T ( f 0 + f^) At, (A. 4-4) 

where o Q , f Q are the phase and phase rate before the time step 
At, and f^ are the phase and phase rate after the time 

step. The new terms Ar and Av refer to changes in the radial 
distance and in the total inertial velocity during the interval. 
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The addition of these simple compensation equations 
to the GPS R/PA navigation state propagation algorithm will 
reduce relativistic effects to a negligible level and will 
preserve the inherent capability of the onboard atomic time 
standard. 
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APPENDIX B 

GPS DILUTION OF PRECISION 


This chapter is a tutorial introduction to the subject of 
measures of dilution of precision (DOP) for the NAVFTAR Global 
Positioning System. Geometric (GDOP) , position (PDOP), time 
(TDOP), horizontal (HDOP), and vertical (VDOP) dilutions of 
precision are discussed. After defining the DOP's for the 
normal case of measurements to four GPS satellites, the case 
of more than four measurements is treated. Coordinate invari- 
ance of GDOP, PDOP, and TDOP is demonstrated, and the singular 
situation when all satellites lie on a cone with apex at the 
user is introduced. To obtain a quantitative feel for the 
DOP's, a particular one-parameter geometric situation is exam- 
ined in detail, and minimum values for the DOP's are established. 

Section B.l summarizes the motivation for performing 
DOP analysis. Section B.2 defines several measures of dilution 
of precision and also reviews the manner in which the GPS re- 
ceiver generates and processes ranging measurements. Measure- 
ments to more than four satellites may be made, and Section 
B.3 addresses this issue. Some measures of dilution of pre- 
cision do not depend on the coordinate systems used to compute 
them, a fact demonstrated in Section B.4. Certain geometrical 
situations are singular - they do not yield enough information 
to allow the user to solve for his position -- and Section B.5 
identifies such situations. Section B.6 is devoted largely to 
an examination of one family of geometrical situations, a fam- 
ily with interesting properties for minimization of dilution 
of precision. 
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B.l INTRODUCTION 

The NAVSTAR Global Positioning System (GPS) is a sat- 
ellite-based radio ranging navigation system scheduled to be 
operational in the mid-1980's. The transmitters are located 
in satellites in circular orbits of about four earth radii, 
having orbital periods of 12 hr and ground tracks that repeat 
once per day. Twenty- four satellites are planned for the oper- 
ational GPS constellation, with fewer (five or six) to be ini- 
tially available. Ranging signal transmission is in the L 
band, so signal reception is limited to "line-of-sight" visi- 
bility between user and satellite. (Clouds and weather wil 1 
not block the signals, but the earth, the user's own vehicle, 
and other obstructions may.) 

GPS users will need to measure one-way range (offset 
by user clock errors) to four or more satellites to determine 
their position. Since generally more than four satellites 
will be visible, a choice of satellites usually must be made. 
This choice, and the user position and time errors that result 
from it, depends upon the geometry of the satellites with re- 
spect to the user. The effect of the geometry is generally to 
dilute the precision available in any single ranging measure- 
ment, producing amplified user position and time errors. Errors 
along particular spatial orientations may be reduced instead 
of amplified, but the effect is still known as dilution of 
precision ( DOP ) . 

B . 2 DEFINITION OF DILUTION OF PRECISION 

A GPS receiver measures what are termed "pseudo ranges" 
to a number of GPS satellites (usually four, but possibly more) 
and then solves for the receiver position and time. Errors in 
the pseudo range measurements cause errors in the determined 
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position and time that are functions of the receiver- to-satel- 
lite geometry. Geometric dilution of precision (GDOP) and its 
relatives are scalar measures of the effect of measurement 
geometry in translating the pseudo range errors into receiver 
position and time errors. The dilution of precision values 
are useful because they reduce the effect of a complex geo- 
metric situation into a single number. The various OOP's have 
nothing to do, however, with the pseudo range measurement er- 
rors themselves, and they accurately describe the effect of 
measurement geometry only when certain restrictive assumptions 
concerning the statistics of the pseudo range errors are satis- 
fied. This chapter is a tutorial attempt to explain the OOP's 
and demonstrate some of their properties. 

A good place to begin this discussion is with an ex- 
planation of the measurement process itself. Figure B.2-1 is 
a conceptual diagram of the pseudo range measurement process. 
Each GPS satellite carries a clock that drives a signal gen- 
erator producing a unique, pre-determined function of time. 

The receiver also incorporates a clock that drives a signal 
generator capable of mimicking any of the satellite signal 
generators. A key conceptual difference between the receiver 
and the satellite is that the receiver is capable of retarding 
the signal generator output by a variable time difference, p. 
The retarded, self-generated signal is compared in the receiver 
to the received satellite signal, and the time di f ference is 
continuously adjusted to keep the two signals in synchronism. 

At any time when the two signals are sufficiently in synch, a 
measurement of pseudo range may be made by simultaneously not- 
ing both the time difference, p, and the receiver time, t f . 

The GPS pseudo range measurement is essentially a 
measure of the time difference between transmission (from the 
satellite) and receipt (by the user) of the GPS signal. The 


B-3 




THE ANALYTIC SCIENCES CORPORATION 


1 


satellite 





TRANSMISSION. 


r 


L 



~i 




Figure B.2-1 Conceptual Diagram of GPS Pseudo-Range 

Measurement Process 

measured time difference is the time it took the signal to 
travel from the satellite to the receiver, plus any offset in 
the receiver clock with respect to the satellite clock (plus 
any uncompensated signal processing delays in either the trans- 
mitter or the receiver). Assuming that the signal travels at 
the constant speed of light (making various corrections for 
ionospheric and tropospheric effects) and that the speed of 
light has been used to scale time intervals so they are mea- 
sured in units of distance , then 


*The speed of light is nearly one foot per nanosecond, so a 
time interval of one nanosecond is measured as a pseudo- 
range of about one foot. 
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p i 3 V(x r - Xi ) 2 + (y r -y.) 2 + u r -z i> 2 + 6t ri . (B.2-1) 

t* h 

In this equation, p. is the pseudo range to the i satellite; 

A 

(x.. y^ , 2 .) are the position coordinates of the i satellite 
at the time, t r -p., the signal was "sent"; (x^, y , z r > are 
the position coordinates of the receiver at the time, t^, the 
signal was "received"); and 6t . is the offset of the receiver 

r l ri 

clock with respect to the i satellite clock. The radical 
term is, of course, the true range from the satellite to the 
receiver. 


Information transmitted along with the satellite sig- 
nals allows the receiver to estimate both the position coordi- 
nates of any satellite at any time and the offset of any sat- 
ellite clock with respect to GPS system time. The satellite 
clock offsets may be used to adjust the measured pseudo ranges 
so that they reflect the offset of the receiver clock from GPS 
system time, effectively removing the subscript i from 6t r ^ in 
Eq. B.2-1. 

Equation B.2-1 thus contains four unknowns related only 
to the receiver position and clock offset: (x r , y , z r , 6t r ). 

Simultaneous measurements of pseudo range to four satellites 
will then allow the receiver to solve for its position and 
clock offset, provided that the measurement geometry is such 
that the four measurements supply four independent pieces of 
information . 

The question remains: How are errors in the measured 

pseudo ranges related to errors in determined receiver position 
and clock offset? Differentiation of Eq. B.2-1 leads to an 
approximate linear relationship for th • errors 
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-Ap . 5 u. Ax + u. Ay +■ u. A 2 + At, (B.2-2) 

K i ix r ly * r iz r r’ 

where Ap^ is the error in the i pseudo range, (Ax r , Ay r , 

Az r > are the coordinates of the resultant error in receiver 

position, and At r is the resultant error in receiver clock 

offset. The new symbols, (u. , u. , u. ). of Eq. B.2-2 are 

ix ly 12 

the components of the unit vector from the receiver toward the 
i th satellite 


A 

U . 
1 



where 


R. 

l 




(x.-x ) 2 + (y.-y) 2 + (z.-z 


Clearly 






1 . 


(B.2-3) 


<B. 2-4 ) 


(B.2-5) 


The four necessary equations of the form of Eq. B.2-2, 
corresponding to pseudo range measurements to four satellites, 
may be combined into one matrix equation 


where 


A 


Ax = Ap , 


Ax = 



(B.2-6) 


(B.2-7) 
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Ap 



-Ap, 


-ap 3 

\-ap 4 


(B. 2-8) 


and 


r lx 

u ly 

u lz 

| u 2x 

u 2y 

u 2z 

1 u 3x 

u 3y 

u 3z 

Wx 

u 4v 

u 4z 



(B.2-9) 


iXote the inclusion of the minus signs in the definition of 
the vector of measurement errors, Ap , for notational conveni- 
ence . ) 


The statement that the four pseudo range measurements 
are independent is equivalent to the statement that the matrix 
A is invertible. Assuming such to be the case, the receiver 
errors are given by 


Ax = A -1 Ap . (B.2-10) 

The discussion now becomes statistical. Assuming 
that ihe expected (or mean) values of the pseudo range mea- 
surements are all zero , then the same is true of all elements 
of the receiver error vector, Ax. Furthermore, the covariance 
of the receiver errors is given bv 

P = E{ Ax Ax T } = A * 1 E{Ap Ap 1 } A~ T , (B.2-11) 

_T 

where E is the ensemble expectation operator, and A is the 
transpose of the inverse of A. Assuming that the pseudo range 
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errors are independent of each other and have a common ariance , 
o~ , or 


E{Ap Ap T } = o 2 I (B.2-12) 

(with 1 the 4x4 identity matrix), then the receiver error co- 
variance matrix is given by 

P = o 2 (A T A ) _1 . (B.2-13) 

Assuming, further, that the pseudo range errors are normally 
distributed . then so are the receiver position and time errors, 
and Eq. B.2-13 provides a complete statistical description of 
the resultant receiver errors due to pseudo range errors. 

To focus on the effect of the measurement geometry 

instead of on the pseudo range errors, the receiver error co- 

variance matrix will be normalized by replacing the pseudo 

2 

range variance, o , in Eq. B.2-13 by unity. Thus, 

P = (A T A) _1 . (B.2-14) 

To reduce further to a single number, define the geometric 
dilution of precision as 

GDOP = % 'tr(f>) , (B.2-15) 

where the symbol tr stands for the trace of its matrix argu- 
ment. that is, for the sum of the diagonal entries. GDOP is 
thus the root sum square (rss) of the errors in receiver posi- 
tion (on all three axes) and time offset, normalized to unit 

variance pseudo range errors. 
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GDOP is a useful descriptor of the geometric situa- 
tion when all four receiver errors are of importance. Often 
this is not the case, and related DOP's are more useful. When 
only the three position errors are important, the useful de- 
scriptor is position dilution of precision 

PDOP = V p u + P 22 + P 33 , (B.2-16) 

which is the (normalized) rss radial position error. If only 
the time is of importance, then time dilution of precision is 
given by 


TDOP = VP 44 , (B.2-17) 

and is simply the (normalized) one-sigma clock offset error. 

To go further requires particularization of the cal- 
culations to a local coordinate system at the receiver. If 
the coordinate system is chosen so that the z-axis is along 
the local vertical (putting the x- and y-axes in the horizon- 
tal plane), then horizontal and vertical dilutions of preci- 
sion are defined as 


HDOP = >/P n + p 2 2 * (B.2-18) 

VDOP = VP^ . (B. 2-19) 

B . 3 MORE THAN FOUR MEASUREMENTS 

The discussion to this point has proceeded as if exact- 
ly four independent pseudo range measurements were always made. 
With fewer than four measurements, Eqs. B.2-2 cannot be solved 
due to insufficient data. Conversely, with more than four 
measurements, there is no exact solution to the overdetermined 
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set of these equations. With any number of measurements, the 
equations may be put in the form of Eq. B.2-6, by altering 
Eqs. B.2-S and B.2-9 so that the number of rows of A and Ap 
equals the number of measurements. With other than four meas- 
urements, the matrix A is not square. With fewer than four 
measurements, additional information of some sort is necessary 

With more than four measurements, the standard procedure is to 

T 

multiply Eq. B.2-6 by A to obtain 

A T A Ax = A T Ap . (B.3-1) 

T 

The matrix A A is a 4x4, and, if it is invertible (if the more 
than four measurements contain at least four independent ones) 
then Eq. B.3-1 may be "solved" to yield 

Ax = (A T A) -1 A T Ap . (B. 3-2) 

This solution for the receiver error vector will, in general, 
not satisfy any of the original equations, but it represents 
the well-known "least squares" solution to the overdetermined 
set of equations, minimizing the sum of the squares of the 
amounts by which the equations are not satisfied. This least 
squares method of solving for the errors in determined posi- 
tion due to errors in pseudo-range measurements is only ap- 
propriate if the position estimate itself was obtained from 
the pseudo range measurements by the least squares process. 

Making the same statistical assumptions about the 
(now more than four) pseudo range errors as were previously 
made yields the same formal results; namely, that Eq. B.2-13 
and its normalized form Eq. B.2-14 are still valid, and the 
definitions of the DOP * s are unchanged. Thus, the dilution of 
precision concept is useful in describing the situation involv 
ing more than four measurements, with the basic definitions 
virtually unchanged. 
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B.4 COORDINATE INVARIANCE OF GDOP, PDOP , AND TDOP 

As was remarked earlier, in order to compute the hor- 
izontal or vertical dilutions of precision an appropriate co- 
ordinate system must be used. It might be expected that GDOP, 
PDOP, and TDOP do not depend upon which coordinate system is 
chosen for the computations. This coordinate invariance is 
demonstrated here. 

Assume that the computations have already been per- 
formed in one coordinate system with the notation previously 
given. Assume a second coordinate system related to the first 
by the orthonormal 3*3 rotation matrix M, so that the unit 
vectors from the receiver toward the satellites are given by 
in the new coordinates, where 

v. = M u. , (B.4-1) 

1 1 ’ 

and M satisfies 

M" 1 = M T . (B.4-2) 

What was termed the A matrix in the original coordinates will 
become the B matrix in the new coordinates. Clearly 




v iy 

v lz 

x \ 

B = 

i 

v 2y 

• 

• 

v 2z 

1 1 

. , (B.4-3) 


iv'- 

V ’ 

v* ' 

i 


' nx 

ny 

nz 

9 


where the general case of measurements to n satellites is con- 
sidered. Reflection shows that 

B = AMj , 


B-ll 


(B.4-4) 
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where the new, 4*4 matrix is given in partitioned form by 



: 0 1 

I . 

(B.4-5) 

\ 0 

» 1 i 

1 



Like M, is orthonormal and satisfies 


M -1 _ m T 
\ ’ M 4 * 


(B.4-6) 


The normalized receiver error covariance matrix in 
the new coordinates becomes Q, where 


Q = ( B T B ) * 1 


(B.4-7 ) 


Clearly , 


Q = (M^AMj )" 1 , 

= p m l • 


(B . 4-8 ) 


If P is partitioned as 

P 


(•!, . : ) . 

\ P I P44 / 


(B.4-9) 


then Q is given as 



Mp \ 

U T M T j 

P 44 / 


(B.4-10) 


As the 4,4 element of Q is identical to P,, , it is 
clear that TDOP is the same when computed in the second coor- 
dinate system as it is when computed in the first. As for 
PDOP . in the first system 
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PDOP 2 = tr(P 3 ) , 


(B.4-11) 


while in the second. 


PDOP 2 = tr(M P 3 M T ) . (B.4-12) 

The trace function has the general property that 

tr(CD) = tr(DC) (B.4-13) 

for any matrices C and D whenever their dimensions give sense 
to both products, (l.e., if C is n*m, then D must be m*n. ) 

Thus 

tr(M P 3 M T ) = tr(P 3 M T M) = tr(P 3 > , (B.4-14) 

and PDOP is the same in both coordinate systems. 

The identical argument applied to the full 4x4 P ma- 
trix shows that GDOP is the same when computed in either system. 
Or, the fact that 


GDOP 2 = PDOP 2 + TDOP 2 , (B.4-15) 

combined with the invariance of PDOP and TDOP will demonstrate 
the coordinate invariance of GDOP. 


B . 5 SINGULAR MEASUREMENT GEOMETRY 

It has been remarked that for certain geometrical 
conditions the pseudo range measurements might not yield four 
independent pieces of information. In this singular case, it 
is not possible to solve Eqs. B.2-1 for the receiver position 
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and clock offset. Also, Eq. B.2-6 will not be solvable for 
the receiver error vector due to errors in the pseudo range 
measurements. The A matrix will be of rank three or less, and 
A A will be singular and have no inverse. In the case of exact- 
ly four measurements, A itself will be singular. The question 
is: For what geometries does this singular condition occur? 

The often heard answer is that if all the satellites 
lie on any (circular) cone with apex at the receiver then the 
geometry is singular. This fact, and its converse, will be 
demonstrated here. 

It is easy to show that if the satellites all lie on 
a cone with apex at the receiver, then the A matrix is of rank 
three (or less), and A A is singular. By the results of the 
previous section, it is clear that it makes no difference what 
coordinate system is used to express the A matrix; if it is of 
rank three in one coordinate system, then it is of rank three 
in all coordinate systems. Supposing all the satellites to 
lie on a cone with apex at the receiver, select the axis of 
the cone as the z-axis of the coordinate system. Then the 
unit vectors from the receiver toward the satellites will all 
have identical z components equal to the cosine of the half- 
cone angle, a. This geometry is illustrated in Fig. B.5-1. 

The A matrix now looks like 

Ujy cos a 1 

U 2 y cos a 1 

• *4 

• *4 

4 • 4 

u„„ cos a 1 

ny 

Clearlv the last two columns of A are scalar multiples of one 

T 

another, so A can be of rank three at most, and A A must be 
singular . 
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Figure B.5-1 


Singular Geometry -- Satellites Lie on a Cone 
With Apex at the User 


The converse statement, that if A is of rank three 
(or less), then all the satellites lie on a cone with apex at 
the receiver, is not much harder to show. If A is of rank 
three, only three of its rows may be linearly independent. 
(The case when fewer than three rows are linearly independent 
will be an obvious extension. ) As it makes no difference how 
the rows are numbered, move these three rows to the top. The 
remaining dependent row (or rows) may then be expressed as 
linear combinations of the top three rows. That is, for k = 
4,...,n there are constants ^k2’ an< * ^k3 suc ^ 
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<u kx u ky u kz 11 * A kl (u lx u ly “lz 1J 

* *k2 (u 2x u 2y u 2z 1J 

* *k3 (u 3x u 3y u 3z 1! • (B.5-2) 

Examination of the fourth component of this equation reveals 
that 


X kl + X k2 * X k3 = 1 • 


'B.5-3) 


Any three vectors must lie on a cone with apex at the 
origin. (This is equivalent to the statement that any three 
points on the unit sphere must lie on some small circle.) So 
the unit vectors to the satellites involved in the top three 
rows of the A matrix must lie on some cone. Take the axis of 
this cone as the z-axis of. the coordinate system, then the z 
components of the first three rows must equal the cosine of 
the half cone angle, or 


u lz 5 u 2z = u 3z = cos “ • (B.5-4) 

Equations B.5-2 and B.5-3 now show that the z compo- 
nents of all the rows of A must be equal, 


x kl u lz * *k2 u 2z + A k3 u 3z = u kz = cos 0 (B ' 5 - 5 

hence all the satellites must lie on a cone with apex at the 
receiver. 


When the satellites do all lie on a cone, then the 
various DOP's do not exist. When the satellites are nearly on 
such a cone, then the DOP's exist but will, in general, be 
large. It might happen that VDOP, HDOP, or TDOP will be of 
reasonable size in this case, but GDOP and PDOP will be large. 


B-16 




THE ANALYTIC SCIENCES CORPORATION 


B . 6 'N IMPORTANT EXAMPLE SITUATION 

The reason for defining the various DOP's is to re- 
duce complex geometrical situations to a single number, or to 
a few numbers, with direct relevance to GPS navigation and 
time accuracy. Equations B.2-15 through B.2-19 defining the 
DOP’s accomplish that objective, but they give no immediate 
insight into the relationship between the DOP's and the 
geometry. 


One way to use the DOP's is to calculate them for 
every alternative geometry in a given situation and choose the 
satellites for tracking which yield the best DOP values. For 
example, if nine satellites are visible to a receiver but only 
four are to be used, the DOP's may be calculated for each of 
the 126 ways four satellites could be selected from the nine 
visible, and the set of four to be used could be taken as the 
one with the appropriate minimum DOP. No geometrical insight 
is needed for this exhaustive process. 

To decide the acceptability of a single DOP in a non- 
exhaustive situation, or to otherwise shorten the computation- 
ally burdensome process of exhaustively computing DOP's, it is 
useful to have a quantitative feel for the relationship between 
geometrical situations and the DOP's, and for the potential 
ranges of the DOP's. 

It would be ideal if there were some simply under- 
standable. closed- form expression for the DOP's in terms of 
obvious parameters describing the geometrical situation, with- 
out reference to unnecessary, complicating details such as the 
orientation of the coordinate system in which the DOP is com- 
puted. In the simplest case for which the DOP’s exist, that 
of measurements to four satellites, there are eight independent 
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parameters necessary to describe the orientations of the unit 
vectors toward the satellites. If one is considering GDOP, 

PDOP , or TDOP, where the coordinates of computation are not 
relevant, then the number of independent parameters may be 
reduced to five. Describing a single unit vector by two inde- 
pendent parameters, however, is not simple; it involves either 
trigonometric functions of two angles or square roots of direc- 
tion cosine parameters. After inserting several of these ex- 
pressions into the matrix multiplication and inversion proc- 
esses used to compute DOP's, the results would be too complex 
to be easily understandable. So the process of gaining a quan- 
titative feel of the relationship of the DOP's to the geometry 
will begin with the close analysis of a constrained geometrical 
situation having only one free parameter. 

This situation is most easily described with refer- 
ence to a locally level coordinate system. Let one satellite, 
the first one, be directly overhead (along the z-axis). The 
other three satellites are separated from one another in azi- 
muth by 120 deg, and all have the same, unspecified coeleva- 
tion (which may vary from zero to 180 deg). The common co- 
elevation angle, or, of the three satellites is the only vari- 
able parameter of this situation. For completeness, assume 
that the vector from the receiver to the second satellite has 
an azimuth of zero, the third satellite has an azimuth of +120 
deg, and the fourth satellite has an azimuth of -120 deg. 

Figure B.6-1 illustrates this geometry. 

Of course, the geometry need not be referenced to any 
local level frame to consider GDOP, PDOP, or TDOP. It may 
simply by described by saying that the vectors toward three of 
the satellites are evenly spread around a cone of half-angle 
a. while the vector toward the fourth satellite lies on the 
axis of the cone (pointing either into or out of the cone). 
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Figure B.6-1 Single Parameter (Constrained) 

Geometry Situation 


The unit vectors toward the four satellites are given 


u x = ° , u 2 = 


/ sin a 
0 

\ cos a 


(1/2) sin a\ 


u^ s | (-v'3/2)sin a 


cos a 


-(l/2)sin or 
0 4 s | -(V3/2)sin or 
cos a 


(B.6-1) 
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The A matrix of Eq. B.2-9 becomes 


A = 


0 

sin a 

-(l/2)sin a 
- ( 1/2 )sin a 


0 

0 

W3/2)sin a 
-(V3/2)sin a 


whereupon 


1 

cos a 
cos or 
cos or 



» 


(B. 6-2) 



[(3/2) sin 2 a 

0 


0 

° \ 

T 

0 

(3/2) sin 2 a 


0 

0 

a t a = 







i ° 

0 

1 + 3 

2 

cos a 

1+3 cos a 


\ 0 

0 

1 + 3 

cos a 

4 / 


(B. 6-3) 


It follows that the normalized receiver error covariance matrix, 

T 

the inverse of A A, is given by 


o • 2 

3 sin a 


P = 


3 sin a 
0 


3(l-cos a)' 


1 + 3 cos a 
3(l-cos a)^ 


Q+3 cos a) 
2 

3( 1 -cos a) 


1+3 cos a 
3(1-cos a) A 


(B . o-4) 


The various DOP ’ s are given by 
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®° p - dhr Vi • (B - 6 ‘ 5) 

VD0P * (1-cos a T Vi • (B.6-6) 

PD0P ' (1-cos T) Vl(l*cos a) • 

TD0P * (1-cos a) . (B.6-8) 


GDOP 


1 _ 

(1-cos 


J ~ 

®Tl( 


cos a 


♦3 cos a + 3 
3(l+cos a ) 


cos'* a 


(B.6-9) 

An issue of immediate interest is the situation in 
which these DOP's are minimized and the minimum values thus 
obtained. It is readily apparent that HDOP is minimized when 
a is 90 deg, putting three of the satellites in the horizontal 
plane . The values of the DOP’s for this situation are 

HDOP = 2 : 1.1547, 

VDOP = 2 ^ : 1.1547, 

HDOP = 2 : 1.6330, 

TDOP = VI : °- 5774 * 

GDOP = S3 : 1.7321. 

The situation in this constrained case which "minimizes 
VDOP would have a be 180 deg, putting three satellites directly 
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below the receiver and one above. This situation is a singu- 
lar one, however, causing HDOP, PDOP, and GDOP to blow up. 

A single situation simultaneously minimizes PDOP, 

TDOP . and GDOP. It has a at about 109 . 47 deg (so that cos or = 
-1/3), meaning that the three satellites are all about 20 deg 
below the horizon. In this situation the four unit vectors 
toward the satellites all terminate at the vertices of a regu- 
lar tetrahedron , illustrated in Fig. B.6-2. This situation is 
completely symmetric in that each unit vector makes the same 
angle, 109.47 deg, with each of the other unit vectors. While 
the minimization taking place in this example is only that of 
selecting the single parameter in a constrained, one-parameter 
geometrical situation, the resultant minimum situation illus- 
trated in Fig. B.6-2 happens tj be a global minimum for PDOP, 
TDOP , and GDOP . That is. no better values for these DOP’s can 
be obtained by any other arrangement of four satellites, a 
fact stated here without proof. The values of the DOP’s in 
this situation are 

HDOP = : 1.2247 

VDOP = : 0.8660 

PDOP =4 = 1.5000 \ 

. I Absolute 

TDOP - j = 0 . 5000 \ minimum for 

I 4 satellites 

GDOP = * 1.5811 ) 

Of course, this global minimum situation would be 
impossible to obtain near the surface of the earth, for tht 
earth itself would block visibility to the three satellites 
below the horizontal plane. A better picture of how the DOP’s 
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Figure B.6-2 Global Minimum for Four-Satellite GDOP, 

PDOP , and TDOP 

vary with the coelevation angle is given in Fig. B.6-3 (a and 
b), which plots the DOP’s versus the coelevation angle. 

Leaving, for a paragraph, the one-parameter geometry 
and taking up the topic of global minimums for the DOP’s. a 
global minimum for HDOP is approached as a singular limit of a 
situation where the unit vectors to the four satellites approach 
the plus and minus x and y axes in the horizontal plane. The 
limit of HDOP in this situation is 1.0, but the limiting situa- 
rion is singular, so HDOP does not exist. Similarly, a global 
minimum for YDOP would be the limit of a singular situation 


' ; *The DOP's could be redefined using matrix pseudo- inverses in- 
stead of ordinary inverses. Then, in some singular situations, 
certain of the DOP’s would still be well defined, as would 
HDOP in the situation indicated. 
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Figure B.6-3(a) Variation of VDOP, HDOP. and PDOP with 

Coelevation Angle for a One-Parameter 
Geometrical Situation 





Figure B. 6- 3(b) Variation of TDOP, PDOP. and GDOP with 

Coelevation Angle for a One-Parameter 
Geometrical Situation 
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with two satellites along the positive vertical (z) axis and 
two along the negative vertical axis. The limit of VDOP in 
this situation is 0.5. 

Returning again to the one-parameter geometry, Fig. 
B.6-3(b) has utility beyond its example purposes. It frequently 
happens that satellite visibility to a user is restricted to 
some cone with apex at the user. This could be because the 
earth masks satellites, leaving the user with an 80 to 85 deg 
half-angle cone of visibility, corresponding to a 10 or 5 deg 

•Jf 

mask angle above the horizon. Or, it could be because the 
user has an antenna with a conical reception pattern, as, for 
example, a hemispherical antenna. In any event, if the user 
is restricted to a cone of satellite visibility of half angle 
109.5° or less, then the best possible four satellite geometry 
is the single parameter geometry shown in Fig. B.6-1, and Fig. 
B.6-3(b) shows the best obtainable TDOP, PDOP, and GDOP. The 
values quoted earlier for the half-cone angle of 90° are thus 
lower bounds for users on the surface of the earth, or with 
hemispherical antennas, making measurements to four satellites. 


*This cone could be larger if the user were high above the 
earth, such as for the Space Shuttle. 
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